Mahla Sarfaraz Khabbaza,
Anupma Thakurb,
Dipak Maity
c,
Sepideh Biabanialitappeha,
Bhargava Morthaa,
Joshua P. McKaya,
Gerardo Paradaa,
Mangilal Agarwalc,
Honghu Zhang
d,
Xin Zhang
e,
J. David Bazak
e and
Xiaoliang Wei
*ab
aSchool of Mechanical Engineering, Purdue University, 723 West Michigan Street, Indianapolis, IN 46202, USA. E-mail: wei304@purdue.edu
bSchool of Materials Engineering, Purdue University, 701 West Stadium Avenue, West Lafayette, IN 47907, USA
cIntegrated Nanosystems Development Institute, Indiana University Indianapolis, 625 West Michigan Street, Indianapolis, IN 46202, USA
dNational Synchrotron Light Source-II, Brookhaven National Laboratory, Upton, New York 11973, USA
ePhysical and Computational Sciences Directorate, Pacific Northwest National Laboratory, 3340 Stevens Drive, Richland, WA 99354, USA
First published on 23rd July 2025
The inexpensive sulfur raw material is promising to enable cost-effective redox flow batteries for long duration energy storage. But the catastrophic through-membrane crossover of polysulfides remains a severe challenge resulting in irreversible performance degradation and short cycle life. In this work, we demonstrate that use of a permselective cation exchange membrane yields a two orders of magnitude enhancement in polysulfide retention compared to the benchmark Nafion membrane. Combined physico-chemical, spectroscopic, and microscopic analyses suggest more disordered sidechain structures, which lead to the more hydrophobic nature and smaller hydrophilic domains in the membrane. The microstructural features contribute to the effective mitigation of polysulfide crossover. As a result, the cycle life of polysulfide/ferricyanide flow cells is boosted over a substantially extended test time. This finding sheds light on the fundamental membrane factors that cause polysulfide permeation and can provide feasible directions in the development of permselective membranes for polysulfide flow batteries.
Elemental sulfur is an inexpensive material sourced mainly from petroleum refining and natural gas processing. Sulfides including monosulfide and short-chain polysulfides (Sx2−, x = 1–4) are among the classical inorganic anolyte redox materials used in aqueous RFBs.10–12 These sulfides generally have high solubilities (>2 M) with two transferred electrons and are highly promising for achieving high energy density in their RFBs. However, sulfide-based RFBs have limited success in commercialization. As the only example with commercial-scale (12 MW/120 MWh) installations, Innogy's polysulfide/bromine flow battery (PSB) projects were all terminated before commissioning.13 One of the near-fatal drawbacks is the irreversible crossover of sulfide species. Even the benchmark cation exchange membranes (CEMs) such as Nafion have failed to prevent sulfide crossover to the catholyte side, resulting in sulfide loss from the anolyte, solid sulfur precipitation at the catholyte, and eventually performance failure. With single or directly stacked Nafion membranes, stable capacity could be achieved for short cycling only and long-duration operation was often unsuccessful.14–18 To address this issue, advanced membrane approaches have been investigated to increase permselectivity against sulfides. Lu et al. developed a charge-reinforced ion-selective (CRIS) membrane formed with porous carbon coatings on a Nafion membrane.19 This CRIS membrane took advantage of the strong sulfide adsorption to form an electro-repulsive permeation barrier to retain sulfides, which enabled a KI/K2S2 flow cell to demonstrate stable cycling for ∼3 months. Chen et al. adopted a similar strategy using an electrocatalyst embedded membrane electrode assembly (MEA) to achieve extended stable cycling for a Na2S2/Na4Fe(CN)6 flow cell.20,21 A dual membrane configuration with parallel aligned CEM and AEM was used in an alkaline polysulfide-air flow cell, but seriously fluctuating efficiency and an increased cell overpotential were observed.22 Chiang et al. employed a Li+-exclusive superionic conducting (LiSICON) ceramic membrane in an air-breathing polysulfide RFB system,23 but the membrane suffers from a critical drawback of limited scalability because of low conductivity, poor mechanical strength, and high cost. Considering these drawbacks, the development of advanced sulfide-blocking membranes is an urgently needed area for enabling long cycling polysulfide-based RFBs.
In this work, we report extended stable cycling for an aqueous near-neutral RFB system using a Na2S2 anolyte and K3Fe(CN)6 catholyte (denoted as S/Fe), enabled by a cost-effective Fumasep F1850 membrane with enhanced permselectivity over polysulfides. The ferro/ferricyanide redox couple is one of the most commonly used catholyte materials for aqueous RFBs with well-defined electrochemical properties.24,25 Combined physico-chemical, spectroscopic, and microscopic characterizations of this membrane indicate higher hydrophobicity and smaller sizes of hydrophilic regions, potentially imparted by the greater extent of structural disorder in the sidearm compared to the benchmark Nafion. This microstructural feature contributes to the favored improvement in polysulfide retention by two orders of magnitude, better than that of Nafion. The S/Fe flow cells demonstrated near-constant capacity for as long as 300 cycles, which corresponds to 49 days, with the fade rate as low as 0.176% per day.
The setup consisted of a 0.5 M source solution (15 mL) on one side and inert permeate solution (15 mL) on the other side, sandwiching the tested membrane. For Na2S2, 1 M NaCl was used as the permeate solution to balance ionic strength and minimize osmotic pressure. For K3Fe(CN)6, 1.5 M KCl was used as the permeate. The effective area of the disk-shaped membrane was 0.78 cm2. Both sides of the H cell were magnetically stirred. At specified time intervals, 3 mL samples were taken from both sides and the permeate side samples were subjected to ultraviolet-visible (UV-vis) measurements using a Thermo Scientific Evolution 600 spectrometer. Following Beer's law, their absorbance peak intensities were interpolated to the pre-determined intensity-concentration calibration plot to obtain the evolving crossover concentrations. The through-membrane permeability, P, was determined according to eqn (1) below:19
![]() | (1) |
![]() | (2) |
![]() | (3) |
![]() | (4) |
MAS-NMR was primarily conducted at 30 kHz, with additional spinning rates of 15 kHz and 20 kHz used to confirm the isotropic peak assignments. 19F spectra were acquired with a Hahn echo, with the net echo delay (194 μs) set to achieve a refocusing time of four rotational periods. Whole-echo acquisition was used with a brief dead time (6.5 μs) and the resulting free-induction decay was left-shifted to the refocusing point using ssNake v1.4 (which was also employed for peak deconvolution, including all spinning sidebands).27 All echoes were recorded with 256 scans, a spectral width of 200 kHz, an acquisition time of 10.2 ms, and a recycle delay of 30 s. Pseudo-Voigt line shapes were used for all deconvolutions, with both isotropic shift and Lorentzian/Gaussian linewidth components for all peaks regressed against the experimental data. Signal fractions were based on addition of all associated spinning side bands to the corresponding isotropic peak areas. 19F–13C cross-polarization spectra were also acquired, with the Hartmann–Hahn match condition determined empirically using N115 and a contact time of 1.5 ms, utilizing rotor-synchronized π-pulse decoupling on the 19F channel during the acquisition period.28
Kinetic measurement of K3Fe(CN)6 was performed using its CV curves at different potential sweeping rates v ranging from 10 to 100 mV s−1. The diffusion coefficient (D) was determined via the Randles–Sevcik equation (eqn (5)):29
ip = 0.4463 (nF)3/2AC(vD/RT)1/2 | (5) |
Ψ = (−0.6288 + 0.0021ΔEp)/(1 − 0.017ΔEp) | (6) |
Ψ = k0(πDnF/RT)−1/2v−1/2 | (7) |
The bolted S/Fe flow cells consisted of graphite plate enclosures, thermally treated GFD2.5 at the catholyte, Co/GFD2.5 at the anolyte, and a membrane in between with an active area of 6.6 cm2. The K+ form membrane was soaked in deionized water for 24 h before cell assembly. The catholyte was 10 mL of 1 M K3Fe(CN)6 and the anolyte was 15 mL of 1 M Na2S2. Based on the loaded electrolytes, the polysulfide utilization was low corresponding to only 33% even at 100% state of charge (SOC) of the flow cells. The electrolytes were circulated at a flow rate of 20 mL min−1 using a Masterflex® L/S® peristaltic pump (Cole-Parmer, Vernon Hills, IL). Electrochemical impedance spectroscopy (EIS) was used to measure the Nyquist plots and derive the area specific resistance (ASR) of the flow cells with a Solartron SI 1287A electrochemical interface with an SI 1250E frequency response analyzer (Scribner Associates, Southern Pines, NC). A constant current mode with voltage control was used to cycle the flow cells at room temperature with a Neware CT-4008 battery tester (Belleville, IL).
To evaluate its effectiveness in improving the S/Fe RFB, the permeation rates of Na2S2 and K3Fe(CN)6 through this membrane were measured using H-type permeation cells (Fig. 1a). The evolving crossover concentrations in the permeate side were quantified by UV-vis spectroscopy. The UV-vis spectra as a function of concentration and the obtained peak intensity-concentration calibration curves for the two redox materials are shown in Fig. S1a, b, S2a and b in the ESI.† Fig. 1b and c show the measured UV-vis spectra of Na2S2 permeating through N115 and F1850 membranes, respectively. It is noted that there are two major absorbance peaks in the UV-vis spectra of dilute Na2S2 solutions. The dominant one at 300 nm belongs to S22− and the minor one at 372 nm can be assigned to S42−. This agrees well with the previous UV-vis reports on species speciation in polysulfide solutions.16,19 The main reason for the coexistence of multiple polysulfide ions is their similar Gibbs free energies resulting in complicated chemical equilibria of mixed disproportionation/comproportionation reactions.36,37 The electrolyte composition is concentration dependent and other polysulfide species such as S32− were also detected using Raman spectroscopy in more concentrated polysulfide solutions.20,21
As plotted in Fig. 1d, the F1850 membrane demonstrated a significantly improved polysulfide blocking ability compared to N115. The permeability of S22− through F1850 was determined to be 3.15 × 10−6 cm2 per day (according to Fig. S2c and d†), which is a nearly 100-fold reduction compared to that through N115 (2.12 × 10−4 cm2 per day). These values are comparable to, or even better than, those of the newly developed CRIS and semi-fluorinated polymer membranes.19,38 In addition, both membranes showed near-zero crossover for the Fe(CN)63− catholyte anion (Fig. S1c†), presumably because of its bulky size and high charge density creating a strong Donnan exclusion effect with the membranes. Therefore, the polysulfide permeation will be the dictating factor that controls the stability of S/Fe flow cells and enhanced cycling performance is expected for F1850.
The amount of water adsorbed by the membrane and the extent of dimensional swelling upon hydration reflect the size of water clusters. As shown in Table 1, the K+-form F1850 exhibits a water uptake (WU) of 3.9% and a swelling ratio (SR) of 4.1%, both of which are substantially lower than those of N115, indicating a lower volumetric presence of hydrophilic domains. As a result, cation transport through F1850 encounters a higher areal resistivity (Fig. S3†) and so does the polysulfide crossover (Fig. 1). Fig. 2a indicates the measured contact angles of 90° and 54° for K+-form F1850 and N115, respectively. Thus, F1850 exhibits an inferior surface wettability by water due to the higher hydrophobicity, which helps maintain low dimensional expansion upon hydration. FE-SEM images of dry membranes shown in Fig. 2b and c exhibit substantially different morphologies. N115 has a uniform, smooth surface, while F1850 exhibits a rough texture with randomly distributed circular bumps of ca. 2–10 microns. However, the formation mechanism for these bumps remains unclear but is speculated to be associated with the phase separation behaviors in F1850's microstructure. Due to the large scan areas, FE-SEM is unable to image nm-sized ion clusters even with an expanded view (Fig. 2c inset). Taking advantage of the good spatial resolution, AFM topography was used to analyze the membrane surface morphologies over smaller scan areas (Fig. 2d and e). A direct observation is that F1850 still has a higher surface roughness than N115 within the much smaller areal dimensions, which agrees with the FE-SEM analysis. The dark and bright colors are considered to represent the hydrophilic and hydrophobic domains, respectively.40 In N117, the hydrophilic regions are well interconnected aggregating into continuous ion channels of relatively large sizes. In contrast, in F1850, the hydrophilic domains are more substantially separated from each other, resulting in reduced interconnections and smaller dark regions. The surface microphase characteristics can explain the pronounced differences in ionic conductivity and polysulfide permeation between these two membranes.
To obtain fundamental insights into the superior polysulfide retention of F1850, solid-state MAS-NMR was used to characterize the microstructure of this membrane. Fig. 3 shows the 19F NMR spectra as well as the computed deconvolutions with chemical shift assignments to the various C–F groups in polymer structure (inset; see also Tables S1 and S2†). The nature and number of the peaks, including the linewidths, for N115 are broadly consistent with past, detailed NMR studies of Nafion polymer structure,26,28 with primary resonances for the {CF2}n peaks (overlapped with the S–CF2–C signal), the O–CF2–C signal (overlapped with the pendant –CF3 signal), and resonances for the –CF– signals at comparatively more negative ppm values (separate signals occur for the backbone and sidechain cases). Conversely, the F1850 structure exhibits a far greater variety of signals, particularly for the shift region associated with the signals from the pendant ionomer sidechains (which are deshielded relative to the backbone signals owing to the proximity to the ether linkages and the terminal sulfonyl groups). Peak assignments are made by a combination of contrasting with prior literature on Nafion's 19F MAS-NMR, comparing across MAS rates to identify spinning sidebands (Fig. S4†), and comparison of linewidths, which derive from a combination of conformational disorder and local dynamical mobility (i.e. 19F peaks of –CF3 rotors are narrow; 19F peaks near the bulky sulfonyl groups are comparatively broader). A similar proliferation of the observed signals for F1850 relative to N115 also occurs from the 19F–13C cross-polarization (CP) spectra (Fig. S5†).
![]() | ||
Fig. 3 19F MAS-NMR: (a) N115 and (b) F1850. Insets give the structural assignment of N115 following ref. 28 and the spectra at full zoom without truncation of the {CF2}n peaks. Deconvolution components are color-coded to the assignments, as given in the legends. Note for F1850 the appearance of additional signals in the sidechain –CF2– region (cyan), with broader linewidths consistent with the identified sidechain –CF2– resonances for N115. |
When contrasting the backbone and sidechain signals for N115, the signals attributed to the former yield a total fraction of 0.75 of the total 19F signal, including sidebands (Sbackbone = S−123 + (4/6)S−118 + S−139), versus 0.25 for the sidechains (Ssidechain = S−81 + (2/6)S−118 + S−145). In contrast, applying the same calculation for assignable signals in F1850 yields Sbackbone = 0.51 and Ssidechain = 0.21, leaving a considerable unassigned fraction, which is also difficult to directly link to the putative sidechain structure of F1850, whether a high degree of cross-linking42 is assumed or otherwise. The 19F–13C CP also indicates, based on the chemical shift range and in comparison with N115, that much of these extra signals are associated with the sidechains, and may involve some degree of branching. A detailed mapping of this unassigned fraction by spectroscopic methods would require extensive use of correlation and recoupling NMR experiments,28 and is beyond the scope of this study. If, instead, all of the signals are grouped by narrow versus broad linewidths, noting that for N115 the linewidths readily partition into a group with <1500 Hz and >1500 Hz (see Table S1†), N115 has a proportion of 0.85:
0.15, narrow-to-broad, whereas F1850 yields 0.57
:
0.43 for the same. The linewidth can therefore be employed as a measure of relative configurational disorder (i.e. an approximate metric for the density/tortuosity of the corresponding membrane), furnishing a means of linking the microstructural characteristics of the membrane to its ability to ensnare large, solvated species and thereby mitigate crossover. By this lens, F1850 intrinsically differs from N115, particularly with regard to its sidechain moieties (based on the typical chemical shift range of the unassigned signals for both 19F and 13C). This added sidechain disorder and potential branching/pseudo-crosslinking may restrict the freedom of microphase separation, congest the ion channels, and limit their permeability to polysulfide species. Correspondingly, this will also limit the volumetric extent of the hydrophilic domains in the membrane via the additional complexity of the sidechain interactions/tangling, which will reduce the mean free path of diffusing polysulfide species and attenuate their traversal of the full breadth of the membrane.
ATR-FTIR analysis was then carried out to identify the different structural features in F1850. The ATR-FTIR spectra and the assignment of major peaks of F1850 and N115 membranes are displayed in Fig. 2f. Nafion membranes have been extensively characterized by this technique and our result agrees well with the previously obtained examples,43–45 which assisted with the peak assignments. First, the major peaks of N115 also occur in F1850 with similar patterns, including vasCF3 at 1319 cm−1, vsCFbackbone at 1304 cm−1, vasCF2 at 1197 cm−1, vsCF2 at 1143 cm−1, vasC–O–C at 1134 cm−1, vsSO3− at 1056 cm−1, vsCFbranch at 982 cm−1, and vsC–O–C at 967 cm−1. Second, a few new peaks are observed in F1850, such as those highlighted by green boxes at 1403, 1073, 878, and 837 cm−1. These new features point to the major structural differences between F1850 and N115, which may be associated with the more disordered sidechain of F1850. However, it remains a challenge to precisely identify these functionalities in the polymer.
![]() | ||
Fig. 4 CV curves of 50 mM K3Fe(CN)6 and 50 mM Na2S2, respectively, in 1 M supporting salts at 20 mV s−1. |
The permeation tests have offered a preliminary indication of the membranes' polysulfide retention abilities under simplified conditions where the concentration gradient across the membrane is the major driving force for crossover. Flow cell cycling involves additional operation factors that may exaggerate polysulfide crossover, such as the presence of an electric field, active Na+/K+ transport, electrolyte circulation, different counter-electrolyte settings, etc. To demonstrate the practical feasibility of F1850, the S/Fe flow cells assembled with this membrane were evaluated under RFB-relevant conditions. Again, N115 was tested as the baseline.
The sluggish redox kinetics of polysulfide conversion requires the use of an electrocatalyst-decorated porous electrode in flow cells. The cobalt-decorated SIGRACELL® graphite felt electrodes (Co/GFD2.5) were prepared by thermally reducing the Co(NO3)2-impregnated GFD2.5 in a 4% H2/argon atmosphere. A Co loading of 11.6 mg cm−2 was obtained. The morphology of the Co coating is revealed with FE-SEM and is shown in Fig. S8.† Irregularly shaped Co clusters sized at several tens of microns are randomly distributed on carbon fibers. Micron-scale etches on the fiber surface are also observed in the FE-SEM image. These functional modifications boost the surface area of the electrode structure and increase the active sites for electrocatalysis.
The S/Fe flow cells with F1850 and N115, respectively, were assembled in the fully charged state with 1 M Na2S2 anolyte and 1 M K3Fe(CN)6 catholyte (Fig. 5a). The use of mixed Na+ and K+ cations has been reported capable of boosting the solubility of K4Fe(CN)6 from ∼0.7 M up to ∼1.5 M.49,50 In the flow cell, the through-membrane transport of Na+ and K+ ions will result in mixed cations. Therefore, cycling with 1 M catholyte will not cause precipitation of K4Fe(CN)6, which was demonstrated in a previous report under similar conditions.14 Fig. S9† shows the Nyquist plots and the fitting with the classical Thevenin equivalent circuit model. F1850 yielded a higher ohmic ASR (RΩ, 6.34 Ω cm2) than N115 (2.77 Ω cm2) to the S/Fe flow cell, and the charge transfer resistivities (RCT) were similar (∼0.4 Ω cm2) for these two membranes. Thus, the RΩ appeared to dictate the flow cell efficiency. As shown in Fig. 5b, F1850 exhibited an inferior rate capability compared to N115 under the same test conditions. For the F1850-based flow cell, as the current density increased from 10 to 80 mA cm−2, the voltage efficiency (VE) was decreased from 75.9% to 25.6% and the energy efficiency (EE) from 75.0% to 25.5%. The catholyte material utilization, based on the discharge capacity, was generally high and slightly dropped from 98.0% to 85.3% with the current density. In contrast, the N115-based flow cell showed decreases in VE from 80.6% to 49.1%, in EE from 79.1% to 48.8%, and in catholyte utilization from 96.1% to 89.7%. In addition, F1850 generally afforded higher coulombic efficiencies (CEs) than N115 at these current densities, which agrees well with the better polysulfide blocking ability of the former membrane. Well-defined single plateaus were observed in the voltage curves of both flow cells at these tested current densities (Fig. 5c and d), although multiple redox peaks were seen in the CV curve of polysulfides (Fig. 4). This is because stoichiometrically excessive polysulfide was used to compensate for its probable crossover loss during cycling. Such electrolyte settings have been commonly used for the testing of polysulfide RFBs in most of the reported studies.
Long-term cyclability is particularly important for bestowing polysulfide RFBs with commercialization relevance. The better polysulfide blocking ability of the F1850 membrane is expected to lead to more stable cycling. To test this hypothesis, continuous galvanostatic cycling at a constant current density of 20 mA cm−2 was carried out on the S/Fe flow cells. Fig. 6a displays the cycling efficiencies and capacities as functions of cycle number and test time for the F1850-based S/Fe flow cell. With an average CE of 99.2% and EE of 62.5%, the flow cell demonstrated a stable capacity profile for approximately 300 cycles over a test time of 49 days, leading to an average capacity fading rate as low as 0.029% per cycle and 0.176% per day. In contrast, the parallel test with N115 yielded a CE of 98.5% and remained stable for only 120 cycles and 19 days (Fig. 6b). Such an extended time duration with the well-maintained capacity at ∼97% SOC reflects the inherent polysulfide retention attribute of F1850. This performance is significantly better than that of most reported polysulfide RFBs and places F1850 as an exceeding candidate for permselective membranes. The cycling results confirm our claim of the favorable properties of F1850 toward advancing polysulfide-based RFBs.
![]() | ||
Fig. 6 Long cycling of S/Fe flow cells using (a) F1850 and (b) N115 at 20 mA cm−2. The anolyte was 1 M Na2S2 (15 mL) and the catholyte was 1 M K3Fe(CN)6 (10 mL). |
It is noted that the VEs of both flow cells were continuously decreased as the cycling proceeded. From the voltage curves of selected cycles, gradual increases in the cell overpotential were observed over the course of the test (Fig. S10†). The major reason for this phenomenon is polysulfide crossover. The leaked polysulfide reacted with the K3Fe(CN)6 catholyte forming an insulating elemental sulfur solid on the surfaces of the GFD2.5 electrode and membrane. The sulfur deposition occurred way before the capacity started to fade. According to the FE-SEM images and EDX sulfur mapping shown in Fig. S11,† sparsely distributed sulfur deposits were already non-negligible at the 120th cycle for the F1850-based S/Fe flow cell. The accumulation of sulfur on the electrode surface eventually resulted in fluctuations or fading of the capacity and efficiency. After a 30% loss of the capacity at the 530th cycle (i.e., the ∼80th day), sulfur particles were detected at the catholyte side – on the membrane and inside the electrode (Fig. S12†). For the N115-based flow cell, sulfur deposition was even more serious. Thick layers of sulfur were visually observed at the catholyte side after only 39 days of cycling when the capacity faded to ∼30% of the initial value (Fig. S13†). Clogging of the flow path can be easily envisaged.
An additional note is that the cycling stability is moderately affected by the polysulfide ion size, which was demonstrated in N115-based S/Fe flow cells. When polysulfides of different chain lengths were used (in large excess) under otherwise identical conditions, different cycling results were obtained. As shown in Fig. S14,† the S42−/Fe(CN)64− flow cell maintained stable capacity for a moderately longer duration of 24 days (i.e., 125 cycles). In stark contrast, rapid capacity degradation was identified for the S2−/Fe(CN)63− flow cell. Intuitively, this result is primarily because the smaller polysulfide size allows easier access into the transport channels in the membrane. Yet, other factors may simultaneously play a role in causing the unstable capacity. It's been reported that S2− is more vulnerable than longer polysulfides to hydrolysis forming HS− or even gaseous H2S, which is a known origin of the chemical instability of monosulfides in near-pH-neutral solutions.51–53 This adverse effect may explain the drastically different capacity behavior of the S2−/Fe(CN)63− flow cell, i.e., the continuous capacity fading even from the first cycle. This possibility also highlights the necessary electrolyte control to maintain the chemical stability of polysulfides.12 Nevertheless, a potential implication from the only mild improvement in cyclability is that electrolyte engineering aimed at enhancing polysulfide retention should rely more on other electrolyte factors than on polysulfides themselves.
We should point out that this F1850 membrane is limited by the well-known conductivity–selectivity trade-off in ion exchange membranes. Because slow kinetics is also a critical challenge for polysulfide redox, the reduced membrane conductivity aggravates the low efficiency issue of polysulfide RFBs. Therefore, advanced membrane approaches need to be developed to achieve both selectivity and conductivity simultaneously.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d5ta05404d |
This journal is © The Royal Society of Chemistry 2025 |