# Quantum-Chemical Study of the FeNCN Conversion-Reaction Mechanism in Lithium- and Sodium-Ion Batteries

Kaixuan Chen, Marcus Fehse, Angelica Laurita, Jeethu Jiju Arayamparambil, Moulay Tahar Sougrati, Lorenzo Stievano, Richard Dronskowski

► To cite this version:

Kaixuan Chen, Marcus Fehse, Angelica Laurita, Jeethu Jiju Arayamparambil, Moulay Tahar Sougrati, et al.. Quantum-Chemical Study of the FeNCN Conversion-Reaction Mechanism in Lithium- and Sodium-Ion Batteries. Angewandte Chemie International Edition, 2020, 59 (9), pp.3718-3723. 10.1002/anie.201914760. hal-02514507

HAL Id: hal-02514507

https://hal.science/hal-02514507

Submitted on 31 May 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Quantum-Chemical Study of the FeNCN Conversion-Reaction Mechanism in Lithium- and Sodium-Ion Batteries

Kaixuan Chen, Marcus Fehse, Angelica Laurita, Jeethu Jiju Arayamparambil, Moulay Tahar Sougrati, Lorenzo Stievano,* and Richard Dronskowski*

**Abstract:** We report a computational study on 3d transition-metal (Cr, Mn, Fe, and Co) carbodiimides in Li- and Na-ion batteries. The obtained cell voltages semi-quantitatively fit the experiments, highlighting the practicality of PBE + U as an approach for modeling the conversion-reaction mechanism of the FeNCN arch type with lithium and sodium. Also, the calculated voltage profiles agree satisfactorily with experiment both for full (Li-ion battery) and partial (Na-ion battery) discharge, even though experimental atomistic knowledge is missing up to now. Moreover, we rationalize the structural preference of intermediate ternaries and their characteristic lowering in the voltage profile using chemical-bonding and Mulliken-charge analysis. The formation of such ternary intermediates for the lithiation of FeNCN and the contribution of at least one ternary intermediate is also confirmed experimentally. This theoretical approach, aided by experimental findings, supports the atomistic exploration of electrode materials governed by conversion reactions.

## Introduction

With the fast advancement in science and technology, and even faster growing demands by society, the field of energy production and storage is becoming critically important. Regarding high-density energy-generating approaches, there has been continuous progress in advanced nuclear power such as breeder/burner technology, thorium reactors, generation IV, and even small modular fission concepts.[1] With respect to low-density energy production and storage, extensive attention has been directed at research concerning solar

cells,[2] supercapacitors,[3] thermoelectrics,[4] Li- and Na-ion batteries (LIBs/NIBs),[5] as well as related fields. For reasons of convenience and the perception of the general public,[6] LIBs are widely used in various portable devices and electric vehicles,[7] but this only works as long as cheap primary electric energy is produced to feed them. Technically speaking, the performance of a Li/Na-ion battery, determined, for example, by its cell voltage, capacity, and reversible cycling number, essentially depends on the assembly of the electrodes and electrolytes.[8]

When viewing the frontier research of these electrodes from a somewhat wider perspective, two kinds of mechanisms seem to dominate, namely intercalation and the conversion reaction.[5b] Typical commercial electrodes such as LiCoO₂ and graphite follow the intercalation mechanism in which Li ions intercalate into or de-intercalate from the electrode host without significantly altering its original structure. Hence, the battery capacity is essentially restricted by the available vacancies in an electrode host to accommodate foreign Li ions. In the conversion-reaction mechanism, however, Li ions occupy not only the vacancies but also the sites of the metal cations, accompanied by a deposition of elemental metal particles. Because of this, drastic structural transformations may occur upon charge/discharge cycling, even leading to the forming of new phases. Compounds such as transition-metal oxides, nitrides, sulfides, and fluorides are generally governed by a conversion reaction.[5b,9]

In recent years, the family of transition-metal carbodiimides (TMC)[5a,10] has attracted increasing interest. The TMC family adopts the generic formula Mₓ(NCN)ᵧ, where M

[*] Dr. K. Chen, Prof. Dr. R. Dronskowski
Chair of Solid-State and Quantum Chemistry
Institute of Inorganic Chemistry, RWTH Aachen University
52056 Aachen (Germany)
E-mail: drons@HAL9000.ac.rwth-aachen.de

