Homogeneity of lithium distribution in cylinder-type Li-ion batteries

Spatially-resolved neutron powder diffraction with a gauge volume of 2 × 2 × 20 mm3 has been applied as an in situ method to probe the lithium concentration in the graphite anode of different Li-ion cells of 18650-type in charged state. Structural studies performed in combination with electrochemical measurements and X-ray computed tomography under real cell operating conditions unambiguously revealed non-homogeneity of the lithium distribution in the graphite anode. Deviations from a homogeneous behaviour have been found in both radial and axial directions of 18650-type cells and were discussed in the frame of cell geometry and electrical connection of electrodes, which might play a crucial role in the homogeneity of the lithium distribution in the active materials within each electrode.

Li-ion batteries are currently dominating the field of electrochemical energy storage especially in portable electronic and electric vehicles applications. Their outstanding position is primarily due to their high power and energy densities, good cycle life and excellent storage characteristics, thus making Li-ion battery technology a favorite on the broad selection of different electrochemical storage media. However there is still space for optimization especially when safety, energy and power, price, weight or operating temperature are concerned. This caused active research on Li-ion batteries in recent years in line with the rapid development of portable electronics and the general trend substituting fossil fuel by renewable energy. As a consequence cell prices are rapidly falling and an annual increase of ca. 8.5% for the accessible cell capacity (from ca. 1.2 Ah for first 18650-type cells from Sony in 1991 to about 3.4 Ah in 2013) is achieved.
Strongly application-oriented lithium-ion batteries are now built in different cell designs with the aim to achieve highest cycle life, energy and/or power density. Besides the fact, that a modern Li-ion cell is a sophisticated electrochemical device possessing numerous and intricately coupled degrees of freedom, it may have a complicated shape with different ways of electrode layer layout, electrical connection, electrolyte wetting etc. Simple in principle but complicated in practice these cell design features may result in spatial inhomogeneities of current, lithium or electrolyte distribution, which are often difficult to quantify but will surely affect performance, cycling stability and safety 1 .
A typical lithium-ion battery is an electrically and environmentally isolated stand-alone system, whose thermodynamic or chemical state changes immediately after any interference. This creates numerous challenges for studies of Li-ion batteries and restricts the broad range of experimental and analytic techniques to the use of in situ and non-destructive characterization methods. Hereby risks of possible materials oxidation, electrolyte evaporation, battery state-of-charge (SOC) changes etc. are eliminated. Bearing in mind that a high penetration depth and a sensitivity to light atoms is needed, the applicable methods out of the accessible experimental toolbox become quite limited.
In such instance neutron scattering is already a well-established experimental technique for the characterization of complex Li-ion batteries. It has numerous advantages 2 : the high penetration depth of thermal neutrons permits non-destructive in situ and in operando studies of different kinds; the capability to localize light elements (e.g. hydrogen, lithium) and to resolve isotopes provides the necessary contrast; the neutron scattering length not dependent on momentum transfer, enabling to achieve accurate structure factors at high 2θ angles etc. Neutron scattering methods undergo an increasing relevance for studies of lithium-ion batteries on different length scales, e.g. neutron imaging 3-5 , reflectometry 6 , small-angle neutron scattering 7-10 , quasielastic neutron scattering 11 and powder diffraction 5,[12][13][14][15][16][17] . In situ experiments with neutrons are performed on self-developed test cells and commercial Li-ion cells of different designs depending on the research needs. When neutron powder diffraction is considered the adapted coin-type cells 18 , pouch-back cells 19 , Swagelok-based constructions 20,21 , cylinder winded cells [22][23][24] or commercial design of 18650 5 /prismatic 12 types are often utilized with the main focus on the development of atomic structure of electrode materials upon various electrochemical and environmental states.
Despite the increasing popularity of neutron scattering studies of batteries at their operating conditions (especially neutron powder diffraction where typically large cell cross sections are exposed to the neutron beam) the problem of cell homogeneity and its effect on the acquired data, e.g. unrealistic peak intensities, reflection shift or broadening is often not properly accounted in literature. However the fact that characteristic peaks at plateaus on differential capacity plots (cycled under similar conditions) for cells in 2032 standard (coin-type) are narrower than their counterparts from 18650-type cells indirectly points out a possible inhomogeneity issue in larger format Li-ion cells 25 . This can be associated with the development of different states of charge in large Li-ion cells due to the voltage drop along the current collector during charge 26 . Direct optical measurement of Li transport in an operating Li-ion cell 27 , indicated a development of lithium gradient through the graphite electrode with complicated microstructure 28 .
To our best knowledge structural studies of cell homogeneity by diffraction methods are limited to the work of Cai et al. 29 . They revealed an inhomogeneous deterioration at the edges of a pouch-back cell after active cycling. Recently 30 we reported results of spatially resolved monochromatic neutron powder diffraction applied to commercial Li-ion cells of 18650-type and LCO cathode along the radial axis where the uniformity of Li distribution in the cell volume was determined along the studied axis. However further experiments performed in other radial directions and various types of cells showed that this not always holds true and requires more systematic characterization.
In the current manuscript we report a further study of the lithium distribution in the graphite anode of different cells of 18650-type by using spatially resolved neutron powder diffraction. Four types of commercial cells based on graphite anodes and different cathode materials were chosen for the experiment. Their properties are listed in Table 1 including the applied voltage windows and nominal capacities.
The layout of cells 1-4 was non-destructively probed using X-ray tomography, where results of the tomography reconstruction for cell 1 are shown in Fig. 1a,b,d-f. Due to poor sensitivity of X-ray based methods to light atoms the lithium distribution in the anode was studied using spatially resolved neutron diffraction. By translating and rotating the sample with respect to the fixed gauge volume a spatially-resolved scattering response from different sample sections can be obtained (see Methods section and Fig. S1 for details). All investigated cells (Nr. 1-4) were studied in charged state. Scans at the center of the cells were performed in a base cylinder plane (Fig. 1a). The distribution of the studied gauge volumes (red rhombs) is shown in Fig. 1b. A diffraction pattern similar to that displayed in Fig. S1c was collected at each point. The lithium concentration in Li x C 6 was determined using relative intensities of 001 LiC 6 and 002 LiC 12 Bragg reflections. Surfaces of spatial x distribution were constructed by the interpolation of experimental data using biharmonic spline interpolation 31 (see Fig. S2 for details) and are shown for cell 1 in Fig. 1c using false color presentation with x ranging from 0.5 (blue) to 0.95 (red).
In order to derive the lithium concentrations the analysis of relative peak intensities has a number of advantages compared to a previously reported evaluation of absolute intensities of single peaks altered by neutron attenuation 29,30 , e.g. the effect of attenuation can be neglected in first approximation, lithium inhomogeneities will become apparent in form of deviations from the flat surface and no signal is expected when no material is in the beam (within the studied gauge volume). From Fig. 1 one can see that for cell 1 the 2D map of the lithium distribution shows homogeneous character besides the narrow regions at the center and close-to the surface. For the sake of comparison two lithium concentrations inside the anode will be analyzed, namely x P given by the plateau with homogeneous lithium distribution and the mean lithium content < x> obtained after averaging over all investigated gauge volumes. For cell 1 a fair agreement between x P = 0.88(1) and < x> = 0.85(4) was observed (see Table 1). One has to ensure that in areas with reduced lithium concentrations the couple of 001 LiC 6 and 002 LiC 12 reflections are clearly observed (Fig. S2) indicating the presence of both stage I and stage II lithiated carbons in the corresponding gauge volume. Otherwise the determination of the lithium concentration x from the ratio of the two peak intensities is not possible.
The existence of regions with low lithium concentrations can be explained in the frame of a typical design of cylinder cells adopting double-coated electrodes wound around a center pin. A reconstructed slice form the X-ray tomography corresponding to the diffraction plane is shown in Fig. 1d-f. In case of cell 1 the typical sequence of the electrode layers is "…/Cu/C| S| LiCoO 2 /Al/LiCoO 2 | S| C/Cu/C| S| LiCoO 2 /Al/…", where S is the separator, LiCoO 2 / Al/LiCoO 2 corresponds to the cathode coated on both sides of aluminium foil and C/Cu/C is the graphite anode on copper foil, respectively. Due to the necessary connection to current leads the electrode foils are not coated at their ends (Fig. 1d,f). The coating of opposite electrodes (LiCoO 2 /Al/LiCoO 2 and C/Cu/C) usually does not start at the same winding position. This may result in layer sequences like "… S| C/Cu/C| S| C/Cu/C| S …" or "… S| Al| S| C/Cu/C| S …" where there is no cathode in the neighbouring layer and, correspondingly, no lithiated graphite. A similar feature may be observed on the cathode sheet. It has to be noted that no traces of the graphite 002 peak (theoretical position 27.78 deg. 2θ ) were observed in the current study. This along with the structural evidence for lithiated carbons can be associated with sufficiently high intralayer lithium diffusion along the negative electrode or a sufficient distribution and high ionic conductivity of the electrolyte. The observed lithium deficiency in cell 1 at regions close to current leads connections, the center pin and the outermost layer of the cell may be explained by a lack of cathode material, which is only partially compensated by the ionic conductivity of the electrolyte and inside the electrode. Similar experiments were performed on other representatives of 18650-type cells. Cell 2 has been found to possess a similar design to that of cell 1 (as confirmed by X-ray computed tomography) and delivered similar spatial distribution of lithium in the graphite anode (Fig. 2a). A lower grade of graphite lithiation was observed in cell 2, which can be attributed to different cell balancing having a higher excess of graphite. A flat concentration plateau with x P = 0.82(1) was observed in a (narrow) range of the cell with radii between 3 mm and 7 mm around the centre pin (Fig. S3). At the outer part of the cell a lower lithium concentration in the anode can be identified (Fig. S6a), which causes a higher discrepancy between x P and < x> = 0.77(3) and accordingly indicates a less homogeneous lithium distribution in cell 2 compared to cell 1.
Cell 3 was characterized by the highest capacity in the series of studied cells. From Figs 2b, S4 and S6c it can be seen that the plateau value x p = 0.93(1) is the highest in the series, which indicates a well chosen cell balancing. Cell 3 also possesses a difference in design (Fig. 2b) -its positive current lead (connecting to the cathode) is located in the middle of the electrode and therefore positioned between the layers of the rolled stack of electrodes. This geometrical detail has a pronounced effect on the homogeneity of the lithium distribution in the graphite anode. At the location of current leads the lithiation level has been found definitely lower: 0.58(1) at the negative (anode) current lead (outer part) and ca. 0.77(1) at the positive (cathode) current lead located in between the rolled electrodes at a radius of approximately 6 mm from the center of cell 3. (Note: The lithium concentration in the anode has been found slightly perturbed at the opposite side (red plot in Fig. S6c) of the cell at ca. 6 mm radial distance from the center of cell 3 (marked by a blue arrow in Fig. S6c). This can be attributed to an inhomogeneous rolling due to a thickness misfit of the current lead and the electrode layers.) This results in a high mean lithium concentration in the anode of cell 3 with < x> = 0.87(5) and in a stronger inhomogeneity compared to cells 1 and 2.
Even more pronounced deviations from a "flat" lithium distribution were noticed for cell 4 (Fig. 2c, S5, S6d) possessing the lowest cell capacity among all investigated cells. A flat "plateau" with relatively low x P = 0.77(1) occurs in a narrow range of radial distances from 4 mm to 7 mm from the centre pin, besides which the lithium concentration rapidly decreases to × ~ 0.50 -0.55. In cells 1 and 2 × ~0.55 was observed locally in regions close to the current leads, whereas in cell 4 a radial distribution can be concluded.
A large drop of the lithiation level in the inner and outer regions of cell 4 causes a relatively low mean lithiation state of graphite < x> = 0.68 (5). The cell is therefore characterized by a high standard deviation of the lithiation level x, corresponding to a strong inhomogeneity of lithiation.
Comparing the < x> and x p values obtained for cells 1-4 and taking their relative difference as a measure of lithium homogeneity one may notice that the degree of lithium inhomogeneity is not directly correlated to the cell capacity and decreases in the series 1> 2 ≈ 3> 4, i.e. the cell 1 is characterized by the best value of lithium homogeneity, the overall homogeneity of cells 2 and 3 is nearly similar and the lithium distribution in cell 4 shows the largest deviations from the plateau.
Differences in the lithium concentration profile observed in all four cells of 18650-type can be attributed to the scheme of the electrode connections leading to a non-uniform current distributions. Similar effects were already discussed in literature formerly 25,26 . Depending on the position of the positive and negative tab and the contribution of the current collectors one may expect different resistive paths for the current within a cell. Given that the current collectors possess a similar resistivity per length there is a linear increase of the total resistance of a current path within the current collector with larger distance from the tab position. Three different schemes were found for the investigated cells: Scheme 1 (realized in cells 1 and 2): The positive tab is positioned at the center of the cell corresponding to the end of the cathode current collector, while the negative tab is placed at that end of the anode current collector ending at the outer cell section (Fig. 3). For an unrolled electrode stack one would find the positive tab on one end and the negative tab at the other end. For this configuration the total resistance of any current path from the positive tab to the negative tab is independent of the position l, where the electrode layer is crossed. From this one would expect a homogeneous current distribution/density i versus l. Likewise this configuration should result in a homogeneous lithiation of the graphite during charge and therefore in the charged state of the cell. Scheme 2 (cell 3): In contrast to scheme 1 the positive current lead is located at the middle of the cathode layer (Fig. 3). For the part right of the positive tab in the sketch of the unrolled electrodes a configuration similar to scheme 1 is achieved, where the resistance is independent of the current path from the positive to the negative tab. Therefore one expects a constant current distribution for this part of the cell, from the position of the positive tab to the outside of the cell. However, on the left side of the positive tab the resistance is increasing for positions closer to l = 0 or in other words towards the center of the cell. Due to an increasing resistance with the length of the current path reaches a maximum at the center of the cell (left end of the electrodes). Thus a decrease of the current density i is expected, when moving from the positive tab to the left end of the electrodes. This scheme suggests some lithium inhomogeneity towards the center of the cell. Actually a reduced lithiation is observed at  (Fig. 3). On the one hand this configuration seems to be the most appropriate for high power cells, as it offers a minimum overall resistance between the positive and the negative leads. The maximum length of a current path is only reached for positions at the end of the electrodes. On the other hand it causes the largest deviation for the total resistance of current paths at the outer ends of the electrodes (left and right) or close to the positive and negative tab (middle). During cell operation (either charge or discharge) it may result in an inhomogeneous current distribution, with a maximum at the middle of the electrode layers and minima at both ends of the electrodes. Accordingly a certain lithium inhomogeneity can be expected proportional to the current distribution, and, therefore, higher litihiation levels close to the leads and lower ones at the outside and center of the cell. Although one would expect some kind of equilibration during the rest phase after charging the cell, the measured lithiation levels inside the cell still reflect a behaviour matching the modeled current distribution.
In summary, the homogeneous behavior of cell 1 and cell 2 is then consistent with the contact Scheme 1, the inhomogeneities observed in cell 4 are explained by adopting contact Scheme 3. For cell 3 a behavior resembling a mixture of cell 1 and cell 4 can be realized. This does not hold to be completely true: cell 2 shows remarkable "wings" at the outer cell regions, whereas cell 3 is characterized by a flat (plateau-like) lithium distribution where inhomogeneities are primarily caused by current leads. On the other hand the proposed contact model initially predicts homogeneous lithium distribution in vertical cell direction (along the cylinder axis), whilst authors 32 reported a pronounced (up to 9%) SOC-distribution in LiFePO 4 |C cell nonlinearly increasing from the cell centre (both vertical and horizontal) towards cell edges\surface.
In order to prove this the lithium distribution in cell 4 it was probed at different heights of the sample with respect to its center (− 20 mm, − 10 mm, 0 mm, + 10 mm and + 20 mm) using unchanged beam geometry and collimation. A smaller slit (2 × 10 mm 2 instead of 2 × 20 mm 2 ) for the incoming neutron beam was used in order to increase the spatial resolution in z-direction (along the cylinder axis). The scan at 0 mm is consistent to the dataset shown in Fig. 2c. Results of the z-scans in form of radial lithium distribution at different heights and z-x (z-y) slices are presented in Fig. 4 applying the previously used color scheme. In general the scans performed at different heights obey a similar behavior to that at 0 mm: a narrow (flat) plateau for radial distances from 4 mm to 7 mm with respect to the centre (pin), beyond which the lithium concentration x in Li x C 6 rapidly decreases to x ~ 0.50 -0.55. However, the observed magnitude of the concentration plateaus has been found clearly depending on the cell height (Fig. S7): plateaus at ca. 0.76(1) have been noticed for scans at cell heights of − 10 mm, 0 mm and + 10 mm. Slightly lower graphite lithiation with a maximum concentration of ca. x = 0.74 was observed for the slice at a cell height of − 20 mm, whilst the corresponding scan performed at + 20 mm revealed a maximum lithiation of x = 0.68 only. A significantly smaller lithium concentration in the graphite anode has been found at the top and bottom parts of cell 4, which might be attributed to a mismatch between electrode (cathode and anode) stripes. This indicates an asymmetry of the lithium distribution in the anode occurring in cell 4 at top ( + 20 mm) and bottom (− 20 mm), where lithium concentrations x in Li x C 6 lithiated graphites differ by 0.06. On the first glance it can be attributed to a gravity effect on the electrolyte concentration in the working cell, i.e. its tendency to accumulate at the bottom of the cell. But the genuine reason for this phenomenon requires further research.

