Zhi
Zhu
*ab,
Shuanglong
Xu
a,
Zhenjie
Wang
a,
Xiaohui
Yan
a,
Guiyin
Xu
b,
Yimeng
Huang
c,
Yuping
Wu
a,
Yin
Zhang
*b and
Ju
Li
*bc
aConfucius Energy-storage Lab, School of Energy and Environment & Z Energy-storage Center, Southeast University, Nanjing, Jiangsu 211189, China. E-mail: zhizhu@seu.edu.cn
bDepartment of Nuclear Science and Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA. E-mail: yinzhang@pku.edu.cn
cDepartment of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA. E-mail: liju@mit.edu
First published on 25th July 2024
Charging LiCoO2 (LCO) to above 4.5 V induces crystal cracking and seriously deteriorates the battery cycle life. Decreasing the range of the LCO misfit strain during deep de-lithiation is useful for preventing cracks, but this is not always achievable. Here, we demonstrate that the limited electrochemical contact area between electronically conductive carbon and the LCO crystal causes “electrochemical indentations” (ECIs) during charging and discharging. Particularly in fast charging, the high local ΔcLi gradient in LCO would cause a local volume of the surficial lattice to shrink while the rest of the crystal is still under stretching, and hence, drive the ECI to cause cracking. Increasing the electrochemical contact area would reduce the ECI and cracking risk. Therefore, we developed a free-standing CNT-LCO electrode in which all of the LCO particles were intimately wrapped with a dense CNT cocoon to establish a larger true electrical contact area. The simulations demonstrated that the radial ΔcLi and ECI decreased significantly in the cocooned LCO particles. The cocooned LCO electrode maintained good morphology and retained 94% of its energy density after 400 cycles when charged to 4.55 V. By removing the need for a current collector and binder, the volumetric energy density of the CNT-LCO cathode reached 3200 Wh L−1 (electrode).
Broader contextUtilizing the extra capacity of the LiCoO2 (LCO) cathode above 4.5 V is enticing, but it is impeded by crystal cracking and rapid cycling decay for unknown reasons. This work demonstrated that the limited true contact area between LCO and electronic conductive agents is the factor that causes “electrochemical indentations” (ECIs) during high-voltage charging and triggers crystal cracking. We develop a numerical model based on electrochemomechanics and show that avoiding ECI is key to preventing LCO cracking and improving the high-voltage cycling life. Subsequently, a free-standing carbon nanotube-cocooned LCO electrode (CNT-LCO) was developed, where only a 0.5 wt% CNT weaved dense elastic cloth fully wrapped all of the LCO particles, like silkworm cocoons, to maximize the electronic contact area. The CNT-LCO retained 94% of its energy density after 400 cycles when charged to 4.55 V, in contrast to the catastrophic capacity decay using 3 wt% carbon black, with everything else remaining the same. The theory of electrochemical indentation and how a low electronic contact area leads to ECI and cracking are essential for all mixed ionic electronic conductors (MIECs) and even for solid electrolytes, which is a limiting case of MIECs. |
In previous works, it has been demonstrated that charging LCO to above 4.5 V involved oxygen redox and triggered particle cracking, which caused the reversible capacity of LCO to decay quickly.8,10,11 Crack formation in LCO was understood from the repeated lattice-volume changes during cycling.12 Thus, many material modifications were conducted on LCO, e.g., doping and coatings,4,8,13,14 to inhibit the phase transformation and decrease the amplitude of the LCO-lattice volume change (ΔV = Vmax − Vmin) during the deep de-lithiation process. However, based on the experimental observation that the crack always originates from the specific sites at the LCO particle surface,6,10 we should be aware that the dynamical spatial gradients of Li-ion concentration in the large LCO single crystals can be a more critical factor for crack formation.
In common cathode electrodes, it was believed that the synergic Li-ion diffusivity and electronic conductivity in the active materials (σ1) and the electronic conductivity (σ2) in carbon black (CB) determined the kinetics of the cathode reaction.15–18 However, as shown in Fig. 1a, the lean CB (<3 wt%) conductive agent that was used in the regular LCO electrode (R-LCO) for high energy density only had a narrow and limited electrical contact area (aA) with the LCO particles, so the local current density (iA ≈ I/aA) would concentrate at the electrical contacting sites (site A), as shown in Fig. 1b. Since the conduction of electrons and Li ions is always coupled in the LCO particles, the high current density resulting in an overpotential (ΔϕA) brings in a high gradient of Li concentration (ΔcLi), causing inhomogeneous lattice-parameter change and dynamic stress inside the single crystal during charging.
Fig. 1 Formation and prevention of electrochemical indentation in LCO particles during high-voltage charging. (a) SEM images of LCO particles and carbon black in R-LCO. (b) Scheme of electrochemical indentation and crack formation in LCO crystals when charging R-LCO to a high voltage. (c) SEM images of CNT-cocooned LCO particles in CNT-LCO. (d) Scheme for preventing electrochemical indentation and crack formation in the LCO crystal when charging CNT-LCO to a high voltage. The fundamental equations for mixed ionic electronic conductor (MIEC) are presented in ref. 27. |
This dynamic stress could be bearable when the LCO crystal was charged to an upper voltage of less than 4.5 V.13 However, when charging the LCO beyond 4.5 V at a high rate, the over-delithiated surface (csurface) at site A would first get oxygen redox (O2− → Oα−, α < 2) involved and rapidly tend to shrink the local lattice (LicsurfaceCoO2), while the less delithiated bulk (Li(csurface + ΔcLi)CoO2) is still under stretching,13 then the high shear stress near site A would act like an electrochemical indentation (ECI) at the crystal surface by the limited-area carbon black, eventually triggering cracking (Fig. 1b). Moreover, the cracks at site A would continue decreasing the electrical contact between the LCO particle and CB (decreasing aA and increasing ΔcLi) and aggravate ECI with continuous crack growth, until the LCO particles lost electrical contact with the CB and eventually terminated the high-voltage battery cycling. Note that we use c in this paper to denote alternatively concentration (unit #Li per m3) or mole fraction (#Li per formula unit, which is dimensionless).
Eliminating ECI could be the key to preventing LCO crystals from cracking; thus, decreasing the dynamical radial overpotential (Δϕ) during the high voltage de-lithiation process is important to improve the cycling stability of LCO. As the electronic contact between LCO and CB is significant for Δϕ, compared with singly improving σ1 with material doping or singly improving σ2 with a new conductive agent, increasing the electrical contact area for the LCO particles to disperse the current density (i) can be more effective. In this work, we established a CNT cocoon cloth to amply wrap the LCO particles in the electrode, named CNT-LCO (Fig. 1c). Compared to R-LCO, CNT-LCO has the following attributes:
(1) Due to the high aspect ratio of CNT (Fig. S1 and Table S1, ESI†), a small weight percentage of CNT can weave dense elastic cocoon cloth for the LCO particles, like silkworm cocoons, which would not only entwine and cover all LCO particles but also result in thorough electron percolation in the electrode, as shown in Fig. 1c. Such intimate contact between the LCO particles and CNT cocoon will maximize the true electronic contact area (aB → asurface, where asurface is the total surface area of all LCO particles) to decrease iB, ΔϕB and ΔcLi (Fig. 1d). Then, the ECI and cracking in the LCO particles can be greatly prevented during high-voltage cycling, as revealed by our finite-element diffusion-stress modeling of the “contact electromechanics” of mixed ionic electronic conducting (MIEC) LCO with a conductive agent.
(2) CNT can form a complete matrix to establish a robust free-standing electrode without any chemical binder or current collector. Thus, the free-standing CNT-LCO electrode has a high volumetric energy density of 3200 W h L−1 (electrode), which would, in turn, ideally improve the energy density of the LCO full cell from the current 780 W h L−1 (cell) to above 1212 W h L−1 (cell).
Though carbon nanotubes or nanobelts have been used in previous works for LCO electrode preparations,19–21 the mechanisms are not fully understood. Creating a CNT network to enhance the electrical conductivity is far from sufficient to stabilize high-voltage LCO cycling. As we will show, the key is to maximize the true electrical contact area to eliminate the electrochemical indentation. In this work, finite-element simulations and experimental analysis indicated that the plane-to-plane contact between the LCO particles and CNT cocoon efficiently eliminated the electrochemical indentation during high-voltage charging so that CNT-LCO had well-maintained morphologies without cracks and retained 94% of its energy density after 400 cycles when charged to 4.55 V. In addition, the CNT-LCO cathode had a high volumetric energy density of 3200 W h L−1, and the cycled CNT-LCO could be facilely regenerated/recycled without any complicated mechanical or chemical treatment that must be involved in the traditional battery material regeneration/recycling processes.
The R-LCO and CNT-LCO electrodes were first tested in half cells at 200 mA g−1. The charge/discharge profiles in Fig. 2a indicate that while R-LCO was charged to 198.1 mA h g−1 and discharged to 195.7 mA h g−1 in the 10th cycle between 3.0 and 4.55 V, CNT-LCO was able to be charged to 200.5 mA h g−1 and discharged to 199.1 mA h g−1. After 300 cycles, the charge/discharge profiles of R-LCO were seriously deformed, as shown in Fig. 2a, which could only be charged to 99.3 mA h g−1 and discharged to 97.8 mA h g−1. In contrast, the voltage profile of CNT-LCO remained very stable, and it could still be charged to 192.7 mA h g−1 and discharged to 192.3 mA h g−1. It is quite amazing that without any chemical binder or extra conductive agent and no changes to the active material LCO itself, using just 0.5 wt% CNT (which may be considered both a physical binder and a conductive agent) versus the standard method of electrode preparation can lead to not only better initial performance but also much better cycling durability. This illustrates that the hypothesis of a low electronic contact area leading to ECI is a very important one to check out, and we aim to develop a theory about this later in this paper.
The cycling performances of R-LCO and CNT-LCO at 200 mA g−1 between 3.0 and 4.55 V in half-cells are compared in Fig. 2b. While the discharge energy density of R-LCO decreased to 210 mW h g−1 after 400 cycles, that of CNT-LCO remained as high as 757 mW h g−1, 94.0% of that of the 10th cycle (0.015% decay per cycle). Additionally, while the average polarization resistance (Rp ≡ (Vc − Vd)/2I) of R-LCO increased to above 120 Ω in the 400th cycle, that of CNT-LCO decreased to 21 Ω. The suppressed R indicated that the CNT-LCO cathode had well-maintained electrochemical kinetics during high voltage cycling, which not only promoted the high retention of discharge energy density but also reduced heat generation during fast charging/discharging.
We also note that the Coulombic inefficiency (CI ≡ (100 − CE)%) of CNT-LCO was as low as 0.16%, while that of R-LCO increased to 1.3% after 400 cycles. The CI increase during the R-LCO cycling is likely a result of particle cracking during high-voltage charging, which would not only aggravate the irreversible electrolyte decomposition at the newly generated crack surfaces and create new cathode–electrolyte interphases (CEI) but also impede electronic/polaronic conduction in the following cycles. In contrast, the steady CI during the CNT-LCO cycling indicates that the CNT cocoons greatly suppressed cracking and irreversible side reactions from the electrolyte at high voltage. The cycling performance at higher rates is shown in Fig. S4 (ESI†). Fig. S5 (ESI†) shows that CNT-LCO retained 83.6% and 75.8% of its capacity after 280 cycles when charged to higher upper voltages of 4.60 V and 4.62 V.
To better understand the mechanism for the different electrochemical behaviors during R-LCO and CNT-LCO cycling, GITT tests were performed on the two cathodes. As shown in Fig. 2c and Table S2 (ESI†), from the ε that was mainly contributed by the Ohmic resistance (IRΩ) during the 10th charging, it can be inferred that the electronic resistance in the 0.5 wt% CNT-cocooned CNT-LCO electrode was slightly lower than that of the 3 wt% CB-added R-LCO electrode in the initial cycles. However, while εR-LCO increased by about three times after 300 cycles, εCNT-LCO remained almost the same as in the 10th cycle. This indicates that the matrix for electrical conduction in R-LCO may have been broken during high-voltage cycling, but the CNT network in CNT-LCO maintains long-range electronic percolation throughout the electrode.
More importantly, the GITT results showed that , which was attributed to the polaron diffusion from the interface to the bulk of the LCO crystals, was surprisingly different between R-LCO and CNT-LCO. Fig. 2c and d and Table S2 (ESI†) shows that the of CNT-LCO (e.g. at an OCV of 4.0 V) was much lower than that of R-LCO in the 10th cycle; moreover, the of R-LCO increased by more than 3 times after 300 cycles, whereas that of CNT-LCO changed very little. For cathode crystals, σ1 can be calculated22,23 using eqn (1):
(1) |
(2) |
eqn (1) can be rewritten as
(3) |
Considering that R-LCO and CNT-LCO have the same LCO particles and should have the same σ1, the value of can be proportioned to L in the particle; therefore, the lower in CNT-LCO indicates that the L in CNT-LCO is much shorter than that in R-LCO. In other words, the full contact between the CNT cocoon and the LCO particles has greatly optimized the polaron diffusion path in LCO and shortened L in the large LCO crystals during cycling. From the GITT results in Fig. 2c–f and Table S2 (ESI†), the L in CNT-LCO can be estimated as 1/1.13 to that in R-LCO in the 10th cycle, while L significantly increased in R-LCO after 300 cycles, it remained much stable in CNT-LCO and could be estimated as less than 1/3.17 of the former.
The shortened L in CNT-LCO can be a result of the maximized electrical contact interface between the LCO particles and the CNT cocoon (aB → asurface), which would erase the huge Li concentration gradient in LCO and eliminate the ECI. To comprehensively investigate the phenomenon of electrochemical indentation of the MIEC under fast high-voltage charging, we conducted finite-element simulations to accurately model the de-lithiation process in the LCO crystals when charging R-LCO and CNT-LCO. To model the LCO particles in R-LCO, the lithium fluxes were imposed on selected twelve rectangular patches with a size of 9 μm2 each. In contrast, the LCO particles in CNT-LCO experienced a uniform lithium flux over the whole surface. The total current contributed by each LCO particle (I) was maintained at 200 mA g−1 for both electrodes.
As shown in Fig. 3a, the LCO particles in R-LCO exhibit a larger radial lithium concentration gradient compared to that in CNT-LCO during charging. Specifically, when the lithium mole fraction dropped to 0.55 near the electric contact area with CB (site A), the lattice at the LCO surface would trigger the phase transition from a hexagonal structure to a monoclinic structure. This low Li concentration at the surface can lead to an irreversible structural transformation that significantly affects battery life. Since the de-lithiation process introduces a volume change of around 2.5%,13 the nonuniform de-lithiation resulting from a small contact area produces a high Li gradient and nonuniform volume change. This incompatibility would lead to an extremely high local mechanical stress of up to 2000 MPa, as shown in Fig. 3b, which approaches the yield strength of the LCO crystal at around 2700 MPa.
ECI is similar to mechanical indentation,24 where a small mechanical load is applied to the materials through a sharp indenter. Here, the partial contact produces a nonuniform Li concentration and volume change, which further leads to a localized stress field near the contact area. Therefore, we can call this phenomenon electrochemical indentation, which is a contact electromechanics problem analogous to the contact problems in mechanics.25 Electrochemical indentation imposes a high local mechanical stress near the surface of LCO particles, which could lead to early cracks and degrade the LCO particles in R-LCO. However, Fig. 3c shows that the LCO particles in CNT-LCO have a much smoother and spherically symmetric Li concentration field because of the fully wrapped CNT cocoon on the LCO particle surface. Note that, under the same charge/discharge current, the ΔcLi between the surface and centroid within the LCO particle in CNT-LCO was only 7% of the initial concentration, while that of R-LCO reached 45%. Due to the relatively low Li concentration gradient, the mechanical interior stress within the LCO particles in CNT-LCO was much lower, as shown in Fig. 3d, and was only around 900 MPa at the outer surface.
More importantly, upon further charging, the deep de-lithiation would involve oxygen redox and cause the lattice at the LCO surface to shrink, which is contrary to the initial volume expansion in bulk. As shown in Fig. 3e, this abnormal shrinkage would apply an additional load to the surface and lead to the inelastic deformation of the LCO particles in R-LCO at around 100s of charging. On the other hand, the LCO particles in CNT-LCO remained elastic and undamaged until 1500s of charging, when the abnormal shrinkage flipped the sign of the pressure field near the surface from positive (compression) to negative (tension). Therefore, the time to the initial damage (td) in the LCO particles was significantly prolonged from 100s in R-LCO to 1500s in CNT-LCO after the introduction of the CNT-cocoon lamination. The finite element simulations demonstrated that an enlarged electronic contact area between the LCO particles and CNT cocoon is greatly favorable for erasing the high local Li concentration gradient and eliminating the electrochemical indentation in the LCO particles during high-voltage cycling.
In Fig. 3f, sharp pressure spikes and transitions from positive to negative pressure are illustrated, corresponding to the non-linear volumetric changes depicted in Fig. S6 (ESI†). During the de-lithiation phase, the unit cell volume of LixCoO2 initially increases from 96.50 Å3 to 98.90 Å3 as the lithium content x approaches 0.6. This volumetric expansion is primarily due to the increased electrostatic repulsion among the CoO2 layers as lithium ions are progressively removed. As de-lithiation continues beyond this point, the unit cell volume subsequently contracts to 95.31 Å3. This contraction is linked to the repositioning of the CoO2 layers in response to the significant lithium removal, resulting in a denser structural configuration. The abrupt change in the partial volume of lithium at the critical x of 0.6 leads to these sharp pressure spikes and the transition from a compressive to a tensile load.
ECI during high-voltage charging would damage the LCO crystal and cause particle cracking; therefore, the elimination of the ECI can be key to maintaining the integrity of LCO particles during high-voltage cycling. Fig. 4a and b compare the morphologies of the LCO particles in R-LCO and CNT-LCO, respectively, after 120 cycles. The SEM images show that in the cycled CNT-LCO electrode, the CNT cocoon still fully covered and tightly wrapped all of the LCO particles, maintaining intimate electrical contact between the LCO particles and the CNT cocoon, whereas much of the LCO particles lost electrical contact with the CB in R-LCO. More importantly, while severe corrosion and deep cracks were found everywhere on the LCO surface in R-LCO, the LCO particles in CNT-LCO still had smooth surfaces.
Fig. 4 Particles after cycling and full-cell cycling performance. (a) and (b) SEM images of the LCO particles in R-LCO and CNT-LCO after high-voltage cycling. (c) TOF-SIMS analysis of the surface of the cycled R-LCO11 and CNT-LCO electrodes. (d) Cycling performance of discharge energy retention and cell resistances of R-LCO and CNT-LCO full cells. (e) and (f) EDS mapping of Co on the graphite anode from the cycled R-LCO (e) and CNT-LCO (f) full cells. |
The cracks at the LCO particle surface would introduce and expose new particle surfaces to the electrolyte and increase the interface area between the cathode and electrolyte. The increase of the free surface area would seriously aggravate the side reactions with liquid electrolytes at high voltages, including the oxygen release and electrolyte decomposition, so that the maintained morphology of LCO particles in CNT-LCO would suppress these side reactions and prolong high-voltage cycling life. To better compare the side reactions that occurred on the surfaces of the two cathodes, we performed TOF-SIMS to analyze the surficial compositions of R-LCO and CNT-LCO after cycling. Fig. 4c and Fig. S7 (ESI†) show the distribution of the CoF3− component at the surface of the two cathodes, which is a representative indicator of cathode corrosion from HF attack in the electrolyte. While the CoF3− layer was thick at the surface of R-LCO and in a penetrating fashion towards the deep,11 it was much thinner at the CNT-LCO surface. In other words, the surface of R-LCO was seriously corroded, and the damage spread into the particle bulk (along the cracks in the crystal) after high-voltage cycling, but such corrosion was very weak and only occurred at the outermost surface of the CNT-LCO.
Fig. 4c and Fig. S7 (ESI†) also compare the distribution of C2HO− species to evaluate electrolyte decomposition. It was quite clear that the C2HO− layer on the CNT-LCO surface was also much thinner than that of R-LCO,11 indicating that the electrolyte decomposition on the CNT-LCO surface was greatly prevented during cycling. The suppressed surface corrosion and electrolyte decomposition at the CNT-LCO surface confirmed the prevention of LCO cracks.
We then tested R-LCO and CNT-LCO in pouch full cells that were paired with commercial graphite anodes and a lean electrolyte (2 g A h−1) at a current density of 200 mA g−1. As shown in Fig. 4d, the discharge capacity retention of the CNT-LCO full cell was 91% after 300 cycles, whereas that of R-LCO degraded to 37%. The interior resistances were also quite different between the two cells: that of R-LCO increased to 1.65 Ω after 300 cycles, while it was only 0.34 Ω for the CNT-LCO cell. The stable cycling stability of the CNT-LCO full cell not only verified the validity of the CNT cocoon in stabilizing high-voltage LCO cycling, but it also indicated that the cocooning technique has practical application potential. After cycling, we also measured the Co distribution at the surface of the graphite anodes that disassembled from the cycled full cells. The EDS mapping in Fig. 4e and f shows that the Co deposition on the graphite anode from the CNT-LCO cell was much less than that from the R-LCO cell. It can be inferred that the prevention of LCO cracks also greatly suppressed Co dissolution from the cathode, which further maintained the interior resistances and stabilized high-voltage full-cell cycling.
In the CNT-LCO electrode, 0.5 wt% CNT was not only fully wrapped around all LCO particles but also self-intertwined to be a free-standing electrode so that no additional binder or metal current collector was necessary for preparing the electrode, as shown in Fig. 5a. In this work, we compressed this free-standing electrode to 4.45 g cm−3 (Fig. 5b), then the volumetric energy density of CNT-LCO reached ∼3200 W h L−1 (electrode) based on the whole free-standing cathode volume, and the volumetric and gravimetric energy densities of the new designed full cell can be expected to be significantly improved to 1212.5 W h L−1 (cell) and 341 W h kg−1 (cell) (Fig. S8 and the Discussion section, ESI†), respectively.
Besides the improvements in high-voltage cycling stability and energy density, the free-standing CNT-LCO cathode also makes it easier to recycle at the end of life. Owing to the absence of an aluminum current collector and sticky polymeric binder, the cycled CNT-LCO electrode can be efficiently regenerated and recycled with a facile in situ air annealing treatment at 500–600 °C, which oxidizes the CNT-cocoon and organic species at the surface. This does not require the complicated separation and collection processes of the active materials, which must be involved in the regular battery regeneration and recycling process.7,26 The thermogravimetric (TG) analysis of the regeneration process is shown in Fig. S9 (ESI†). The recycled LCO materials are still robust and can be again directly used for a new free-standing CNT-LCO electrode (inset in Fig. 5c) in the cell.
Fig. 5c and d show the charge/discharge profile and cycling performance of the regenerated CNT-LCO cathode between 3.0 and 4.55 V in half-cells. Though the reversible capacity of the regenerated CNT-LCO was slightly lower than that of the fresh electrode (Fig. 2a), it was surprising that the cycling stability even improved, as the decay of energy density was less than 0.013% per cycle after the 150th cycle. The improved cycling stability of the regenerated CNT-LCO can be understood from the appropriate lattice doping of F at the LCO particle surface during annealing (Fig. 5d), where the F element came from the CEI that formed in the previous long cycling. The F doping at the regenerated LCO surface is demonstrated by the XPS analysis, as shown in Fig. S10 (ESI†). An appropriate amount of F doping at the LCO particle surface could greatly prevent lattice-oxygen escape and HF corrosion from the electrolyte, which both favored high-voltage LCO cycling stability. This result indicates that free-standing CNT-LCO can be a promising green battery candidate.
This CNT-cocooned free-standing electrode technology can also be generalized to other electrode materials, such as high-Ni NMC cathodes, LiFePO4 cathodes, graphite anodes, and Si anodes. The CNT cocoon will not only improve the cycling stability of the active materials (especially under fast charging) but also improve the energy density of the cell. In addition to the simple, environmentally friendly and cost-friendly recycling and regeneration process, the CNT-cocooned free-standing electrodes provide important guidance for new sustainable electrode design and battery fabrication. To promote this free-standing electrode technology in larger pouch cells, new electrode structures and fabrication technologies are being developed.
Lastly, we believe that the theory of electrochemical indentation (ECI) and how a low electronic contact area leads to ECI-induced cracking is fundamentally important for all mixed ionic electronic conductors (MIEC) and even for solid electrolytes (SE), which is just a limiting case of the MIEC (electronic transference number telectronic < 10−4).27 With no changes to the active material, using just 0.5 wt% CNT (serving both as a physical binder and a conductive agent) instead of the standard electrode preparation led to better battery performance and durability. This very surprising demonstration and the underlying mechanism should be realized by the broader electrochemical material community.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4ee00722k |
This journal is © The Royal Society of Chemistry 2024 |