Dr. M. Fehse, A. Laurita, Dr. J. J. Arayamparambil, Dr. M. T. Sougrati, Prof. Dr. L. Stievano
Institut Charles Gerhardt Montpellier, CNRS
Université de Montpellier
34095 Montpellier (France)
E-mail: lorenzo.stievano@umontpellier.fr

Dr. M. Fehse, Dr. J. J. Arayamparambil, Dr. M. T. Sougrati, Prof. Dr. L. Stievano
Alistore—European Research Institute, CNRS
80039 Amiens (France)
Dr. M. Fehse
Dutch-Belgian (DUBBLE), ESRF—The European Synchrotron
38043 Grenoble (France)

Dr. M. T. Sougrati, Prof. Dr. L. Stievano
Reseau sur le Stockage Electrochimique de l'Energie (RS2E), CNRS
80039 Amiens (France)

Prof. Dr. R. Dronskowski
Hoffmann Institute of Advanced Materials, Shenzhen Polytechnic
7098 Liuxian Blvd, Nanshan District, Shenzhen (China)

<figure><img src="image_2.png" /></figure>

Supporting information and the ORCID identification number(s) for the author(s) of this article can be found under:
https://doi.org/10.1002/anie.201914760.

<figure><img src="image_3.png" /></figure>

© 2019 The Authors. Published by Wiley-VCH Verlag GmbH & Co. KGaA. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

<figure><img src="image_4.png" /></figure>
denotes a transition-metal cation and the $NCN^{2-}$ carbodiimide is the “divalent nitride” anion, mimicking the $O^{2-}$ anion by a nitrogen functionality. Its members, including $Cr_2(NCN)_3$, [11] $MnNCN$, [12] $FeNCN$, [13] $CoNCN$, [14] and $NiNCN$, [15] have been successfully fabricated by “soft chemistry” (metathesis) routes and furthermore, they turned out to be excellent anode candidates possessing a high reversible capacity upon charge/discharge cycling in Li- and Na-ion batteries.[5,16] Nevertheless, the detailed reaction pathway and unsolved structural transformations that occur upon charge/discharge cycling are difficult to detect experimentally because of the often amorphous nature of these materials. Hence, the main focus of this contribution is to carry out an in-depth density functional theory (DFT) study of the conversion-reaction mechanism of typical FeNCN as a LIB (NIB) electrode material by employing advanced DFT calculations combined with operando X-ray absorption spectroscopy (XAS).

## Results and Discussion

We calculated the average cell voltages[17] of the TMC (TM = Cr, Mn, Fe, Co) vs. pure Li and Na metals using the PBE and PBE + U functionals (see the Supporting Information).[18] Pure PBE produces less reliable voltage data (the largest error being 1.2 V for $CoNCN$ in a LIB) because the U correction is necessary to model the highly correlated TM electrons.[19] As shown in Table S2 (Supporting Information), introducing the U correction[20] effectively decreases the voltage errors (to about 0.8 V for $CoNCN$ in a LIB), so we chose PBE + U for a large-scale study.

Since we focus on the reaction pathways and voltage profiles during the charge/discharge process for the FeNCN phenotype in Li/Na-ion batteries, the crystal structures of the fully delithiated (desodiated) compound FeNCN and the lithiated (sodiated) compound $Li_2NCN$ ($Na_2NCN$) should be analyzed first (Figure 1). FeNCN crystallizes in the hexagonal system with $P6_3/mmc$ symmetry, that is, octahedral coordination for Fe and trigonal-prismatic coordination for NCN, resembling the nickel arsenide motif. $Li_2NCN$ ($Na_2NCN$), however, adopts a tetragonal $I4/mmm$ structure with tetrahedral coordination for Li (Na) and a quadratic-prismatic coordination for NCN. Like the related oxide compounds, transition-metal carbodiimides follow a conversion-reaction mechanism when used as electrode materials. That means that

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Crystal structures of a) FeNCN and b) $Li_2NCN/Na_2NCN$.</figcaption></figure>

drastic structural transformations occur during the charge/discharge process. It is unclear up to which degree the crystal structures of the intermediate compounds may resemble FeNCN or $Li_2NCN$ in a LIB. It is therefore far from trivial to correctly deduce the atomistic structural details of these ternaries denoted as Li-Fe-NCN, in particular because reliable experimental models (for example, from X-ray diffraction) are missing. Because of this, we tentatively used four kinds of structural bases for a thorough survey of possible Li-Fe-NCN ternaries: two from divalent carbodiimides,[21] namely the NiAs-type FeNCN and the NaCl-type $MnNCN$, one from $Li_2NCN$, and one from Li-containing metal oxides such as $CuLi_2O_2$ (where one needs to replace $O^{2-}$ anions with $NCN^{2-}$ groups).

