Rhys
Montecillo
ab,
Cheng-Sao
Chen
*c,
Kuei-Chih
Feng
ade,
R. R.
Chien
ad,
Pin-Yi
Chen
*ade and
Chi-Shun
Tu
*af
aInternational PhD Program in Innovative Technology of Biomedical Engineering and Medical Devices, Ming Chi University of Technology, New Taipei City 24301, Taiwan. E-mail: pinyi@mail.mcut.edu.tw
bDepartment of Physics, Silliman University, Dumaguete City, 6200, Philippines
cDepartment of Mechanical Engineering, Hwa Hsia University of Technology, New Taipei City 23567, Taiwan. E-mail: rickchen@cc.hwh.edu.tw
dDepartment of Mechanical Engineering, Ming Chi University of Technology, New Taipei City 24301, Taiwan
eResearch Center for Intelligent Medical Devices, Ming Chi University of Technology, New Taipei City 24301, Taiwan
fDepartment of Physics, Fu Jen Catholic University, New Taipei City 24205, Taiwan. E-mail: 039611@mail.fju.edu.tw
First published on 22nd April 2024
High energy-storage capability and electric breakdown strength are critical elements in next-generation pulse-power dielectric capacitors. In this report, perovskite (Bi0.7Ba0.3)1−xNax(Fe0.7Ti0.3)1−xTaxO3 relaxor ferroelectric ceramics (x = 0–0.3) were tailored in terms of configuration entropy from a medium entropy of 1.21R to a high entropy of 2.07R to improve energy storage. The integration of paraelectric NaTaO3 into BiFeO3–BaTiO3 results in breaking of the long-range order and formation of multiple lattice distortions toward relaxor ferroelectric characteristics. Excellent recoverable energy densities of 9.6 J cm−3 and 10.3 J cm−3 with efficiencies of 77% and 68% at 350 kV cm−1 and 550 kV cm−1 (at 10 Hz) were achieved for x = 0.15 and 0.20, respectively. Wide operating frequency (1–100 Hz) and temperature (25 °C–150 °C) stabilities were confirmed at 300 kV cm−1. Grain boundaries and nanoclusters play critical roles as electric barriers to suppress charge mobility and increase electric breakdown strength. This study presents a promising scheme to utilize high-configuration entropy BiFeO3–BaTiO3-based ceramics for high energy-density electrostatic capacitors.
The energy-storing capability (ESC) of a dielectric capacitor is characterized by high power density, but relatively low recoverable energy density (Wrec). Thus, the development of dielectric capacitors has prompted great interest in enhancing Wrec. Wrec is defined as where E, P, Pmax, and Pr are the E field, polarization, maximum polarization at maximum E field, and residual polarization after discharging (at zero E field), respectively.2 Thus, a larger Pmax and a smaller Pr (or larger ΔP = Pmax − Pr) can result in a higher Wrec. A large Pmax can be achieved by applying a higher E field, which has prompted researchers to explore strategies for enhancing the electric breakdown strength (Eb). Eb is the E-field threshold that the material can accommodate before it breaks. Another essential indicator in energy storage is the efficiency determined by the relation, η = Wrec/(Wrec + Wloss), where Wtot and Wloss are total energy density stored during charging and the dissipated energy remaining in the material after discharging, respectively.2
Bismuth ferrite (BiFeO3, abbreviated as BF) is a perovskite ferroelectric (FE) system with a large spontaneous polarization (Ps ∼ 90–100 μC cm−2) as proposed by first-principles calculation,3 and Ps ∼ 90–150 μC cm−2 in the BF thin film.4,5 Similar to Pb2+, Bi3+ in BF with the electron configuration of a 6s2 lone pair makes BF a potential substitute for lead-based materials in energy-storage applications.6 However, a pure BF ceramic suffers from current leakage and low breakdown strength due to Bi volatilization and Fe3+ to Fe2+ valence hybridization.7,8 The drawbacks of utilizing BF for energy storage have been overcome by rare-earth doping and/or designing binary materials such as BiFeO3–BaTiO3 (BF–BT).7,9–15 The BF–BT binary system exhibits a temperature-dependent morphotropic phase boundary (MPB) that permits composition-driven structural transition resulting in improved energy storage properties.1 The substantial polarization of Pmax > 40 μC cm−2, moderate remanent polarizations (Pr), and moderate Eb were reported at the MPB in Nd-doped BF–BT ceramics.16 The addition of 3 mol% Nb2O5 in 0.65BF–0.35BT reduced Pr from ∼40 μC cm−2 to 5.33 μC cm−2.17 These results are due to disruption of the long-range order and simultaneous induction of a relaxor ferroelectric (RFE) phase characterized by nanoscale domains and polar nanoregions (PNRs).18 Several studies have shown that (1 − x)BF − xBT undergoes a composition-driven structural transition from rhombohedral R3c at x = 0–0.25 to pseudo-cubic at x = 0.33–0.94.19–21 Although BF-BT- and BT-based ceramics have made some progress in energy storage, they still have relatively lower energy storage than lead-based systems.22–27
Since the local compositional inhomogeneity is linked to the disorder and random E field in the REFs, the design of REFs has been approached by modulating the configuration entropy (Sconfig) to improve energy-storage performance.18–21Sconfig in oxide systems is defined as where xi and xj describe the mole percentage of various elements in cation and anion sites.28R is the universal gas constant. Sconfig is divided into three levels: low entropy (Sconfig < 1R), medium entropy (1R ≤ Sconfig < 1.5R), and high entropy (Sconfig ≥ 1.5R).28 A high-Sconfig scheme has been employed in many studies because it can stabilize a single-phase structure29 and multiphase mixtures30 and enhance Eb.28
The energy-storage performance of dielectric materials mostly depends on the breakdown strength, which is influenced by dielectric characteristic, electronic bandgap, and grain size. Paraelectric NaTaO3 (NT) appears to be a good candidate for adding into BF-BT to increase Sconfig due to its incipient FE nature with a glass-like or relaxor-like phase.2,31 A high Wrec = 6.3 J cm−3 and a high η ≥ 86% were achieved in NT-doped BF-based ceramics with an Eb ∼ 425 kV cm−1 and ΔP = Pmax − Pr ∼ 38.5 μC cm−2.2 An enhanced bandgap was attributed to a high Wrec = 6.5 J cm−3 and a high η = 94% accompanied by temperature stability in a Ta-doped sodium niobate-based relaxor ceramic, (Na0.8Bi0.1)(Nb0.9−xTaxTi0.1)O3, with a high Eb = 450 kV cm−1.32 NT-doped (Bi0.5Na0.5)TiO3 exhibited a large Pmax with low Pr owing to the coexistence of R3c nanodomains and P4bm PNRs.33 These results demonstrated that NT promotes RFE behavior, causing lower Pr and higher η.
A schematic diagram of the design in this study is shown in Fig. 1 for improving energy-storage properties by tailoring Sconfig. In this study, paraelectric NT was added into 0.7BiFeO3–0.3BaTiO3 to produce (Bi0.7Ba0.3)1−xNax(Fe0.7Ti0.3)1−xTaxO3 ceramics with x = 0–0.30, yielding high Sconfig values from 1.21R to 2.07R and Eb from 175 kV cm−1 to 550 kV cm−1.
A scanning electron microscope (SEM, JEOLJSM-7610F plus) was used to image cross-section grain morphologies. A Bruker X-ray diffractometer (D8 ADVANCE ECO) was employed to investigate lattice structures assisted by Rietveld refinement analyses using HighScore Plus software version 3.0.5. A HORIBA LabRAM HR Evolution system was used with a 532 nm laser to acquire Raman scattering spectra. Using a Workstation 2000 (Radiant Technologies Precision LC II), polarization vs. E field (P–E) hysteresis loops were acquired at frequencies of 1 Hz and 10 Hz. For the P–E loop experiment, to construct the structure of a parallel-plate capacitor, samples were carefully polished to 0.05 mm thickness, which is much smaller than the electrode's diameter of 1 mm, as shown in Fig. S2 (ESI).† Au electrodes were deposited using the dc sputtering method. The nanoscale morphology and electron diffraction were investigated with a JEOL transmission electron microscope (JEM-2100 LaB6) with 200 kV-acceleration capability. A phase impedance analyzer (Wayne-Kerr PMA6420A) was used to assess dielectric properties. The electrical potential distribution on the gain matrix was probed by employing a Kelvin probe force microscope (KPFM, Bruker Multimode 8).
![]() | ||
Fig. 2 Thermally etched fractured cross-section SEM micrographs of BBNFTT-x (x = 0–0.30). Insets are the distributions of grain sizes. |
To identify the structural evolution, XRD spectra from 20 to 80 at 2θ degrees are presented in Fig. 3a, showing perovskite structures. The elaborations of (110), (111), and (200) diffraction peaks in Fig. 3b show a shift toward higher 2θ degrees from x = 0 to x = 0.30. This phenomenon indicates an A-site replacement of Bi3+ (∼1.4 Å) and Ba3+ (∼1.6 Å) in the lattice with smaller cations Na+ (<1.4 Å) causing a reduction in lattice d spacing. It is noticeable that the peaks are narrow (without splitting) suggesting an overall pseudo-cubic structure.36 In particular, extra peaks appeared at ∼29°, ∼34.5°, ∼49°, and ∼58° in NT-contained specimens, as observed in Fig. 3a and were identified as a cubic Fdm Na2Ta2O6-like structure by the refinements in Fig. S3 (ESI).† The refinement analyses suggest multiple phases and lattice distortions due to the integration of aliovalent cations, namely, Na+ into the A-site and Ta5+ into the B-site. The lattice distortions, due to the increasing configuration entropy, restrain the formation of secondary phases37 while stabilizing co-existing multiphases.28,34 Rietveld refinement analyses were carried out to quantify the percentages of various symmetries, as displayed in Fig. 3c. The Rietveld refinement reveals the coexistence of nonpolar cubic Pm
m, FE rhombohedral R3c, and nonpolar orthorhombic Pnma symmetries at x = 0. The variation in phase percentages, as presented in Fig. S4 (ESI),† shows a rapid decline in R3c symmetry accompanied by an increase in orthorhombic Pnma symmetry with integration of an increasing amount of NT, while cubic Fd
m symmetry exhibits rapid growth from x = 0.20 to x = 0.30.
To investigate the local structural alteration, Raman spectra were acquired, as displayed in Fig. 4a, showing broad profiles with wavenumbers ranging from 50 cm−1 to 900 cm−1. The 13 Raman-active modes of pure BF in the rhombohedral R3c space group are indicated in Fig. 4a. The solid curves in Fig. 4a show the fitted peaks of the Raman spectra using Peakfit software. The frequency trends of Raman modes with respect to NT content are plotted in Fig. 4b to clearly visualize gradual shifts. The broad Raman-active modes indicate that the aliovalent integration of the cations (Bi3+, Ba2+, and Na+) in the A site and (Fe3+, Ti4+, and Ta5+) in the B site has increased the configuration entropy to favor cubic symmetry.38 Typically, three zones that correspond to various cationic interactions may be distinguished in the Raman spectra of a perovskite ABO3 structure. Wavenumbers < ∼150 cm−1 were assigned to the A-site cationic vibrations (E1 and A1-1) and slower lattice motion. Wavenumbers between ∼150 cm−1 and 300 cm−1 assigned to A1-2, A1-3, and E2 are related to B-site cationic vibrations (B–O bonds), and wavenumbers > ∼300 cm−1 assigned to E3, E4, A1-4, E5, E6, E7, E8, and E9 are related to the bending/stretching modes of the BO6 octahedron.39,40 E1 and A1-1 modes show an upward trend, indicating that the lighter Na+ cations have successfully replaced Bi3+ cations in the A site of the unit cell. The downward trend in the A1-2, A1-3, and E2 modes signifies the B-site occupancy of the Ta5+ cations, which are heavier than Fe3+ and Ti4+ cations. An upward trend was obtained in the E3, E4, A1-4, E5, E6, E7, E8, and E9 modes, implying an increase in interatomic bonding strength.
The polarization vs. E field (P–E) hysteresis loops were assessed using various applied E fields at 1 Hz, as shown in Fig. 5a, to ascertain the influence of integration of NT on the ESC. As NT (x) increases, the P–E hysteresis loop evolves from a typical FE loop characterized by high Pmax and Pr at x = 0 to one with an RFE nature at x = 0.05, 0.15, then towards an almost linear loop at x = 0.20 and 0.30. Noticeably, Eb significantly rises from 175 kV cm−1 at x = 0 to 550 kV cm−1 at x = 0.20 and drops to 425 kV cm−1 at x = 0.30. The grain size is a crucial factor affecting the dielectric breakdown strength.41 The evolution of Eb is consistent with the trend in average grain sizes in Fig. 2. The larger Eb at x = 0.3, which has a larger average grain size (compared to x = 0.1 and 0.15), may result from a collective effect. Eb can be affected by intrinsic or extrinsic factors. Eb is intrinsically affected by local defect structures (such as dipolar defects and vacancies),42 secondary phase, band gap, space charge, and nanostructures (such as nanodomains and nanoclusters). Eb can be also changed by extrinsic factors, including grain size, sample thickness, and porosity. Moreover, declining R3c symmetry occurring with rising nonpolar Pnma and Fdm symmetries, as shown in Fig. S4 (ESI),† is accountable for the reduced Pmax and Pr as NT increases. The trends in Wrec (J cm−3) and η (%) vs. E field in Fig. 5b show that Wrec increases from 2.5 J cm−3 to 8.7 J cm−3, and η increases from 25% to 66% at x = 0 to x = 0.20 at Eb (or maximum E field). However, Wrec and η at Eb drop to 4.6 J cm−3 and 57%, respectively, at x = 0.30, which is attributed to the development of a nonpolar characteristic by a high configuration entropy.30 Additionally, ΔP (= Pmax − Pr) and Ec were plotted as functions of E field from 100 kV cm−1 up to Eb, as shown in Fig. 5c, which shows an increasing trend in ΔP for x ≤ 0.15 and later a decline for x ≥ 0.20 due to the decreasing Pmax. This confirms a diminishing FE R3c structure accompanied by increasing nonpolar Pnma and Fd
m symmetries as NT increases, as evidenced in Fig. S4 (ESI).† In addition, P–E loops of BBNFTT-x (x = 0.15 and 0.20) were measured at 10 Hz under various E fields, as shown in Fig. 5d, with curves of Wrec, η, ΔP, and Ecvs. applied E field in Fig. 5e and f, respectively. Recoverable energy densities of 9.6 J cm−3 and 10.3 J cm−3 with efficiencies of 77% and 68% under E fields of 350 kV cm−1 and 550 kV cm−1 were achieved at x = 0.15 and 0.20, respectively.
The evolution trends in Pmax and Pr at Eb based on Fig. 5a are plotted in Fig. S5 (ESI).† The Pmax ∼ 110.6 μC cm−2 and Pr ∼ 71.6 μC cm−2 are initially high at x = 0, as a result of the FE character. The significant reductions in Pmax (∼32.6 μC cm−2) and Pr (∼7.2 μC cm−2) at x = 0.30 indicate that the integration of NT into BF–BT is favorable for disrupting the long-range FE order toward an RFE characteristic.43 The enhanced relaxor characteristic with increasing NT content confirms that a high Sconfig, can contribute to disruption of the long-range FE order, facilitating PNR formation.29 The collapse of long-range polar order is responsible for the reduction in Pmax and Pr. Fig. S5b (ESI)† summarizes Eb, Wrec (at Eb), and η (at Eb) with increasing NT content. The energy storage parameters of BF- and BF–BT-based ceramics reported by other researches together with this work are presented in Table S1 (ESI),†14,15,44–51 which demonstrate that the results from this work are superior to reported parameters of Eb and Wrec. The energy-storage parameters are also higher than those in BNT-BT-based ceramics.52 To assess the ESC per unit E field, recoverable energy-storage intensities (ρ) were calculated based on Fig. 5b, c, e and f using the relation ρ = Wrec/ΔE (ref. 53) and are presented in Table S2 (ESI).† ΔE is Eb − Es, and Es is the initial E field. The calculated maximum ρ of 27.4 × 10−3 J kV−1·cm−2 at 10 Hz was attained at x = 0.15 and is larger than those reported in BT-, BF-, NN-, KNN- or BNT-based ceramics.53 Furthermore, a high Wrec of 10.5 J cm−3 and η of 87% were achieved in 0.62BiFeO3–0.3BaTiO3–0.08Nd(Zn0.5Zr0.5)O3 multilayer ceramic (7-active layer) at an E field of 700 kV cm−1.45 An ultrahigh Wrec of 18.5 J cm−3 with an Eb of 995 kV cm−1 was also reported in the relaxor–antiferroelectric 0.90NaNbO3–0.10BiFeO3 ceramic.54
Fig. 6a–d exhibit frequency- and temperature-dependent P–E loops for BBNFTT-0.2, showing good operating frequency (1–100 Hz) and thermal (25 °C–150 °C at 1 Hz and 10 Hz) stabilities under 300 kV cm−1. The results indicate excellent operating frequency and thermal stabilities. To confirm the reliability of experimental Eb fields, Weibull distributions have been plotted in Fig. 6e for BBNFTT-x (x = 0.1, 0.15, 0.20, and 0.30) based on relations, Xi = ln(Ei) and where i, Ei, and n are the ordinal of the sample, the Eb of the ith sample, and the total number of samples, respectively.2 The fitted Weibull moduli (β) are close to 20, suggesting the high reliability of experimental Eb values. In addition, cycling P–E loops (3 × 103 cycles) under 300 kV cm−1 at 10 Hz for BBNFTT-0.2 in Fig. 6f demonstrate that Pmax and Pr vary insignificantly with continuously increasing cycle number. The curves of Wrec and η vs. cycle number shown in Fig. 6g reveal reliable fatigue resistance with insignificant changes in Wrec or η.
To further verify the RFE characteristic and thermal stability induced by the integration of NT, the temperature-dependent dielectric permittivity (εr) and loss (tanδ) were collected at varying frequencies of 1 kHz–1 MHz, as presented in Fig. 7a–f. The dielectric permittivity maxima become broad with declining εr, as NT increases from x = 0 to x = 0.30, which is characterized by a dispersed frequency shifted towards a lower temperature. The dielectric behavior suggests that the integration of an increasing amount of NT induces RFE characteristic evolving towards a linear-like dielectric.2,43 The broad dielectric maximum can be associated with the nanoscale relaxor nature, originating from compositional disorder due to high configuration entropy. The increase in the loss (tan
δ) at higher temperatures (>300 °C) possibly originates from the ionic conductivity.55–57 The flattening and compression of the εr curve with rising temperature imply that energy-storage properties are thermally stable. Enhancement in Eb is made possible by the decline in εr as NT increases, because of the transformation of FE domains into nanodomains and/or PNRs and nanoclusters attributed to the enhanced relaxor characteristic induced by integration of an increasing amount of NT.38 The temperature coefficients of capacitance (TCC) were used to evaluate the thermal stabilities at 25 °C–125 °C, as shown in Fig. S6 (ESI).† TCC was determined using the equation
where C25 is the capacitance at 25 °C, and CT is the capacitance at temperature T. Ceramics with x = 0.20 and x = 0.30 remain thermally stable with TCC < 25%, reflecting an improved temperature stability via integration of NT.28,56
Nanoscale TEM images and selected area electron diffraction (SAED) patterns were examined in BBNFTT-0.1 and BBNFTT-0.2, as shown in Fig. 8. The typical granular microstructures of nanoclusters and dark/light-contrast structures in the grain matrix can be identified in Fig. 8a. A contrast in bright-field (BF) TEM images usually results from various intensities of diffraction. A non-uniform distribution of elemental composition inside a grain can yield dark/light contrasts in the BF TEM images. A magnified image of the nanoclusters from light-contrast region “b” in the pink box in Fig. 8a is displayed in Fig. 8b. The spot pattern of the cubic Pmm space group shown in the inset of Fig. 8b from the [100]c-zoned SAED suggests that the nanoclusters are segregated in the grain with the same crystalline orientation. The morphologies of the nanoclusters are different from nanograins, which are typically identified by ring patterns.58 The high-resolution TEM image of the nanoclusters with stripes in Fig. 8c from the dark-contrast region “c” in Fig. 8a shows a parallel-striped superlattice perpendicular to the [001] direction. The inset of Fig. 8c shows a corresponding [100]c-zoned SAED pattern with superlattice spots, suggesting a multi-symmetry structure of cubic Pm
m and orthorhombic Pnma space groups. The SAEDs along the [110]c-zone axis in Fig. 8d and e indicate the coexistence of R3c, Pm
m, and Pnma symmetries accompanied by
and
superlattice diffractions in regions “d” (light-contrast) and “e” (dark-contrast) in Fig. 8a. The antiphase oxygen octahedral tilt along the [111]c direction, designated as the a−a−a− tilt system is the origin of the
superlattice reflection marked by green circles, indicating FE R3c symmetry.55,57 Here, subscript “c” refers to the pseudo-cubic notation, while “o” stands for an odd Miller index. The octahedral tilt twin disorder in the [001]c direction gives rise to the
superlattice diffraction pattern, which reflects Pnma symmetry with the a−a−c+ mixed tilt system.59–61 According to first-principles calculation, the antiferroelectric Pnma structure is accompanied by the presence of two sublattices with opposite polarizations.62 The
superlattice ordering of cation polarization arrangements in the orthorhombic Pnma structure is connected to the stripes seen in Fig. 8d. Fig. S5a–c (ESI)† show EDX spectra, revealing that the light-contrast and dark-contrast regions from “d” and “e” areas in Fig. 8a are Ba/Ti-rich and Fe-rich, respectively. High-configuration entropy ceramics typically consist of disordered multi-cation sublattices originating from fluctuations in composition and charge.63
![]() | ||
Fig. 8 (a–c) TEM images and (d and e) SAED patterns of BBNFTT-0.1. The insets in (b) and (c) are corresponding SAED patterns. (f and g) TEM images and (h and i) SAED patterns of BBNFTT-0.2. |
The TEM images and SAED patterns for BBNFTT-0.2 were also probed, as displayed in Fig. 8f–i, revealing smaller grains, dark-contrast precipitates along the grain boundary (GB), and nanoclusters. The overall sizes of the nanoclusters are smaller than those in BBNFTT-0.1, indicating more interfaces of nanoclusters, as shown in Fig. 8g. Cubic Pmm, Fd
m, and orthorhombic Pnma space groups are revealed from [110]c and [111]c-zone-viewed SAEDs in Fig. 8h and i from “h” and “i” regions in Fig. 8f. Fig. 8h presents the SAED pattern from the smaller grain with dark contrast. The streak-like diffraction pattern indicated by yellow arrows in Fig. 8h could result from incommensurate structures along the [001]c direction in the Pnma symmetry, possibly due to disordered multi-cation sublattices resulting from fluctuations in composition and charge.63Fig. 8i shows the [111]c-zone-viewed SAED from the GB region in Fig. 8f, indicating a cubic Fd
m phase. The cubic structure precipitated along GBs may increase electric insulation to enhance Eb. The EDX spectra of BBNFTT-0.2 are displayed in Fig. S7d–f (ESI),† reflecting a stoichiometric distribution in the grain interior, an Fe/Ta-rich GB, and an Fe/Bi-rich triple junction.
To assess the electric characteristics of GBs, potential mapping images from BBNFTT-0.0 and BBNFTT-0.2 were acquired across the GBs in a 15 μm × 15 μm area, as shown in Fig. 9a and b. The electrical potential gradient (or change per unit distance) is steeper at the GBs than in the grain interior. According to the relation E = −dV/dx (potential difference per unit distance), this result suggests that the GB has higher local E fields, which can act as barriers to suppress the charge transport and thus enhance Eb. A recent study shows that interfacial polarization (between GBs and grains) plays a significant role in local E-field redistribution and affects energy-storage performance.64 Interfacial polarization (also known as Maxwell–Wagner polarization) results from space charges in inhomogeneous systems and causes non-uniform E fields along the GBs. The potential peaks suggest accumulation of positive charges at the GBs and are consistent with EDX mapping, as shown Fig. S7e and f (ESI).†
From the perspective of nanoscale structures, randomly distributed compositions with increasing Sconfig can give rise to the formation of a local random field, nanodomains, and nanoclusters.65,66 As demonstrated in Fig. 9c and d, the sizes of the nanoclusters in BBNFTT-0.2 are considerably smaller than those in BBNFTT-0. Fig. 9e shows the high-resolution TEM image of the nanoclusters in BBNFTT-0.2. The small nanocluster size in BBNFTT-0.2 indicates a higher-density nanocluster boundary that will cause local lattice distortion and elastic strain. Nanoscale segregations (such as nanoclusters) can result in a scattering effect of charge carriers to decrease the probability of electric breakdown and improve Eb.67 Thus, the superb Eb and relaxor nature can be linked to the reduced sizes of grains and nanoclusters as the amount of NT integrated increases.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4ta00921e |
This journal is © The Royal Society of Chemistry 2024 |