1. Introduction
Lithium metal holds promise as an anode material for next-generation batteries due to its high specific capacity and low electrochemical potential, which could lead to energy storage devices with unprecedented energy densities. However, a significant challenge during charging of these batteries is the formation of lithium dendrites, which can pierce the separator, posing the risk of short-circuiting and thermal runaway [
1].
Lithium dendrite growth in lithium metal batteries, can be primarily attributed to the uneven deposition of lithium during the charging process. This process occurs at the anode, where lithium ions are reduced and deposited as metallic lithium. The concentration polarization phenomenon is crucial in the context of dendrite growth. During charging, lithium ions are consumed faster at the anode surface than they are replenished from the bulk electrolyte, leading to a localized drop in ion concentration. This drop increases the overpotential required for lithium deposition, promoting uneven lithium growth and dendrite formation [
2]. Strategies like using electrolytes with higher lithium salt concentrations, adding specific electrolyte additives, and designing advanced anode structures are employed to mitigate the effects of concentration polarization and thus reduce dendrite growth [
2,
3,
4,
5,
6,
7,
8,
9].
In situ analytical techniques, such as NMR, MRI or X-ray tomography, have become invaluable for the advancement of battery research. These techniques provide crucial insights into the dynamics of dendritic growth, offering the potential to guide the development of rechargeable lithium metal batteries [
10,
11,
12]. Chemical shift NMR is particularly valuable for detecting dendritic structures, allowing to distinguish between different electrochemical environments within the cell [
9,
13,
14,
15]. Moreover, the use of
6Li and
7Li isotopes in MRI facilitates direct imaging of dendritic growth. The
7Li MRI, in particular, provides real-time visualization of dendritic formation, crucial for studying the effects of various operating conditions on dendrite development [
12,
16,
17,
18]. Recently, indirect MRI was used in detecting dendrite growth via their effects on the surrounding electrolyte [
19]. The effects studied were the displaced volume of the electrolyte, the magnetic susceptibility effects on the electrolyte, and the RF field effects on the electrolyte signal. As the electrolyte is typically composed of proton-rich material, the
1H signal is much stronger, compared to
6Li and
7Li MRI. This method allows for fast 3D MRI experiments with high resolution, providing a more comprehensive understanding of dendrite growth dynamics. In addition to imaging the evolution of the morphology of deposited lithium in a cell, MRI has also been used to characterize electrical currents within cells [
20]. NMR and MRI techniques provide a multifaceted approach to understanding and mitigating dendritic growth in lithium batteries. The integration of direct and indirect imaging methods, alongside the capability to characterize internal currents, positions these techniques as indispensable tools in the quest for safer battery technologies.
This study focuses on the comparison of dendrite formation in different circuits of battery cells – single cell, cells in parallel and cells in series. Suppose there is a cell with more dendrites than other cells. This cell would then have a lower electrical resistance than the other cells, which would cause an even faster growth of dendrites in this cell in the case of cells connected in parallel due to a constant and equal voltage on all cells, while it would not cause a change in the growth of dendrites compared to the other cells in the case cells connected in series due to the same current in all cells. Our hypothesis is therefore that the growth of dendrites in series-connected cells is self-regulated, i.e., the same in all cells and there is no cell with excessive dendritic growth, while this is not the case for parallel-connected cells. This study was performed on symmetrical Li-metal electrode cells in three different electrical circuits that were sequentially imaged by 3D MRI during the charging process. The measured 3D images of the electrolyte in the cells were then analyzed to extract the morphological development of dendrites and identify the effects of the circuit on the growth patterns. This approach allowed to evaluate the dendrite growth mitigation performance for each type of circuit, which could facilitate the design of safer and more efficient lithium metal batteries.
2. Materials and Methods
2.1. Sample Preparation
The battery was assembled in one of three different configurations: as a single cell, cells in parallel, or cells in series, as shown in
Figure 1. The battery body components were fabricated using stereolithography (SLA) 3D printing technology using a Prusa SL1S printer with Prusament Tough resin (Prusa Research, Prague, Czech Republic). For the electrodes, two square 0.3 mm thick lithium metal plates were stamped from a roll of lithium foil. The electrolyte was a solution of 1 M LiPF
6 in a volumetric mixture of ethylene carbonate (EC) and dimethyl carbonate (DMC) in a 1:1 ratio. Strips of copper foil were used as current collectors. These strips were stretched through specific slots in the cell housing and then sealed. The lithium electrodes were inserted into the corresponding slots of the housing during the assembly of the battery. This was then filled with electrolyte and sealed with a screw top. Battery assembly was performed in an argon-filled glove box (Vigor Gas Purification Technologies, Marktheidenfeld, Germany). Before inserting the battery into the MRI magnet, it was connected to a constant current source.
2.2. Magnetic Resonance Imaging
The MRI experiments were performed using a 9.4 T (400 MHz proton frequency) wide-bore vertical magnet (Jastec, Tokyo, Japan), equipped with a Bruker Micro 2.5 gradient system (Bruker, Ettlingen, Germany) and controlled by a Tecmag Redstone spectrometer (Tecmag, Houston, TX, USA). Data acquisition was performed with a 30 mm quadrature Bruker RF probe operating in linear polarization mode.
For optimal imaging, the cell was aligned inside the magnet so that the B0 and B1 magnetic fields were parallel to the electrode faces. 3D 1H MRI was performed using a 3D RARE sequence with parameters: inter-TE of 5.8 ms, TR of 2030 ms and ETL of 8. The field of view was set to 24×24×12 mm across the x, y, and z axes, with a 128×128×64 signal acquisition matrix in directions: read (x), first phase (y) and second phase (z). This setup yielded a 3D image that encompassed both electrodes and the electrolyte region, with an isotropic resolution of 187.5 μm.
Dendrite growth began concurrently with the start of 1H MRI image acquisition, maintaining a constant charging current of 1 mA throughout the experiment, with an upper voltage limit of 4.0 V. Each image acquisition lasted approximately 2 hours and 20 minutes, followed by a 1-hour pause before starting the next image acquisition. In total, for each battery configuration (single cell, cells in parallel, cells in series), 20 sequential images were acquired, which corresponded to 67 hours of experiment time.
2.3. Image Analysis
The first step in image analysis was to determine an intensity threshold indicating the boundary between dendritic and pure electrolyte voxels. This was done by image histogram analysis, which involved calculating the histograms of all images in the sequence, followed by subtracting the histogram of the base (first) image from these. These 'difference histograms' were analyzed to identify a transition point - where the number of low-intensity voxels increases and high-intensity voxels decreases - which represents the threshold intensity. After determining the threshold intensity for each image, inverted thresholding was performed, where voxels with values below the threshold were assigned value of 1 and other voxels were assigned a value of 0. These masks were then manually corrected for possible inclusion of artifacts due to gas bubbles or other image imperfections.
Further analysis was focused on extracting the apparent and more accurate dendrite volume as a function of battery charge time. The apparent volume of dendritic structures
was simply equal to the sum of all dendritic voxel volumes, i.e., voxels with intensity of 1 in inverted-threshold MR images. As this results in an overestimation of the amount of lithium transferred between the electrodes during charging, this estimate must be further calibrated. In this study, this was done using a simple model in which all dendritic voxels are assumed to contain the same volumetric fraction
f of lithium at a given charging time. This fraction was then determined by equating the product
, where the sum includes all cells of all battery configurations, with the corresponding sum of the dendritic volume transferred between the electrodes according to the laws of electrolysis in all cells of all battery configurations
. From here follows
Note that the
f-factor thus determined is calibrated for each time point and can in principle change with charging time. Knowing the
f-factor makes it possible to estimate the amount of lithium transferred between electrodes for any cell of any battery configuration
3. Results
Figure 2 (Columns A, C, E) displays sequential 3D MR images of dendritic growth in symmetrical lithium-metal electrode cells over time. The growth was due to a constant total current of 1 mA through the circuit and was the same in the three cases studied. These images show a representative slice (
xy plane) through the cells, which best displays the morphology of the dendrite structure development across the cells. In the MR images, the electrolyte is bright, while the dendrites correspond to signal voids. In the case of a parallel circuit of cells (Column C), dendritic growth patterns are observed within the compartments formed by the dividing walls. The third cell shows notably less dendrite development compared to the other cells. In the case of a series circuit of cells (Column E), dendritic formations appear more uniform across the cells, with cells 1 and 2 showing initial rapid growth. Gas pockets are visible as dark circles on MR images and can thus be easily distinguished from dendrites, which have more irregular shapes. The presence of gas pockets and their size vary according to different types of circuits. In a parallel circuit, these pockets remain small and are confined within individual compartments. The signal intensity variation in proximity of the metal electrodes, which is best seen in the single cell and in parallel cells, is mostly due to the eddy currents induced by the
B1 filed and the susceptibility effects of the
B0 field. These effects were already minimized by proper orientation of battery cells in the MRI magnet, as described in [
21,
22], but their complete elimination is not possible.
In order to optimally determine the image threshold required for the segmentation of dendritic structures from MR images, the temporal evolution of the image intensity distribution was also analyzed. For the single-cell experiment, this analysis is shown in
Figure 3A, with an image of stacked histograms in which each color-coded row corresponds to an image histogram at a given time point, i.e., the row corresponds to the time point and column to the bin index (image intensity), while brightness (color) to the voxel count of the bin. From this image can be seen that the image histograms at the beginning (
t < 20 h) have one dominant line around voxel intensity 100, corresponding to electrolyte-rich regions, and one emerging line around intensity 10, corresponding to the deposited lithium on the anode. This line becomes much more intense with time, which can be explained by the increasing amount of deposited lithium. In addition to this dominant line, two significantly weaker lines appear around intensity 70 and 120. These two lines later merge into one broader line that reaches an intensity of 220 and narrows to an intensity of 180 after
t = 40 h. The image threshold was then determined from the difference between the current and baseline (first) image histograms as the transition point (intensity) between intervals of increasing and decreasing voxel counts. For the histograms in
Figure 3A, this analysis is shown in
Figure 3B with a color-coded image of stacked histogram differences in rows. This image has two distinct regions, one with an increasing voxel count (in red) and the other with a decreasing voxel count (in blue), that are separated by a transition region with a stationary voxel count at intensity of about 50, corresponding to the threshold. Note, that the transition between these two regions is less clear with first few images, as there is little change in the development of the image then.
Once the optimal thresholds were determined for each MR image of the experiment, dendritic structures were segmented using inverted thresholding (
Figure 2, columns B, D, F), in which dendrites and electrodes were assigned a binary value of 1, and the rest of the image was assigned a binary value 0. This allowed quantification of the apparent dendrite volume, which was calculated as the sum of the volumes of all dendritic voxels (binary value 1) of all slices in the region between the electrodes. The graphs in
Figure 4 show the resulting apparent dendrite volumes (
) per electrode surface area (
) as a function of time for cells connected in parallel (panel A) and for cells connected in series (panel B). For reference, each graph also contains this dependence for a single-cell experiment (blue circles). Plot of the growth for a single-cell shows a deviation from the linear growth, representing a sigmoidal trend, with a steep rate of growth in the first 19 hours, then slowing down and increasing again around the 40
th hour of measurement. Plots of dendrite growth for cells in parallel show different rates of growth, with cell number 3 showing minimal increase. A plot of the average growth of cells in parallel (red triangles) shows the growth with the best linearity. In this case, the growth is similar to that in the single-cell experiment. Plots of dendrite growth for cells in series were more similar to each other and also to growth in a single-cell experiment. The exception is cell number 4, which showed a rapid initial growth rate that latter approached the growth rates of the other cells. As expected, the growth uniformity was much greater for cells in series than for cells in parallel. Note that because in the series cell configuration, dendrites started to reach the opposite electrode after about 28 hours. For this reason, experimental points for cells in series do not exist after this time point. The dashed line in the graphs corresponding to the linear growth model is derived from Faraday’s law of electrolysis for constant charging current of 1 mA.
Figure 5 shows the development of the anode plating in a representative slice across a single cell with images taken at 20-hour intervals. The images are the result of processing that involves subtraction of the current image from the baseline (first) image, followed by normalization of this result to the to the average intensity of the electrolyte and final multiplication with a mask created by the inverted thresholding (
Figure 2). This results in an image with intensity 0 in the electrolyte region and intensities from the threshold to 1 in the dendritic voxels. Here, an intensity of 1 (white voxel) indicates that a significant part of the voxel is filled with dendrites, so that in this voxel, no MR signal was produced due to the physical presence of the dendrite or its influence on the magnetic fields around it and thus disturbed signal reception. A dendritic voxel with intensity less than 1 (gray voxel) is in a smaller part filled with dendrites. Apparently, there is a positive correlation between the intensity of a dendritic voxel and the proportion of dendrites in it, but a more precise relation between the two is difficult to find without a theoretical analysis of the MR signal around electrically conductive structures of the dendritic shape [
19].
A comparison between the
plots (different symbols) and the theoretically predicted volume of deposited lithium per electrode surface area (dashed line) in the graphs of
Figure 4 shows that the latter is significantly smaller for any cell of any battery configuration. This discrepancy is due to dendritic voxels being only partially filled with lithium. We propose a simple solution to this discrepancy, namely a correction based on the assumption that all dendritic voxels are filled with lithium with the same fraction
f. This
f-fraction was determined for each time point as the ratio between the deposited lithium per electrode area according to the theory divided by the average of
over all cells in all experiments (Equation 1). With known
f the apparent volume of dendritic structures
can be replaced by a more precise measure for the volume of dendrites
equal to
. With this correction/calibration the plots of
from
Figure 4 transform to their calibrated equivalents of
as a function of time for all cells and battery configurations that are shown in
Figure 6. From these plots can be seen that in average these values are close to the theoretical values for the volume of deposited lithium per electrode surface area. In addition, it can be clearly seen that plots of
are much more uniform for different cells in case of cells in series (panel B) than in case of cells in parallel (panel A).
Graph in
Figure 7 depicts dependencies of the
f-factor and of the signal
S against charging time
t. The signal
S corresponds to an average MR signal of dendritic voxels normalized to the signal of the electrolyte. Form the graph it can be seen the voxels with lower filling factor
f produce higher signal
S, such voxels are found mostly at the beginning of the experiment, while the voxels with higher filling factor f produce lower signal
S, which happens mostly in the second half of the experiment. It can also be seen that both, the
f-factor plot as well the signal
S plot flatten after
t = 40 h.
4. Discussion
The dark, branched shapes in MR images in
Figure 2A,C,E represent the growing dendritic structures on the anode of symmetric lithium cells. As can be best seen from MR images of a single cell (
Figure 2A), initially the deposits are dense, moss-like, but in the later stages of the experiment they develop a more branched structure. This is also seen to some extent in the parallel cell experiment (
Figure 2C), but interestingly not so much in the series cell configuration (
Figure 2E). This may be the result of the development of the lithium salt concentration profile in the electrolyte, which in turn dictates the start time of dendritic growth, as per the Sand's model [
8]. As the distance between the electrodes shortens with time, the concentration gradient established with time is steeper, which may explain the grater branching of the later structures. Dendrites tend to grow towards regions of higher ion concentration, as these regions provide more material for deposition [
8,
23]. A steeper concentration gradient can thus cause more pronounced and faster growth of dendrites. Another reason for the faster growth may be the increase of the electric field in regions with smaller interelectrode distances. For example, Chazalviel’s model states that the growth rate of dendrites is proportional to the factor of anion mobility and the electric field (
) [
8].
The distinct transition between the dense and branched morphology can be explained by the fact that high current densities lead to inhomogeneous lithium-ion flux at the electrode-electrolyte interface, leading to uneven lithium deposition and the formation of dendritic structures. At lower current densities, lithium ions have sufficient time to diffuse evenly across the electrode surface, promoting the formation of a smooth, dense lithium layer. However, as current density increases, the rate of lithium-ion deposition surpasses the diffusion rate, resulting in localized areas of high lithium-ion concentration and subsequent dendritic growth. Furthermore, the concentration gradient of lithium ions in the electrolyte is a dynamic factor during battery operation. As lithium ions are plated onto the anode during charging, their concentration near the electrode surface decreases, potentially leading to uneven deposition if the ions in the electrolyte at the anode are not replenished quickly enough [
24].
Assessment of the development of dendritic growth in cells in series cell configuration reveals remarkable uniformity in growth rates. In a representative MRI slice, the volume of dendrites is approximately equal in all cells at each timepoint, with all dendritic structures reaching the opposite electrode at approximately the same time (after approximately 28 hours). This suggests the existence of some kind of self-regulatory mechanism in which the existence of a runaway cell is unlikely. In this configuration, the current through all the cells is the same, and so is the dendrite growth rate, which means that the amount of lithium deposited must be the same in all cells at all times. Note, that this can still not explain the uniformity of dendrite growth across the anode surface, as this is still a very uncontrollable phenomenon.
Since the nucleation rates are heavily dependent on the current density, any surface imperfections are expected to increase the probability of nucleation. This can be seen, for example, in the series cell configuration, second cell from above. Due to cell assembly imperfections, the middle electrode is slightly bent, resulting in an uneven distance between the electrodes in the second and third cells, which has as a consequence preferential growth at the point of shortest distance in the second cell, as well as a lack of growth at the same (horizontal) location in the third cell.
In the case of parallel cell configuration, all four separate cells are connected to the same source and have thus the same electric potential. This basically means that if there is a cell with a faster dendrite growth rate, that cell has a relatively lower resistance than the rest of the cells and thus a proportionally higher current. This further accelerates the growth of dendrites in this cell compared to the rest of the cells, resulting in a runaway cell. Another factor that indicates that this configuration is disadvantageous compared to the single cell, is the fact that ion diffusion is limited from the neighboring cell, which means that in theory, local concentration gradients should be much higher in the case when cells are narrower. Note, that in our experiments, the cells were not completely isolated from each other, therefore some amount of lateral ion diffusion is still permitted in the vicinity of electrodes.
Numerous studies have focused on the morphology of dendritic structures and found their typical scale to be in the range of 1–3 μm [
25]. In contrast, the voxel size used in our experiments was 187 μm. However, this limitation of MRI resolution is not the only limiting factor. Very important is also signal attenuation due to magnetic susceptibility and RF shielding effects. The extent of this attenuation depends not only on the dendrite volume to voxel volume ratio, but also on the shape and orientation of the dendrite with respect to the magnetic fields
B0 and
B1 [
19,
22]. As a result, images generated by the indirect imaging technique are limited in providing direct information about specific properties of dendrites, such as their density.
Graphs of the apparent dendrite volume per electrode area as a function of charging time in
Figure 4 offer several interesting insights. First, we can see that the apparent dendrite volume is up to ten times larger than the corresponding theoretical value of deposited lithium per electrode area, which is about 0.25 mm
3 per 1 mm
2 of electrode after 60 hours. Second, contrary to theoretical expectations, the growth of apparent dendritic volume is non-linear. This growth pattern indicates the presence of different growth regimes at different charging times, resulting in dendrites of different structures and thus different effects on the MR signal reduction. Thus, the growth of the apparent dendrite volume is gradual in denser dendrites, while it is faster in more branched dendritic structures. There is also a marked difference in the apparent dendrite volume growth among different battery configurations. Specifically, in the parallel cell configuration the growth is significantly more pronounced compared to the series cell configuration. This finding highlights the importance of cell configuration in influencing dendrite growth patterns in lithium-ion batteries.
The presence of gas pockets in all configurations and their random changes raise questions about chemical reactions or possible leaks in the battery case. The composition of the electrolyte must be stable in the operating voltage and temperature range [
26]. Although moisture levels were consistently below 1 ppm during assembly, LiPF
6 is sensitive to traces of moisture impurities in the electrolyte and will react to produce a small amount of hydrofluoric acid (HF) [
27]. It could also be that the sealing of the assembled cell was compromised, allowing moisture from the environment to enter the cell during the experiments. As for the effect of the formation of gas bubbles, they will disturb the dendrite formation process by disrupting the flow of electric current. Given the sensitive nature of lithium dendrites, it is likely that forces arising from surface tension may also influence the development and morphology of these dendritic structures. However, this phenomenon requires further investigation for a better understanding.
In the original publication on the indirect lithium dendrite imaging [
19], the authors used a FLASH imaging sequence to calculate the dendrite volume fraction from the image signal intensity. Given that FLASH is a gradient-echo sequence, it is more sensitive to
B0 inhomogeneities and thus to susceptibility effects, which were identified in this study as major contributors to signal attenuation. However, a different MRI sequence was used in this study, namely the RARE imaging sequence, which is based on spin-echo and is therefore less susceptible to magnetic susceptibility effects and thus to dendritic structures than the FLASH sequence. The sensitivity of the used imaging sequence to magnetic susceptibility and to RF shielding may affect the threshold settings for correct segmentation of dendrites or the ratio of MR signal to the dendrite volume fraction. Direct determination of lithium volume fraction from the normalized MR signal remains a possible better alternative to thresholding for our future studies. One such relation can already be discerned from
Figure 7 of this study; i.e., the relation between the
f-factor (lithium volume fraction) and the normalized signal
S. However, this relation only applies to the mean values of dendritic voxels, and not for all voxels and all of their values at any given time in the experiment. For example, at the beginning of the experiment dendrites are less developed and less dense, while latter in the experiment they are more developed and therefore denser. Due to this different structure of dendrites, their effect on reducing the MR signal is different; i.e., a lower amount of lithium in earlier dendrites can cause the same reduction in MR signal as a higher amount of lithium in later dendrites.