During the discharge process in a Li-ion half-battery with Li metal as the anode, the cathode (FeNCN) accommodates both $Li^+$ cations (through the electrolyte) and electrons (through an external loop) to form Li-Fe-NCN ternaries. Note that an oxidation reaction in the positive FeNCN electrode forming $Fe^{3+}$ cations is prevented during lithiation (cathodic conditions), which rather leads to the reduction of the material producing the deposition of elemental $Fe^0$. To maintain electrical neutrality in a Li-Fe-NCN ternary, $Fe^{2+}/Fe^0$ reduction and $Li^0/Li^+$ oxidation must proceed in a 1:2 proportion since one $Fe^{2+}$ acts as two $Li^+$ in the charge compensation. In other words, a cathode host accommodates two Li and precipitates one Fe simultaneously. That is, the intermediate compounds, whatever they may be, share the generic chemical formula $Li_xFe_{(1-x/2)}NCN$ with $0 < x < 2$. We enforced this confinement due to charge balance when establishing the model-type intermediate compounds[22] as derived from the aforementioned four structural bases and as tabulated in Table 1.

<table><caption>**Table 1:** Total number of computations to survey the possible Li-Fe-NCN (Na-Fe-NCN) ternary compounds.</caption><thead><tr><th rowspan="2">structural basis</th><th rowspan="2">space group</th><th colspan="2">no. of calculated compounds</th></tr><tr><th>Li-ion</th><th>Na-ion</th></tr></thead><tbody><tr><td>FeNCN</td><td>P6<sub>3</sub>/mmc (no. 194)</td><td>17</td><td>17</td></tr><tr><td>Li<sub>2</sub>NCN/$Na_2NCN$</td><td>I4/mmm (no. 139)</td><td>23</td><td>16</td></tr><tr><td>MnNCN</td><td>R3m (no. 166)</td><td>10</td><td>10</td></tr><tr><td>CuLi<sub>2</sub>O<sub>2</sub></td><td>Immm (no. 71)</td><td>11</td><td>11</td></tr></tbody></table>

After full optimizations, the structures were used to perform electronic self-consistency calculations. To examine the stability of each intermediate compound, its formation energy $E_{formation}$ was calculated, which is defined as the energetic difference between an intermediate compound $Li_xFe_{(1-x/2)}NCN$ and its corresponding linear combination of $Li_2NCN$ (fully lithiated) and FeNCN (fully delithiated):

$$E_{formation} = E(\mathrm{Li}_x\mathrm{Fe}_{(1-x/2)}\mathrm{NCN}) - \frac{x}{2}E(\mathrm{Li}_2\mathrm{NCN}) - (1 - \frac{x}{2})E(\mathrm{FeNCN}) \quad (1)$$

Figure 2a shows the formation energies of all calculated $Li_xFe_{(1-x/2)}NCN$ compounds at various configurations. By convention, the more negative the formation energies, the
<figure><img src="image_6.png" /><figcaption>Figure 2. a) Formation energies (crosses) of the calculated Li-Fe-NCN ternary compounds together with the energy convex hull as a blue line and b) voltage profile of FeNCN upon lithiation in a LIB.</figcaption></figure>

more stable the compounds. All unstable (or metastable) compounds with positive formation energies were not included in Figure 2a. For convenience, a line connecting the compounds that are most stable at each configuration is also depicted in blue (minimum-energy curve or energy convex hull).[23] Given full thermodynamic control, the reaction pathway during the charge/discharge process would have to follow this line (related phase diagram shown in Figure S3).

Each voltage pair (given by two intermediate compounds with adjacent Li content in the phase diagram) gives a voltage plateau, and a voltage profile can be obtained by combining all voltage plateaus. For example, given a voltage pair of two intermediate compounds written as $Li_{x_1}Fe_{(1-x_{1/2})}NCN$ and $Li_{x_2}Fe_{(1-x_{2/2})}NCN$ (taking $x_1 < x_2$ with FeNCN as a cathode host in a LIB), their reaction reads:

$$ (x_2 - x_1)Li + Li_{x_1}Fe_{(1-x_{1/2})}NCN \rightarrow Li_{x_2}Fe_{(1-x_{2/2})}NCN + \frac{1}{2}(x_2 - x_1)Fe, \quad (2) $$

