A schematic of the experimental set-up is displayed in the upper part of Fig. 1a, where we can see the commercial cell (NMC811 vs. graphite) paired with an US emitter receptor. The cell was placed on a custom sample holder [20] developed by Lyonnard’s group (Fig. 1b) allowing lateral scans with the microbeam. The cycling protocol of the cell is displayed in Fig. 1c where three cycles were performed, a first one at C/3 rate, to verify the cells behaviour, a second one at C/8 rate and the last one cycled by galvanostatic intermittent titration “like” technique (GITT) using C/3 rate charging with 1 h pulses and 1 h relaxation steps.
Prior starting the operando measurement, we first verified that the measurement of the ultrasound signal was not impacted by temperature fluctuations of the hutch (see Note 1, Supporting Information) and that the measurement at the beamline is similar to the one obtained in the laboratory (Note 2, Supporting Information). We also checked that the cell did not suffer from X-ray damage (see Note 3, Supporting Information). Then, scalar quantities are extracted from the ultrasound signal as mentioned in Note 4 (Supporting Information), especially the time of flight (ToF), one of the most relevant parameters when it comes to ultrasound investigation in batteries. However, the interpretation of ToF’s evolution during cycling is far from trivial as it can be seen in Eq. 1. The ToF depends on constants such as the total mass and the section probed, and on two variables being the thickness and the average effective elastic modulus (Young’s modulus) of the entire medium crossed by ultrasound waves.
$$\:\begin{array}{c}ToF=\frac{d}{\sqrt{\frac{{E}_{eff}}{\rho\:}}}=\sqrt{\frac{m}{A}}\sqrt{\frac{d}{{E}_{eff}}}\:\#\left(\varvec{E}\varvec{q}\varvec{u}\varvec{a}\varvec{t}\varvec{i}\varvec{o}\varvec{n}\:1\right)\end{array}$$
with m (kg) being the mass probed, A (m2) the section, d (m) the thickness of the sample probed, Eeff (kg.m− 1.s− 2) the sample’s effective Young’s modulus and ρ (kg.m− 3) the sample’s average density. Both variables depend on the sample’s chemical composition and environmental conditions, thus temperature and pressure.
To understand the ToF evolution and determine which parameter(s) dominate it, a measurement of the cell’s thickness and the global effective elastic modulus (referred to as Young’s modulus in the following) evolution is needed. For the former, several papers denote a thickness increase of several tens of micrometers during the charge of NMC vs. graphite cell (graphite lithiation and NMC delithiation) depending on the number of electrode layers [21], [22], [23]. During charge, NMC particles undergo an anisotropic variation in lattice parameters resulting in a global contraction of the lattice volume [24]. At the NMC electrode level, it is difficult to determine if the volume change of NMC particles has an impact in any of the three dimensions of the electrode, since the electrodes are porous. For graphite particles, the variations in lattice parameters are also anisotropic and induce a large volume change along the c-axis of the crystal structure (corresponding to graphene layers d-spacing), this volume increase reaching 14%. At the graphite electrode level, the increase of the interlayer d-spacing upon lithiation results in a unidirectional expansion, easily visible even at electrode level [25], [26]. Based on this discussion, the global expansion of a full cell NMC | graphite upon charge is attributed to the more dominant volume expansion, i.e. the graphite electrode. To verify our hypothesis, we dismounted, in an Ar-filled glovebox, the same commercial cell as the one brought to the synchrotron, washed the electrodes, gently removed a side of the double coated electrodes and mounted these single-side coated electrodes in a nano-dilatometer allowing the operando monitoring of the working electrode (graphite) expansion in a three electrodes set-up (NMC811 as counter electrode and Li metal as reference) (see Note 5, Supporting Information). In Fig. 2, we compare the full cell’s thickness evolution (evolution at 25 kPa constant pressure) reported from literature [23] with the graphite electrode expansion measured by operando dilatometric measurement as mentioned above. Figure 2a shows the full cell voltage and the negative electrode potential of the graphite vs. NMC cell. Qualitatively, we notice a very similar evolution of the full cell expansion (Fig. 2b) with the graphite electrode expansion (Fig. 2b), supporting the fact that graphite dominates the full cell thickness evolution in the low SoC regime (below 80%). Note that above this value, the positive electrode starts contributing with a negative expansion that is correlated to the shrinking of the NMC c-axis parameter at high voltage. At SoC > 80% it seems that graphite and NMC effects are counter balancing each other, yielding almost unchanged full cell thickness.
Based on Eq. 1, we need to consider the Young’s modulus evolution of the electrode materials (positive and negative). This metric is more complicated to assess, even more upon lithiation and delithiation of active materials, since generally, this kind of measurement should be performed without any current collector, which is far away from a realistic electrode and is hardly performed operando due to liquid electrolyte state. Generally, the Young’s modulus describes the elastic properties of cell components, so in the present case we have three Young’s modulus in total, two coming from each electrode and one coming from the separator, the last one could be attributed to cell housing (generally negligible). From the literature, graphite, NMC811, and most electroactive materials display a Young’s modulus that increases upon lithiation [27]. According to Davies et al., graphite displays an important variation of electrode sound velocity compared to LiCoO2 (LCO) [28], explaining the larger impact of the graphite electrode on the average Young’s modulus of a full cell. A reason for this statement resides in the “softness” of graphite that can deform reversibly upon de/lithiation [29] compared to lithium oxide materials such as LCO [30] or NMC [27] that easily break upon mechanical strains during cycling.
Based on Eq. 1, an increase in cell thickness upon charge should lead to an increase in the ToF signal which is not consistent with our observation (Fig. 3) nor with the literature [22] showing a decrease in ToF upon charge.
The remaining parameter that evolves during cycling and could dominate ToF signal evolution is the effective Young’s modulus of the cell. During cycling, the Young’s modulus of graphite electrode becomes stiffer upon lithiation, while NMC’s Young’s modulus is decreasing upon delithiation which matches our results indicating that the ToF signal evolution of the full cell is dominated by the Young’s modulus evolution of graphite electrode. However, the non-monotonic evolution of ToF measurement highlights a complex correlation between the various cell components and additional correlation can be found with both electrode materials.
To deconvolute the ToF signal and attribute each process to a specific electrode phenomenon, we carried out a full structural investigation through operando XRD measurements while simultaneously collecting the ToF signal. As can be seen in Fig. 4, we can follow the evolution of the main NMC811 (003) Bragg reflection, as well as the one of graphite (002), both being properly separated in the q range. Using a home-made software for data correlation based on MatPlotLib Python library [31] and LMFIT package [32], we fitted both main reflections (raw data plotted in Note 6, Supporting Information) and obtained the average peak position to follow the structural changes through the evolution of the c-axis. The c-axis is indeed interesting for graphite electrodes since the de/lithiation processes modify the interlayer d-spacing [33]. For the NMC811 electrode, a-axis and c-axis are evolving. However, as shown by Quilty et al. [24], the c-axis is the one evolving the most, especially during the H2-H3 phase transition, justifying the choice of this axis too.
Although we can separately access the structures of the positive and negative active materials in given states of lithiation by XRD and access the c-parameter variations, the deconvolution between the positive electrode and the negative electrode reaction mechanisms is far from trivial since, in a full cell, the electrochemical signatures of both electrodes are merged. To address this issue, the operando nano-dilatometric measurement allowed a deconvolution of each electrode potential to understand the global electrochemical signature as well as correlate the volume change of one electrode at a time as a function of the cell’s voltage.
From the three-electrode set-up, the full cell voltage evolution is a linear combination of the two electrodes' potential variations vs. lithium metal as described in Eq. 2. However, the weight of each electrode potential variation in the full cell voltage evolution varies depending on the de/lithiation processes. Then, from the classical galvanostatic curves E = f(t) obtained, we can define weighting functions for both electrodes (Eq. 3), ωgr and ωNMC, which quantify the impact of the electrode potential variation on the full cell voltage. Those weighting functions are plotted in Fig. 5a.
$$\:\begin{array}{c}\frac{dE}{dt}\left(FC\right)=\:\frac{dE}{dt}\left(Gr\right)+\frac{dE}{dt}\left(NMC\right)\#\left(\varvec{E}\varvec{q}\varvec{u}\varvec{a}\varvec{t}\varvec{i}\varvec{o}\varvec{n}\:2\right)\end{array}$$
$$\:\begin{array}{c}{\omega\:}_{Gr}=\frac{dE}{dt}\left(Gr\right)/\frac{dE}{dt}\left(FC\right)\:and\:{\omega\:}_{NMC}=\frac{dE}{dt}\left(NMC\right)/\frac{dE}{dt}\left(FC\right)\#\left(\varvec{E}\varvec{q}\varvec{u}\varvec{a}\varvec{t}\varvec{i}\varvec{o}\varvec{n}\:3\right)\end{array}$$
Additionally, the dt/dE (analogue to Differential Voltage Analysis under galvanostatic conditions) evolution (Fig. 5b) provides information about the nature of the de/lithiation mechanism. For NMC811 material, lithium de/insertion takes place through a solid solution mechanism with some lattice parameter changes inducing specific variations in the dt/dE representation. For graphite, biphasic mechanisms generate peaks in the dt/dE curve, while solid solution ones do not impact significantly the curve evolution. Here, we should note that the three curves are linked through the following Eq. 4:
$$\:\begin{array}{c}\frac{dt}{dE}\left(FC\right)=\frac{1}{\frac{dE}{dt}\left(Gr\right)+\frac{dE}{dt}\left(NMC\right)}\#\left(\varvec{E}\varvec{q}\varvec{u}\varvec{a}\varvec{t}\varvec{i}\varvec{o}\varvec{n}\:4\right)\end{array}$$
With this representation, we can assign each phenomenon to a particular electrode [34]. The region A of the charging cell (Fig. 5) displays a voltage variation imposed by the graphite electrode as showed by the weighting function curve. However, in this region, no specific peak in the dt/dE representation is observed since solid solution like processes of lithiation/delithiation are taking place for both electrodes. In region B (Fig. 5), the full cell voltage variation is due to both electrodes’ potential variation but the similar shape of peak observed for the full cell and graphite curves in the bottom graph indicates that graphite undergoes its first biphasic phase transition (stage IV, formation of LiC36 from LiC72 [35]) while NMC811 is still delithiated through solid solution mechanism with no phase transition. A similar observation can be made in region C (Fig. 5), where the apparent peak of the full cell curve is attributed to the graphite stage III liquid-like transition. In region D (Fig. 5), two processes can be found, in the first one (depicted as D1) graphite undergoes the stage II phase transition (LiC12 phase formation) while the hexagonal to monoclinic phase transition H1 to M of NMC811 takes place (common phase transition for the NMC811 compounds family). The second subregion (D2) contains still the graphite stage II formation, but NMC811 stops its specific phase transition. The last region E (Fig. 5) displays first a local minimum attributed to the change in biphasic mechanism of graphite from LiC12 formation to LiC6 formation (stage I) and the NMC electrode undergoes first the monoclinic to second hexagonal phase transition M to H2 in the E1 subregion. The second characteristic peak in subregion E2 correlates with the second to third hexagonal transition H2 to H3, involving Ni3+ and Co3+ oxidation into Ni4+ and Co4+ [36].
Now that we better understand the electrochemical response of a full cell (graphite vs. NMC811), we correlate the ToF signal to the electrochemical processes of NMC811 and graphite and to their structural evolutions obtained from operando XRD. Figure 6a and Fig. 6c plot the evolution of peak position from (002) Bragg reflection of graphite and (003) Bragg reflection of NMC811, respectively, while Fig. 6b plots the evolution of dt/dE from the full cell voltage variation obtained at C/8 correlated to the ToF evolution. Similarly, as for the dilatometric data, we divided the obtained correlated ToF-XRD data by the five same regions of interest, beginning by the cell discharge (graphite delithiation and NMC811 lithiation). It should be noted that the boundaries of potential are a bit shifted between Fig. 5 and Fig. 6 due to the high separator resistance used in the nano-dilatometric measurement. Starting with the region E2’, we noticed a huge variation of ToF signal correlating with the peak attributed to the H3 ◊ H2 transition of NMC811, inducing an important change in lattice volume of NMC811 estimated at 4% by Quilty et al. [24], which is the most important variation in volume among all the transitions of NMC811. For the rest of the region E’, from 4.1 V – 3.8 V, a linear increase of ToF signal is observed and could be attributed to the LiC6 phase disappearance. No obvious change in slope of ToF signal are noticed when NMC811 is undergoing the H2 ◊ M phase transition (E1’ subregion). For the D2’ region, we noticed a change in slope of ToF signal compared to the E1’ region which correlates with the change in biphasic mechanism of graphite (from LiC6 to LiC12). When NMC811 undergoes the M ◊ H1 phase transition in the subregion D1’, an increase in the slope intensity of ToF is observed. For the third region (C’), a period of pseudo-plateau in ToF signal correlates with liquid-like stage III of graphite. In this range of potential, very small amount of charges is exchanged leading to very small variation in global Young’s modulus of the cell. Concerning the B’ region, again, a specific slope in ToF increase is observed correlating here with the stage IV of graphite disappearance. For the last region of interest (A’), a kind of plateau in ToF signal correlates with very low amount of lithium stored through solid solution like mechanism in both materials.
Regarding the previous discussion, we observed that ToF signal is mainly driven by the Young’s modulus evolution of graphite, leading to a decrease in ToF during charge of the cell and an increase during the discharge. By deconvoluting all the electrochemical processes, we also demonstrated that most of the modulations in ToF are correlated to structural changes of the two electrode materials. We notice that, in general, the ToF changes linearly depend on the potential in the selected regions, with varying positive or negative slopes that might reflect a correlation between specific phase transitions and ultrasounds properties. This effect could be further investigated on simplified systems to obtain a quantified relationship between ToF and electrochemical processes, which is beyond our scope here.
Additionally, galvanostatic intermittent titrations “like” technique (GITT) was implemented to understand the impact of the relaxation in the structure of the material as well as in the ToF signal. This electrochemical characterization is widely used in the field of batteries [37], consisting in the application of regular pulses of low-intensity current. A long resting period is then left between those pulses, allowing the cell to relax and the electrode material to reach a state of quasi-equilibrium. Lithium concentrations are then equilibrated in the cell leading to homogeneous lithium distribution throughout the electrode volume. Here a current pulse equivalent to C/3 rate was selected to generate heterogeneities in the cell, to understand the sensitivity of the ultrasound technique.
As depicted in Fig. 7, ToF signal is sensitive to the relaxation processes occurring in the cell, since during the resting periods of the cell (R1 – R6 on Fig. 7) the ToF signal is not constant. Several parameters could contribute to the relaxation processes, here we tried to discriminate the one dominating this relaxation: i) the temperature, imposed by a cycling at C/3 rate, could generate local heating that the cell needs to dissipate and could impact the ToF signal. From previous experiment performed in the laboratory with similar commercial batteries, an extrapolated dependency of 26 ns/K was estimated. Considering most of the relaxation phases (expect R3 which displays higher relaxation), an average decrease of 12 ns is observed in the ToF signal, corresponding to a temperature decrease of ~ 0.5 K, a realistic temperature change in such system at C/3 rate [38]. ii) Young’s modulus of graphite, which was the dominant factor for the overall evolution of ToF signal upon cycling, could also be sensitive to relaxation processes due to the local inhomogeneous reaction imposed by the C/3 rate. For this purpose, specific attention was given to the diffraction patterns of NMC811 and graphite evolution during relaxation periods (see Note 7, Supporting Information), specifically the last relaxation period related to high state of charge of the cell meaning the transition LiC12 to LiC6. Figure 8 represents the X-ray diffraction patterns of the two electrodes during the R6 relaxation involving H3 phase of NMC811 and a mixture of LiC12 (~ 1.78 Å−1) and LiC6 (~ 1.70 Å−1) graphite phases. Here, we did not notice any visible shift of the peak position for NMC811 leading us to conclude that NMC is not presently affected by a relaxation process at C/3 rate. Concerning the graphite, no changes in LiC6 peak intensity, nor in its shape (constant FWHM) are observed and the small variations in LiC12 contribution are noise-related, indicating a constant structural state during the relaxation process. However, we may note that, since the numerous layers probed in the battery induce large diffraction peaks, small variations in local state of charge of both electrode material could be hard to detect in the present through-plane-averaged XRD setup. Therefore, we cannot strictly exclude an eventual relaxation in Young’s modulus that should impact the ToF. Nevertheless, our observations prove that Young’s modulus is not the main factor driving the ToF relaxation. iii) Mechanical relaxation of the electrode, coming from the electrode volume changes, could also impact the ToF signal. So far, the binder signal was not detectable due to its very small contribution (less than 2wt.% in each electrode). Then, in the present case, temperature relaxation due to ohmic relaxation takes part in the ToF relaxation, and we cannot exclude a contribution from mechanical relaxation, this one coming from the electrode binder or the separator.