To corroborate the CMCD observations in the reciprocal space, we employ optical microscopy to directly image an exposed cathode cross-section during battery cycling. The use of optical microscopy is a low-cost, readily available, low-beam damage approach with chemical and dynamic sensitivities to the electrochemical reactions in layered oxide cathodes. The experimental details and analysis of the inter-secondary-particle behaviors have been reported by some of our co-authors26. It has been demonstrated that the optical intensity is monotonically correlated with the SOC16, 17, 18, 27, which is in consistence with our observations. Therefore, we utilize the normalized optical intensity of the NMC particles as a proxy for the local SOC, and the derivative of which is interpretated as the local current density.
A notable feature of the first charge is that the activation of the cathode particles takes place consecutively, which has been extensively discussed in the earlier publication26. In this work, with identification and segmentation of nearly 100 particles, we investigate the microscopic dynamics with statistical significance. It is interesting to point out that, regardless of the size and shape of the particles, it takes around the same amount of time for an individual particle to reach its maximum SOC after its electrochemical onset (Figure S1), while the entire electrode takes ~ 15 hours to complete the first charge in our experiment. This is quite interesting and motivates an in-depth analysis at the sub-particle level.
To understand the intra-particle dynamics, we select a typical particle for a detailed investigation (Fig. 2). As shown in Fig. 2a, when this particle is being actively charged, it gradually increases its optical intensity. By calculating the differential maps using the images with consecutive time stamps, the spatial distribution of the local currents can be visualized. It is interesting that the broadly conceived shrinking-core model is not observed in our experiment. Instead, the particle undergoes a near-surface point onset, followed by a gradual reaction front propagation through the entire particle (Fig. 2b). This holds true for most of the nearly 100 particles imaged in our experiment, featuring a strong statistical representativeness. In Fig. 2c, we further plotted the time-dependent intensity evolution for several randomly selected pixels, which clearly demonstrate the asynchronicity of these local regons. The averaged behavior of the entire particle is plotted in black in Fig. 2c, and its derivative, the local current variation over time, is shown in Fig. 2d, indicating that this particle takes ~ 8 minutes to complete its charging process after its electrochemical onset. This anisotropic delithiation pattern could potentially be attributed to the irregularity in the particle shape28, non-uniform carbon-binder contact20,29,30, grain structure31, compositional variations32 or electrode heterogeneities33, 34, which are often purposely tuned to adjust the cell behavior.
With high-throughput analysis of nearly 100 particles, we single out an interesting one that demonstrates a distinct dynamic behavior. This particle took significantly longer to be fully charged (Fig. 3a). Its optical intensity evolution features a relatively rapid raise near the beginning, which is followed by a slow and steady increase throughout the entire charging process. We enlarge the time window of 20 to 80 minutes for a closer look and, interestingly, observe three steps in this region (Fig. 3b, gray). The corresponding current density plot (Fig. 3b, orange) suggests that there are three charging pulses in this time window, each goes into a different local domain (D1, D2, and D3, respectively). These three charging pulses last for 7, 8, and 11 minutes, respectively. We further extract the normalized intensity profiles for all the four sub-particle domains separately (colored curves in Fig. 3c). Based on this analysis, we identify several critical time stamps for this particle, namely at 38 minutes, 45 minutes, 54 minutes, 65 minutes, and the end of charge, 891 minutes. The particle’s electrochemical activity maps (Fig. 3d to g, see Methods section in supplementary materials for more details) highlight four different local regions that are being actively charged by different charging pulses. The respective behaviors of domains D1, D2, and D3 are very similar to that of a normal individual particle. Domain 4 (D4), on the other hand, shows a gradual charging behavior over the course of ~ 15 hours, which is quite rare during the first charge. The varying spatial distribution of the local current over this particle is further highlighted by the circular plot in Fig. 3d to g, which demonstrate the sequential activities in the four domains.
To gain a mechanistic understanding of the observed intra-particle asynchronous domain activity, we conducted finite element analysis using a strongly coupled multiphysics theoretical framework26. Detailed information regarding the FEM model and material properties can be found in the Methods section. Figures 3h to j depict that the interaction between electrochemistry and mechanical damage controls the intra-particle charging behavior (see Movie S2 for more details). The seven domains inside this NMC particle represent intra-particle regions that exhibit chemical homogeneity at the initial state, with directionally isotropic transport and mechanical properties (elastic stiffness and strain). At the particle surface, where the porous carbon binder partially covers the particle, interfacial reactions based on the Butler-Volmer equation occur. The interfacial reaction occurs only at the carbon binder and NMC interface, and lithium diffusion across the domain is controlled through Fick’s diffusion law. During charging, as lithium ions deintercalte the host lattice, NMC domains experience net lattice volumetric reduction. Consequently, the shrinking of domains generates stress at the domain boundaries, and mechanical detachment/damage occurs. As shown in Fig. 3h to j, the developed heterogeneous nature of the domains within the particle leads to the creation of non-uniform charging behavior among different regions. This asynchronous and heterogeneous charging behavior triggers lattice mismatch, which results in the progressive development of damage. These damages continually restrict and redirect the flow of lithium within the particle, thereby affecting the charging behavior. In Fig. 3k, we observe a monotonic decrease in net particle lithium concentration, with considerable differences in lithium concentration profiles among all seven domains at different time points. These asynchronous reactions are accompanied by an increase in mechanical damage (Figure S2). The simulation results clearly demonstrate the dynamic nature of the lithium diffusion pathways, which is modulated by intra-particle damage and heterogeneous interfacial reaction with carbon binder. In other words, the relative difference in electrochemical and mechanical properties among domains and the degree of incomplete carbon binder coverage of NMC particles can influence the intensity of such heterogenous intra-particle behavior in composite electrodes. This heterogeneity also leads to further damage of the NMC particles.
We declare that the presented simplistic 2D model does not entirely replicate the experimental configuration. In the real world, the unique sub-particle level dynamic features would be more complicated, involving cracks, defects, impurities, and compositional heterogeneity35. Understanding the origin of these heterogeneities is crucial for improving the performance and reliability of batteries. The identification of distinct sub-particle level dynamic features suggests that battery performance and reliability may be improved by tailoring the internal structure and composition of cathode particles to reduce heterogeneity in electrochemical reactions36, 37.
The co-existence of several local domains with very different electrochemical behaviors within an individual secondary particle echoes our CMCD observations of the asynchronicity in microscopic domain dynamics. Their co-evolution could have a role to play beyond the first charge. To investigate this further, we followed another particle throughout its second cycle (Fig. 4). Direct visual assessment of the electrochemical activity maps at different SOC clearly reveals not only the spatial heterogeneity but also the temporal asynchronicity. Two domains with distinct behavior are identified (D5 and D6). D5 initiates its charging process at the beginning of the cell charge and reaches to its maximum SOC well before the end of the cell charge. Interestingly, while the whole cell is still being charged between time stamps 3 and 4, D5 already reverse its trend and starts to demonstrate a local discharging behavior. D6, on the other hand, shows a considerable delay. Its charging behavior starts around time stamp 3 and persists until time stamp 6, which is already significant into the cell discharging process. The evolution of the probability distribution of the local electrochemical activity is plotted in Fig. 4b, showing that the orange peak (for D6) is lagging and following the purple peak (for D5). These asynchronous activities of D5 and D6 are schematically illustrated in Fig. 4c. Collectively, these domain dynamics govern the cell behavior, which is essentially an integration of all the active domains in the battery electrode. Based on this observation, it can be inferred that diffusion resistance among different local domains could affect the asynchronous domain dynamics and equilibration. A well-designed cathode material that enables a smooth interdomain lithium diffusion and charge transfer and, thus, facilitates a rapid domain equilibration could effectively suppress the stress accumulation and prevent particle cracking.