so that the voltage plateau $\Phi(x_1, x_2)$ results, similar to that from $M_x(NCN)_y$ to $Li_2NCN$, in:

$$ \Phi(x_1, x_2) = -\frac{\Delta G}{(x_2 - x_1)zF} \approx -[E(Li_{x_2}Fe_{(1-x_{2/2})}NCN) + \frac{1}{2}(x_2 - x_1)E(Fe) - E(Li_{x_1}Fe_{(1-x_{1/2})}NCN) - (x_2 - x_1)E(Li)][(x_2 - x_1)zF]^{-1} \quad (3) $$

Figure 2b depicts the theoretical voltage profile for the archetyp FeNCN using a blue solid line, corresponding to the minimum-energy curve in Figure 2a. The staircase shape originates from the limited number of stable compounds obtained from the also limited number of DFT calculations.

For conversion-reaction mechanism materials, there is clearly a large experimental difference between the first and second lithiation, but the subsequent cycles are identical.[24] Indeed, during the first lithiation process, large and well-crystallized particles are transformed into nanosized and amorphous particles. Only after this irreversible morphological change, the lithiated material becomes the true cycling material. Thus, the experimental profile of FeNCN at the second cycle at a rate of C/50 (corresponding to the reaction of one mole of Li per mole of FeNCN in 50 h) is depicted (as a dashed line) for comparison.

In the theoretical data, an abnormally high voltage plateau (about 4.2 V vs. $Li^+/Li^0$) is observed at a low Li content ($x < 0.2$). This effect is a computational artifact since an extremely large supercell would be required to perfectly distribute Li sites and Fe vacancies at very low Li content, which cannot be modeled with our limited computational resources. Fortunately, the requirement is much weaker for a medium Li content, which makes the voltage profile more reliable in that range. Another divergence is the maximally available Li content $x$, which should be 2 (theoretically) when the fully lithiated compound $Li_2NCN$ forms. In the experiment, $x$ can become larger than 2, and actually goes up to about 2.8 because of side reactions. We note that the electrode is made not only of FeNCN but also includes 40 w% of carbon and binder (CMC, carboxymethyl cellulose). Therefore, the reasons for the larger Li content may include the formation of a passivation layer (SEI, solid-electrolyte interface) on the active material and electrolyte degradation. This is also the reason why the capacity of the first lithiation is always higher than that of subsequent cycles. For FeNCN, the SEI (or electrolyte degradation) seems to be largely reversible, and seemingly keeps cycling at more than 2Li. This phenomenon is also observed in transition-metal oxides.[14] Despite the two unavoidable divergences, the calculated voltage profile agrees well with experiment and nicely reproduces a similar voltage plateau for a wide Li range ($0.2 < x < 2.0$).

Figure 3a,b depicts the minimum-energy curve and voltage profile for the FeNCN archetyp in a NIB, obtained with the same computational approach as used for the LIB before. An analogous high voltage plateau (about 3.4 V vs. $Na^+/Na^0$) appears at a low Na content ($x < 0.2$), similar to the LIB. The maximum experimentally achievable Na content, however, is only $x = 1.5$. This is because Na fails to react with all FeNCN in experiments, probably due to kinetic hindrance, so there is only a partial reaction (or partial discharge). The SEI (or electrolyte degradation) effect is the same as in a LIB, but with less reversible capacity to form the active material. If we set $x = 1.5$ as the theoretical maximum Na content, that is, the discharge process ends at $Na_6Fe(NCN)_4$ and not $Na_2NCN$, we end up with a new voltage profile for a partial discharge. It is reassuring that this new profile, depicted as a solid orange line in Figure 3b, agrees nicely with the experimental result.
<figure><img src="image_7.png" /><figcaption>**Figure 3.** a) Formation energies (crosses) of the calculated Na-Fe-NCN ternary compounds together with the energy convex hull as a yellow line and b) voltage profile of FeNCN upon lithiation in a NIB.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 4.** a) Chemical-bonding analysis of the Fe-N and Li-N bonds in three $Li_2Fe_7(NCN)_8$ structural models. The Fermi level is set to zero.
b) Mulliken atomic charges of Fe and ICOHP data of the Fe-N bonds in the energetically favorable Na-Fe-NCN ternaries.</figcaption></figure>