Conclusions
The homogeneity of the lithium concentration in the graphite anode (in charged state) of 18650-type cells with different electrode materials and electrode connection schemes was probed using spatially-resolved neutron powder diffraction. The lithium concentration x in lithiated graphites Li x C 6 was determined from the relative ratio between the peak intensities of 002 and 001 reflections of LiC 12 and LiC 6 , respectively.
For a homogeneous cell no spatial variation of the lithium distribution is expected, which would be characterized by a constant (plateau-like) lithiation behaviour. The magnitude of the plateau might directly be related to the cell balancing. In general deviations from the constant plateau-like shape were observed at the regions close to the center pin, the outside of the cells and around current leads, which was associated with the lack of electrode coating for the neighboring cathode sheet. However in some cells deviations from a plateau-like behavior were noticed also in the regions with sufficient electrode overlap (as revealed by X-ray computed tomography).In cell 4 the constant concentration range was only observed in a very narrow range of radial distances (4-7 mm) with respect to the center of the cell. Similar observations hold for cells 2 and 3, but have been found less pronounced.
Observed findings correlate well with the current density distribution expected from tab configurations in cells 1-4. Therefore, most homogeneous lithiation levels are expected for tab configurations realized in cell 1 and/or cell 2, while locating both tabs at the center of the cathode/anode stripes is supposed to result in the highest gradient for the current density inside the cell. Behaviour of cells 1 and 4 can be fairly explained, whereas details observed for cells 2 and 3 as well as the lithium inhomogeneity at different heights of cell 4 can not be fully accounted. Spatial lithium homogeneity might depend not only on the spatial current distribution, but also on electrode materials with complicated morphology, electrolyte, wetting, state and type of solid-electrolyte interphase, various geometrical factors etc. Information about underlying processes defining lithium distribution is crucial for the manufacturing process of safe, robust and high-performance Li-ion cells. It has to be noted that other contributions like current tabs inside the electrode rolls (as shown in cell 3) cause irregularities of the electrode coating and an inhomogeneous winding. In combination with a missing coating of the electrodes at both ends this has an obvious effect on the lithium distribution.
In summary, spatially resolved neutron powder diffraction proves to be a powerful tool for probing the lithium distribution in graphite and has sufficient potential to probe cell homogeneity on a millimeter length scale in various cell designs and states up to a cell thickness of several centimeters. In the current study the stage nature of lithiated graphites in fully charged state has been explored, but the method can potentially be adapted to determination of lithiation degree in commonly-used cathode materials at various conditions. Based on the evaluated lithiation levels for the plateau the cell balancing could be optimized. In combination with efforts to reduce inhomogeneities, e.g. an altered contact scheme of the current collector, the cell capacity could be increased by better usage of the applied electrode materials.