For further interpretation and understanding of the aforementioned findings, chemical-bonding analysis was employed. In the cathode of a LIB or NIB, the ideal reaction pathway (or minimum-energy curve) starts at FeNCN, proceeds through a variety of intermediate (amorphous) compounds, and ends at $Li_2NCN$ or $Na_2NCN$. For an intermediate compound with a certain stoichiometric configuration, for example, $Li_2Fe_7(NCN)_8$, different crystal-structure models were computationally established. Taking a thermodynamic point of view, the one crystallizing with the lowest formation energy was preferred but we note that other metastable variants with higher formation energies but lower activation barriers are also possible according to Ostwald's step rule.[25] In the following, we rationalize the occurrence of structural polymorphs from the aspect of chemical bonding. Figure 4a depicts the composition $Li_2Fe_7(NCN)_8$ (corresponding to a Li content of $x = 0.25$) for which three structural models ($P6_3/mmc$, Immm, and I4/mmm) were established. Among these, I4/mmm has the lowest formation energy. Previous research indicates that the structural preference in carbodiimides can be determined from chemical-bonding analysis.[21,26] Figure 4a presents a crystal orbital Hamilton population (COHP) analysis of both Fe-N and Li-N bonds in the three $Li_2Fe_7(NCN)_8$ models adopting $P6_3/mmc$, Immm, and I4/mmm symmetry. Negative (or positive) COHP values indicate energy-lowering bonding (or energy-increasing anti-

bonding) interactions. It is apparent that bonding interactions dominate in the low-energy valence-band region whereas anti-bonding states appear in the unoccupied conduction band, but also just below the Fermi level. For a quantitative assessment, all COHP values were integrated up to the Fermi level to yield ICOHP data for the Fe-N and Li-N bonds and averaged over all corresponding bonds in the unit cell. It turns out that the I4/mmm model of $Li_2Fe_7(NCN)_8$ not only holds the strongest bonding interactions for both Fe-N and Li-N bonds but also has the lowest formation energy. Hence, the chemical-bonding analysis is useful to deduce where the structural preference of the intermediate compounds originates from.

As discussed above, the cell voltage for FeNCN originates from the potential difference between the $Fe^{2+}/Fe^0$ and $Li^+/Li^0$ (Na$^+/Na^0$) redox couples in a LIB (NIB). An ideal voltage profile displays a flat plateau with a theoretical voltage of 1.6 V (1.1 V) in a LIB (NIB), as tabulated in Table S2. The actual profile, however, exhibits a monotonous decrease in voltage upon lithiation (sodiation) both experimentally and theoretically. We rationalize this decreasing character from a monotonous decrease in the atomic charge of Fe, expressed as $y$ in $Fe^{y+}$ in the energetically low-lying intermediate ternaries, and evaluated from Mulliken charge analysis in Figure 4b. To zeroth order, the higher the cationic charge, the stronger the reduction power. Indeed, the Mulliken charge of
Fe decreases from 1.1 in FeNCN to 0.9 in Na₂₂Fe₇(NCN)₃₈, seemingly leading to a decrease in the Fe⁺/Fe⁰ redox potential; note that these (basis-set-independent) Mulliken charges were numerically projected from the entire wave-function expanded from plane waves.[27] Figure 4b also demonstrates a more negative ICOHP progression for the Fe–N bonds, but now upon sodiation in a NIB. This indicates an almost continuous increase in covalency for the Fe–N bonds in the energetically low-lying intermediate ternaries, another confirmation of a smaller amount of electron density transferred from Fe to N and, thus, a decrease in the charge of Fe.

At last, the lithiation of FeNCN was followed experimentally by operando Fe K-edge XAS in order to corroborate or falsify the theoretical findings.[28] In particular, the extended fine structure of the XAS spectrum (EXAFS)[29] is very sensitive to the modification of the nearest-neighbor shell of the Fe atoms and is hence used to validate the DFT-predicted ternary intermediates. The complete operando data set depicting the evolution of Fe K-edge XAS spectra acquired during the first electrochemical discharge of FeNCN vs. Li is presented as a contour plot in Figure 5a. The plot shows a decrease of the edge-energy intensity accompanied by a shift to lower edge energy, reflecting the effective chemical reduction of the Fe atoms.

<figure><img src="image_9.png" /><figcaption>**Figure 5.** a) Contour plot of the entire Fe K-edge XAS operando data set revealing the evolution of the XANES and EXAFS of the Fe K-edge upon electrochemical discharge vs. Li; b) Fourier transform of the EXAFS of the Fe K-edge of intermediate component 2 obtained via MCR-ALS and EXAFS fit with a fitting window.</figcaption></figure>

Principal component analysis (PCA)[30] indicates that three independent principal components are needed to reproduce the entire data set, yielding 99.9% total variance,

as shown in Figure S4. Subsequently pure spectral components were reconstructed via multivariate curve resolution alternating least squares (MCR-ALS).[31] The EXAFS spectra indicate that component 1 and component 3 can be easily fitted with the structure of pristine FeNCN and the iron nanoparticles, as shown in Figure S6 and S7, respectively. Component 2 represents the ternary intermediate, which reaches its maximum concentration at about one third of the first discharge (Figure S5b), and can therefore be successfully fitted by the theoretical structure of Li₂Fe₂-(NCN)₃, in line with the amount of reacted lithium. The compound Li₂Fe₂(NCN)₃, structurally derived from Li₂NCN, yields satisfactory results, as shown in Figure 5b and Table 2.

<table><caption>**Table 2:** EXAFS fitting parameters of MCR-ALS component 2 of the Fe K-Edge.</caption><thead><tr><th>shell</th><th>coord. no.</th><th>theor. distance (Å)</th><th>exp. distance (Å)</th><th>σ<sup>2</sup></th></tr></thead><tbody><tr><td>Fe–N</td><td>4</td><td>2.07</td><td>2.14(1)</td><td>0.007(2)</td></tr><tr><td>Fe–C</td><td>3</td><td>2.93</td><td>3.06(5)</td><td>0.01(1)</td></tr><tr><td>Fe–Fe</td><td>2</td><td>3.46</td><td>3.27(1)</td><td>0.007(3)</td></tr></tbody></table>

The quantum-chemical bonding analysis in Figure S8 also indicates that the Fe–N bonds in the Li₂Fe₂(NCN)₃ structure derived from the Li₂NCN supercell possess a more negative ICOHP value (meaning stronger bonding, –1.84 eV) than the ones in the FeNCN supercell (–1.46 eV). In other words, the good agreement of the DFT-predicted structure with the experimental spectra confirms the formation of an intermediate yet amorphous phase compatible with the local structure of Li₂Fe₂(NCN)₃ during the conversion reaction of FeNCN vs. Li. Unravelling such transient phases in complex reaction mechanisms is a valuable asset of the MCR-ALS approach, was recently shown for sodiation of SnSb.[32]

## Conclusion

In conclusion, given the new family of 3d transition-metal carbodiimides which have experimentally turned out as excellent candidates for Li- and Na-ion batteries, we have demonstrated the first computational attempt to model its electrochemical mechanism. A first-principles study on the conversion-reaction mechanism of the archetyp e FeNCN with Li and Na during charge/discharge has been performed by simulating a total of 61 (54) configurations of Li–Fe–NCN (Na–Fe–NCN) ternaries. Based on their structures, spin-polarized PBE + U calculations gave access to voltage profiles which agree nicely with experiment for both full (LIB) and partial (NIB) discharge. The origins of the energetically preferred structures adopted in the intermediate compounds were rationalized by chemical-bonding analysis. Additionally, the continuously lowering voltage profile was also mirrored in the COHP progression and the Mulliken charge transfer. Operando XAS analysis performed during the lithiation of FeNCN substantiates the results of these calculations, confirming the formation of an intermediate species with a structure compatible with the theoretically
predicted transient phase Li₂Fe₂(NCN)₃. Finally, this study offers a new challenge to chemists in the preparation of metastable Li-Fe-NCN intermediates.

## Acknowledgements

K.C. gratefully acknowledges financial support from the Alexander von Humboldt Foundation. The simulation work was supported by the IT center of RWTH Aachen University under the JARA-HPC grant JARA0179. We thank the NWO and beamline DUBBLE at ESRF, Grenoble (France) for providing beamtime, and in particular Dr. Alessandro Longo for technical advice on beamline operation. This work was supported by the public grant ANR-10-LABX-76-01 (Labex STOREX) overseen by the French National Research Agency (ANR) as part of the "Investissements d'Avenir" program. Alistore-ERI is warmly thanked for financing the PhD grant of J.J.A. and the post-doc grant of M.F.

## Conflict of interest

The authors declare no conflict of interest.

**Keywords:** batteries · carbodiimides · computational chemistry · density functional calculations · lithium ions