Methods
Electrochemical characterization. Prior to diffraction and tomography experiments "fresh" (formed) commercial cells (Nr. 1-4 types) were "re-activated" by charging and discharging in the nominal voltage windows ( Table 1). Cell cycling was performed using a CCCV protocol at a constant current of 400 mA; In order to achieve equillibrium cell conditons the cutoff current in the CV phase of discharge was set to 1% of 1 C. The battery charge/ discharge was applied using a VMP3 multichannel potentiostat from Bio-Logic and obtained cell capacities were found in fair agreement with the nominal values listed in Table 1.
Computed X-ray tomography. The internal cell geometry with a focus on the current leads was checked by X-ray computed tomography using a v|tome|x s 240 tomography scanner from GE. A direct tube was used at a voltage of 130 kV and a current of 100 μ A. The measurements were carried out in cone beam geometry with a sensitive detector area of 200 mm × 200 mm corresponding to 1000 × 1000 pixels. The magnification was adjusted to get an effective pixel size of 40 μ m. For every sample two data sets of 1000 projections each were collected, where each projection corresponds to an average of 3 to 4 single exposures of 1000 ms. The reconstruction was performed with the phoenix datos|x software. Two 3D datasets (top and bottom part) were merged for every cell using ImageJ. The final visualization of the data including the separation and colour coding of the different parts was performed in VGStudio MAX from VolumeGraphics.
Spatially-resolved neutron powder diffraction. The spatially-resolved neutron powder diffraction study was performed at the engineering diffractometer STRESS-SPEC 33 equipped with translation and rotation stages for sample positioning with respect to the direct and scattered neutron beam. A sketch of the experimental setup is displayed in Fig. S1a. On the first instance a cylindrical cell of the 18650-type can be approximated by two concentric cylinders with an absorbing and scattering medium in between, where the outer cylinder coincides with the steel housing and the inner one is related to the centre pin. Usually the active electrode layers are isolated from the cell housing and the centre pin by several layers of separator. Assuming an isotropic medium, which holds to be fairly correct 26 for gauge volumes larger than 1 mm 3 , the spatial distribution of structural details of the cell organization can be monitored using neutron powder diffraction on the given length scale.
The shape of the gauge volume was defined by the collimation of the incident and scattered neutron beam and the scattering angle. The incident monochromatic neutron beam (25′ horizontal beam divergence, λ = 1.615 Å) was shaped to 2 mm width and 20 mm height using the slit system in front of the sample. The scattered signal was analyzed by a radial oscillating collimator with acceptance width of 2 mm. The 2D position sensitive neutron detector (300 × 300 mm 2 active area, 1056 mm sample-to-detector distance, central scattering angle 26 deg. 2θ ) was used for data collection. Diffraction measurements were performed against time, where the exposure time per pattern was in the range of minutes. A typical 2D diffraction dataset showing sections of Debye-Scherrer rings for two reflections is displayed in Fig. S1b. A homogeneous intensity distribution along the diffraction rings in vertical direction was observed thus indicating that there is no prominent contribution of preferred orientation to the diffraction intensities. The 2D diffraction data were corrected for detector nonlinearities, geometrical aberrations and curvature of diffraction rings. After that they are transformed into conventional 1D diffraction patterns (intensity vs. 2θ angle, Fig. S1c) by data integration.
Analysis of spatially-resolved neutron powder diffraction data. For powder diffraction data collection a central detector angle of 26 deg. (2θ ) was chosen in order to map the structural evolution of 00l reflections from a mixture of stage I (LiC 6 ) and II (LiC 12 ) compounds. For these chemical states of the anode the structure remains unchanged and the relative phase concentration is enough to determine the lithium content in the graphite anode and, correspondingly, the local state of charge. Being extremely sensitive to state-of-charge, state-of-health and temperature the stage I/stage II ratio is often used for monitoring Li-ion cell performance 14,15,29,30,34 . Molar quantities were determined from the diffraction data as m = I hkl [V/F hkl ] 2 [M/ρ], where I hkl is an observed intensity of the Bragg reflection with Miller indices hkl, F hkl is the corresponding structure factor, V, M and ρ are the cell volume, phase molar mass and density, accordingly. The couple of 001 LiC 6 and 002 LiC 12 reflections (defining the interplanar distances between the graphite sheets in corresponding lithium intercalated carbons) is well suited for the quantization of relative phase fractions of stage I and II due to their high-intensity and sufficient separation (no reflection overlap). Using the structure factors F 001 = 4.197 (LiC 6 ) and F 002 = 8.022 (LiC 12 ) along with the cell volumes V LiC6 = 59.676(8) Å 3 and V LiC12 = 112.284(9) Å 3 the LiC 6 to LiC 12 molar ratio (n LiC6 :n LiC12 ) for the case shown in Fig. S1c can be estimated as 0.863(9):0.137 (9). However, for the sake of comparison it is more comfortable to deal with a single parameter and to switch to a relative (mean) lithium concentration x in the lithiated graphite Li x C 6 defined as x = 6/(12-6n LiC6 ) = 1/(2-n LiC6 ). This will yield x = 1 in Li x C 6 for a pure stage I with n LiC6 = 1.0, whilst for a pure stage II x = 0.5. The stability range for stage II is quite broad (x = 0.25-1.0) 14 thus leading to a poor definition when x ≤ 0.5. In the current study traces of 001 LiC 6 reflection were always present, indicating x > 0. 5 in Li x C 6 hereby. For the LiC 6 to LiC 12 ratio shown in Fig. S1c x = 0.880 (7) in Li x C 6 . The standard uncertainty Δ x in Li x C 6 has been found strongly dependent on the Bragg peak intensities so that Δ x typically ranges from 0.002 to 0.032 (mediated by neutron attenuation) with < Δ x> = 0.01. The angular dependence of neutron transmission may only slightly influence the intensity ratio of neighbouring Bragg reflections not exceeding < Δ x> = 0.01 for a linear attenuation coefficient of μ = 0.221 mm −1 .