**How to cite:** *Angew. Chem. Int. Ed.* **2020**, *59*, 3718–3723
*Angew. Chem.* **2020**, *132*, 3747–3752

<footnote>[1] a) D. Bodansky, *Nuclear Energy: Principles, Practices, and Prospects*, Springer, New York, 2010; b) P. Breeze, *Nuclear Power*, Elsevier, Dordrecht, 2017.</footnote>

<footnote>[2] P. K. Nayak, S. Mahesh, H. J. Snaith, D. Cahen, *Nat. Rev. Mater.* **2019**, *4*, 269–285.</footnote>

<footnote>[3] G. Z. Chen, *Int. Mater. Rev.* **2017**, *62*, 173–202.</footnote>

<footnote>[4] a) K. X. Chen, S. H. Lyu, X. M. Wang, Y. X. Fu, Y. Heng, D. C. Mo, *J. Phys. Chem. C* **2017**, *121*, 13035–13042; b) K. X. Chen, X. M. Wang, D. C. Mo, S. S. Lyu, *J. Phys. Chem. C* **2015**, *119*, 26706–26711; c) K. X. Chen, M. S. Li, D. C. Mo, S. S. Lyu, *Front. Energy* **2018**, *12*, 97–108.</footnote>

<footnote>[5] a) M. T. Sougrati, A. Darwiche, X. Liu, A. Mahmoud, R. P. Hermann, S. Jouen, L. Monconduit, R. Dronskowski, L. Stievano, *Angew. Chem. Int. Ed.* **2016**, *55*, 5090–5095; *Angew. Chem.* **2016**, *128*, 5174–5179; b) J. Cabana, L. Monconduit, D. Larcher, M. R. Palacín, *Adv. Mater.* **2010**, *22*, E170–192.</footnote>

<footnote>[6] The Royal Swedish Academy of Sciences (2019, Oct 9). *The Nobel Prize in Chemistry 2019*. Retrieved from https://www.nobelprize.org/prizes/chemistry/2019/press-release.</footnote>

<footnote>[7] D. Lin, Y. Liu, Y. Cui, *Nat. Nanotechnol.* **2017**, *12*, 194–206.</footnote>

<footnote>[8] a) J. B. Goodenough, K. S. Park, *J. Am. Chem. Soc.* **2013**, *135*, 1167–1176; b) M. V. Reddy, G. V. Subba Rao, B. V. Chowdari, *Chem. Rev.* **2013**, *113*, 5364–5457.</footnote>

<footnote>[9] Y. Zhao, L. P. Wang, M. T. Sougrati, Z. X. Feng, Y. Leconte, A. Fisher, M. Srinivasan, Z. C. Xu, *Adv. Energy Mater.* **2017**, *7*, 1601424.</footnote>

<footnote>[10] a) T. D. Boyko, R. J. Green, R. Dronskowski, A. Moewes, *J. Phys. Chem. C* **2013**, *117*, 12754–12761; b) C. Liu, C. Zhang, H. Fu, X. Nan, G. Cao, *Adv. Energy Mater.* **2017**, *7*, 1601127.</footnote>

<footnote>[11] X. Tang, H. Xiang, X. Liu, M. Speldrich, R. Dronskowski, *Angew. Chem. Int. Ed.* **2010**, *49*, 4738–4742; *Angew. Chem.* **2010**, *122*, 4846–4850.</footnote>

<footnote>[12] X. Liu, M. Krott, P. Müller, C. Hu, H. Lueken, R. Dronskowski, *Inorg. Chem.* **2005**, *44*, 3001–3003.</footnote>

<footnote>[13] X. Liu, L. Stork, M. Speldrich, H. Lueken, R. Dronskowski, *Chem. Eur. J.* **2009**, *15*, 1558–1561.</footnote>

<footnote>[14] J. J. Arayamparambil, M. Mann, B. Fraisse, A. Iadecola, R. Dronskowski, L. Stievano, M. T. Sougrati, *ChemElectroChem* **2019**, *6*, 5101–5108.</footnote>

<footnote>[15] M. Krott, X. Liu, B. P. Fokwa, M. Speldrich, H. Lueken, R. Dronskowski, *Inorg. Chem.* **2007**, *46*, 2204–2207.</footnote>

<footnote>[16] M. T. Sougrati, J. J. Arayamparambil, X. Liu, M. Mann, A. Slabon, L. Stievano, R. Dronskowski, *Dalton Trans.* **2018**, *47*, 10827–10832.</footnote>

<footnote>[17] A. Urban, D.-H. Seo, G. Ceder, *NPJ Comput. Mater.* **2016**, *2*, 1–13.</footnote>

<footnote>[18] G. Kresse, J. Furthmüller, *Phys. Rev. B* **1996**, *54*, 11169–11186.</footnote>

<footnote>[19] Y. S. Meng, M. E. Arroyo-de Dompablo, *Energy Environ. Sci.* **2009**, *2*, 589–609.</footnote>

<footnote>[20] A. Jain, G. Hautier, S. P. Ong, C. J. Moore, C. C. Fischer, K. A. Persson, G. Ceder, *Phys. Rev. B* **2011**, *84*, 045115.</footnote>

<footnote>[21] K. Chen, R. Dronskowski, *J. Phys. Chem. A* **2019**, *123*, 9328–9335.</footnote>

<footnote>[22] S. P. Ong, W. D. Richards, A. Jain, G. Hautier, M. Kocher, S. Cholia, D. Gunter, V. L. Chevrier, K. A. Persson, G. Ceder, *Comput. Mater. Sci.* **2013**, *68*, 314–319.</footnote>

<footnote>[23] H. Kim, D. J. Kim, D.-H. Seo, M. S. Yeom, K. Kang, D. K. Kim, Y. Jung, *Chem. Mater.* **2012**, *24*, 1205–1211.</footnote>

<footnote>[24] M. Fehse, A. Darwiche, M. T. Sougrati, E. M. Kelder, A. V. Chadwick, M. Alfredsson, L. Monconduit, L. Stievano, *Chem. Mater.* **2017**, *29*, 10446–10454.</footnote>

<footnote>[25] R. A. Van Santen, *J. Phys. Chem.* **1984**, *88*, 5768–5769.</footnote>

<footnote>[26] a) S. Maintz, V. L. Deringer, A. L. Tchougréeff, R. Dronskowski, *J. Comput. Chem.* **2016**, *37*, 1030–1035; b) W. L. Li, C. Ertural, D. Bogdanovski, J. Li, R. Dronskowski, *Inorg. Chem.* **2018**, *57*, 12999–13008; c) S. Maintz, V. L. Deringer, A. L. Tchougréeff, R. Dronskowski, *J. Comput. Chem.* **2013**, *34*, 2557–2567; d) R. Dronskowski, P. E. Blöchl, *J. Phys. Chem.* **1993**, *97*, 8617–8624; e) V. L. Deringer, A. L. Tchougréeff, R. Dronskowski, *J. Phys. Chem. A* **2011**, *115*, 5461–5466.</footnote>

<footnote>[27] C. Ertural, S. Steinberg, R. Dronskowski, *RSC Adv.* **2019**, *9*, 29821–29830.</footnote>

<footnote>[28] J. B. Leriche, S. Hamelet, J. Shu, M. Morcrette, C. Masquelier, G. Ouvrard, M. Zerrouki, P. Soudan, S. Belin, E. Elkaim, F. Baudelet, *J. Electrochem. Soc.* **2010**, *157*, A606–A610.</footnote>

<footnote>[29] M. Fehse, A. Iadecola, M. T. Sougrati, P. Conti, M. Giorgetti, L. Stievano, *Energy Storage Mater.* **2019**, *18*, 328–337.</footnote>

<footnote>[30] D. L. Massart, B. G. M. Vandeginste, L. M. C. Buydens, S. De Jong, P. J. Lewi, J. Smeyers-Verbeke, *Handbook of Chemo-metrics and Qualimetrics: Part A, Vol. 20*, Elsevier, Dordrecht, **1997**.</footnote>

<footnote>[31] a) J. Jaumot, R. Gargallo, A. de Juan, R. Tauler, *Chemom. Intell. Lab. Syst.* **2005**, *76*, 101–110; b) J. Jaumot, A. de Juan, R. Tauler, *Chemom. Intell. Lab. Syst.* **2015**, *140*, 1–12.</footnote>

<footnote>[32] M. Fehse, M. T. Sougrati, A. Darwiche, V. Gabaudan, C. La Fontaine, L. Monconduit, L. Stievano, *J. Mater. Chem. A* **2018**, *6*, 8724–8734.</footnote>

<footnote>Manuscript received: November 19, 2019</footnote>

<footnote>Accepted manuscript online: December 12, 2019</footnote>

<footnote>Version of record online: January 22, 2020</footnote>