Gui Lia,
Artem Iakunkova,
Nicolas Boulangera,
Oana Andreea Lazarb,
Marius Enachescub,
Alejandro Grimmc and
Alexandr V. Talyzin*a
aDepartment of Physics, Umeå University, Umeå, SE-90187, Sweden. E-mail: alexandr.talyzin@umu.se
bCenter for Surface Science and Nanotechnology, University Politehnica of Bucharest, Splaiul Independentei 313, Bucharest 060032, Romania
cDepartment of Forest Biomaterials and Technology, Biomass Technology Centre, Swedish University of Agricultural Sciences, Umeå, SE-901 83, Sweden
First published on 12th May 2023
Activated carbons have been previously produced from a huge variety of biomaterials often reporting advantages of using certain precursors. Here we used pine cones, spruce cones, larch cones and a pine bark/wood chip mixture to produce activated carbons in order to verify the influence of the precursor on properties of the final materials. The biochars were converted into activated carbons with extremely high BET surface area up to ∼3500 m2 g−1 (among the highest reported) using identical carbonization and KOH activation procedures. The activated carbons produced from all precursors demonstrated similar specific surface area (SSA), pore size distribution and performance to electrodes in supercapacitors. Activated carbons produced from wood waste appeared to be also very similar to “activated graphene” prepared by the same KOH procedure. Hydrogen sorption of AC follows expected uptake vs. SSA trends and energy storage parameters of supercapacitor electrodes prepared from AC are very similar for all tested precursors. It can be concluded that the type of precursor (biomaterial or reduced graphene oxide) has smaller importance for producing high surface area activated carbons compared to details of carbonization and activation. Nearly all kinds of wood waste provided by the forest industry can possibly be converted into high quality AC suitable for preparation of electrode materials.
Extremely high surface area is also required in sorption applications, e.g. for improved hydrogen storage and for efficient removal of various contaminants from aqueous waste solutions. It is well known by now that physisorption of hydrogen is almost linearly proportional to SSA both at ambient and liquid nitrogen temperature.9 This trend is valid for all kinds of nanostructured carbon materials including AC, carbon nanotubes (CNT) and graphene-related materials.10,11
Many nanostructured carbon materials had emerged over past two decades, but ACs remain to be materials with highest SSA values, in some reports as high as ∼3000–3900 m2 g−1. It should be noted that at least part of reports for extremely high SSA of AC (>3500 m2 g−1) is likely related to incorrect procedures applied to analysis of sorption isotherms. Verified procedure specific for microporous materials has been designed by J. Rouquerol et al.12 The procedure allows reading reproducible values of SSA in microporous materials by selecting specific interval of relative pressures for a BET plot.12 However, there is no doubt that AC with BET SSA exceeding 3000 m2 g−1 is possible to synthesize using carefully tuned activation procedures.5,13
Activated carbons are commonly produced using large variety of bio precursors, most commonly plants, wood or different kinds of organic waste, e.g. old tires,14 food waste,15 spent tea leaves,16 lignin,17 olive pits,18,19 waste tea20 and coffee,19 corn cob,21 coconut shell,13 sawdust.22 In this case converting bio precursors into AC is not only providing useful materials but also contributes to waste recycling.
High surface area is achieved after carbonization of bio precursor by annealing under inert gas and activation of so prepared biochar using variety of methods. Most common in industrial applications is steam activation typically providing AC with surface area of about 2000 m2 g−1.23–25 The highest BET SSA is found in AC produced using KOH activation,26 and the method is based on mixing biochar with KOH in different proportions and annealing the mixture at elevated temperatures (∼700–1000 °C). Possibly the most studied high surface area AC (∼3000 m2 g−1) was commercially produced for period of time using petroleum coke under the name “Maxsorb”.27 AC with similarly high surface area was produced also starting from several bio precursors and, more recently, using reduced graphene oxide (under the name “Activated Graphene”).28–35
It should be noted that achieving BET surface area values on the level of 3000 m2 g−1 requires careful tuning of several parameters of KOH activation procedure. Most important parameters are: activation temperature, KOH to carbon ratio and homogeneous distribution of KOH in the carbon precursor.3,10,11
A very high surface area has been reported recently for AC produced using KOH activation of pinecone biochars.36–40 Pinecone derived AC have been studied previously for removal of contaminants from aqueous solutions by adsorption.22,36,38,40–43 and for the preparation of supercapacitors electrodes.39,44–48 However, most of these studies were performed with AC with SSA (1000–2000 m2 g−1)39,45,48 and single report on material with ∼2500 m2 g−1.44
Recently we demonstrated that KOH activation tuned previously for preparation of “activated graphene” could be used to produce AC with similarly high BET SSA of about ∼3000 m2 g−1 starting from pinecone biochar.49 The pinecone AC showed similar to “activated graphene” pore size distribution and good performance in supercapacitor electrodes produced by both pellet forming and spray deposition49 using aqueous dispersions.50,51
The pinecone derived AC provides an example of inexpensive alternative for sorption applications and electrode preparation for energy storage applications. However, it is not clear how much the properties of the pinecone AC are determined by the precursor biomaterial and by details of KOH procedure. Literally hundreds of studies have been published in recent years claiming advantages of this or other bio precursors for preparation of AC. However, systematic studies focused on using identical activation methods with biochars of different origins are very scarce and not available (to the best of our knowledge) for high surface area AC (∼3000 m2 g−1).
In this study, we produced AC with an exceptionally high BET surface area up to ∼3500 m2 g−1 using several common in Scandinavia by-products of wood industry, typically considered as waste. Identical KOH activation was applied to biochars produced using identical procedure starting from pinecones, spruce cones, larch cones and bark/wood chips. The ACs produced from all tested precursors appeared to be very similar except for somewhat different ash content. Therefore, it is demonstrated that details of KOH activation are a lot more important for synthesis of high surface area AC compared to the kind of bio precursor. The AC was used as electrode material in supercapacitors demonstrating excellent energy storage parameters. Hydrogen sorption was measured for the high surface area AC using volumetric and gravimetric tank methods. It is demonstrated that adding AC to the storage tank provides ∼9% increase in amount of stored hydrogen compared to compressed H2.
The two step KOH procedure used in this study for producing AC with extremely high surface area was earlier optimized using experiments with rGO precursor.3,10,11 It should be noted that precise control of KOH to carbon ratio is required for achieving highest surface area. This control is impossible in one step procedure where KOH is added directly to bioprecursor. The one step procedure results in materials with significantly smaller surface area (<2000 m2 g−1).49,52 Using ethanol as solvent for KOH is justified by rapid evaporation and more homogeneous distribution of KOH in the powder biochar sample.
The activation was performed using following temperature profile: heating up under argon flow from room temperature to 200 °C and maintaining at this temperature for 30 min in order to remove water and volatile impurities; slow heating up to 850 °C (2 hours and 40 minutes); maintaining at 850 °C for 3 hours before cooling down slowly. The activated materials were rinsed with acetic acid solution (10 vol%) in order to remove acid soluble impurities. The last step was vacuum filtration of the solution with 4 L distilled water to wash out remaining KOH and water soluble impurities.
Typically the starting amount of precursor was about 15–16 g resulting after carbonization in about 9–11 g of biochar. Four grams of biochar were typically used in each activation batch resulting in about 1–1.5 g of AC material. The overall yield of two step procedure is therefore, about 7–12%.
The four kinds of activated carbon from the cones of pine (Scots, Pinus sylvestris), spruce (Picea abies), silver spruce (Picea pungens) and larch (Larix decidua) trees were named as PAC, SAC, SSAC and LAC, respectively. Another type of tested precursor was bark with wood chips produced as waste in a wood processing company. The bark chips from Scots pine (Pinus sylvestris) provided by Ekenäs Timber AB, Ekenässjön, Sweden, were ground and sieved before carbonization. The samples with a particle size of 2–3 mm were used for producing the activated carbon, named as SBAC.
XPS spectra were recorded using a Kratos Axis Ultra electron spectrometer with a monochromatic Al Kα source at 150 W and processed by the Kratos software. Raman spectra were obtained from Raman spectrometer (Ranishaw Qontor) equipped with 532 nm solid state laser.
Cyclic voltammetry (CV), Galvanostatic charge–discharge cycling (GD), electrochemical impedance spectroscopy (EIS) and cycling stability tests were performed using Iviumstat potentiostat and symmetric cells with two identical AC-based electrodes, two stainless steel current collectors and a glass fiber membrane separator (Whatsman). 6 M KOH solution and 1 M tetraethylammonium tetrafluoroborate in acetonitrile (TEA TEF) electrolytes were used. For EIS, 10 mV sinusoidal amplitude alternating voltage at frequencies from 1 MHz to 10 MHz was applied. The stability test with 10000 charging and discharging cycles was done at 3 A g−1 in aqueous electrolyte and at 5 A g−1 in organic electrolyte.
The main purpose of gravimetric tank method (GTM) is to evaluate directly how much hydrogen is stored in the given tank with and without hydrogen storage material. The method involves a direct weight measurement of small volume tank. The weight of the tank is measured first in empty state with vacuum inside, next for the same tank filled with hydrogen under required pressure (120 bar) and finally the tank filled with material after dozing with hydrogen at 120 bar. The weight of hydrogen inside a material-filled test tank is compared then to the weight of the same tank filled with compressed hydrogen in the absence of loaded material.
The hydrogen storage capacity of material is then characterized relative to the amount of hydrogen stored in pure compressed gas using parameter introduced in our earlier study under the name hydrogen storage “Gain”. For example, the tank filled with material and Gain value of 10% will provide a 10% increase in the weight of hydrogen compared to material-free tank filled with hydrogen gas the same P–T conditions.
The test tank is sealed using standard valve and disconnected from the hydrogen dozing system for the weight measurements.54 Weight measurements were performed using a Radwag AS 510.3Y balance with a maximum load of 510 g and readability of 0.1 mg. All weight measurements were repeated 10 times and averaged to improve repeatability (see ESI† for details).
The GTM requires relatively large amount of material for reliable evaluation of H2 storage capacity of material (limited by the balance repeatability). The AC sample used in the H2 sorption test was produced in a separate synthesis batch using larch cones as a precursor. The weight of the sample after degassing was 0.9136 ± 0.0004 g. The inner volume of this assembly measured by pycnometry was 10.49 ml with a useful volume (volume occupied by material) of ∼7.2 ml, sufficient to accommodate about 1 g of AC material. Gas dosing was performed using a Hiden Isochema IMI volumetric system with a typical pressure of 120 bar and temperature of 23 °C. The sample was covered by a layer of silica wool in order to prevent outflow of material in the process of gas compression and pressure release.
Total Storage Density (TDS) specific for the given tank is defined as
TSDtank = mtotal/Vtank | (1) |
The material related Tank Storage Density (TSDmat) is defined as a density of H2 gas within the volume occupied by the sample (useful volume Vuseful). The dead volume (Vdead) is the part of the tank not filled with material due to sample size limitations (incompletely filled volume) and includes the volume of connection pipes and the volume inside of the valve. The total weight of hydrogen in the material filled volume (TSDmat) can be calculated by subtracting the amount of hydrogen stored as gas in the dead volume Vdead from the total storage mtotal (1):
TSD = (mtotal − Vdead × ρH2(P,T))/Vuseful | (2) |
Gain value shows the change in the density of hydrogen inside of the test tank due to addition of material relative to the density of hydrogen gas at given P–T conditions. It is expressed in percent. A positive value of Gain will be found for materials, which store hydrogen better compared to compressed gas and negative for materials, which store hydrogen less well. This parameter is calculated as the difference in hydrogen storage between a tank filled with a material and an empty tank with the same volume (3).
Gain = 100% × (TSD − ρH2(P,T))/ρH2(P,T) | (3) |
Fig. 1 (a) N2-sorption isotherms; (b) pore size distribution (left axis) and cumulative surface area (right axis) of PAC, SAC, SBAC and LAC samples (see Table 1 for meaning of abbreviations). |
Rather similar pore size distribution observed in all samples independently on the type of precursor points out to essentially the same nature of AC. Moreover, the pore size distribution and BET SSA of AC prepared in our experiments are rather similar to those of “activated graphene” (Fig. S3†), material prepared by KOH activation of rGO and characterized in detail in earlier studies.10,29,30
The scatter in the SSA values measured from samples listed in Table 1 is found to correlate with amount of non-combustible impurities (ash) determined using TGA traces recorded under air up to 700 °C (example shown in Fig. 2).
The C 1s XPS spectra show main peak from carbon in non-oxidized state (285.1 eV) and small components typically assigned to carbon with oxygen functional groups (Fig. 3a).
Fig. 3 XPS spectra for SAC sample: (a) C 1s and (b) O 1s. Survey XPS spectra can be found in ESI File.† The overall C/O ratio is calculated using integrals from C 1s and O 1s parts. |
Oxygen part of spectra shows one major feature assigned to carbon–oxygen bonds and small feature due to oxygen bound to aluminium impurity (Fig. 3b). Aluminium impurity is likely from the ceramic of boats and tube furnace. Survey XPS spectra do not reveal other impurities (Fig. S8†). The Al impurity was found to be in a form of Al (OH)3 as revealed by XRD. The AC itself showed no diffraction peaks in XRD patterns with some diffuse scattering in the low angle region (ESI File†).
Note that XPS is surface method and needs to be used for evaluation of bulk impurity content with some caution. Since the starting cones were collected outdoors in natural environments, these bio precursors are always contaminated with impurity particles, some non-combustible inorganic contaminants are always present also in the wood itself and finally some contamination is added in a process of carbonization/activation (e.g. materials from ceramic boats and furnace tube). Therefore, it is informative to calculate SSA value for only carbon part of samples excluding the weight of inorganic non-combustible contaminants. The SSA values in the range ∼2800–3500 cm3 g−1 were found for samples listed in the Table 1. The AC produced from bark shows smallest SSA value. The difference in SSA AC produced using cones is at least partly related to the ash content. Excluding non-combustible part from calculation provides “ideal” SSA values of cone derived AC in the range ∼3100–3700 m2 g−1. Some similarity in pore size distribution and surface area could possibly be anticipated for AC produced starting from cones of different kinds of tree species. However, nearly identical properties were found also for AC produced using bark/wood chip mixture, the precursor very different compared to cones and collected directly at the waste disposal locality.
Characterization of AC produced using precursors listed in the Table 1 was performed also using Raman spectroscopy (Fig. 4) and SEM/TEM imaging. The Raman spectra of biochars and AC's are typical for amorphous carbon materials exhibiting broad D- and G-modes. Raman spectra of SAC and spruce biochar are shown in the Fig. 4a for the region from 500 cm−1 to 3200 cm−1. The value of intensity ratio ID/IG is changed from 3.04 in the biochar to 4.75–5.78 as a result of KOH activation (Table 2), indicating higher purity of carbon in the activated materials. Once again, Raman spectra recorded from all types of AC showed similarity with relative intensity ID/IG nearly identical independently on the used precursor.
Samples | D peak (cm−1) | G peak (cm−1) | ID/IG |
---|---|---|---|
PAC | 1336.1 | 1591.5 | 5.78 |
SAC | 1334.0 | 1589.3 | 4.75 |
LAC | 1334.9 | 1590.4 | 4.83 |
SBAC | 1334.5 | 1590.2 | 4.90 |
Carbonized spruce | 1356.3 | 1592.0 | 3.04 |
Morphology of biochars produced by carbonization of cones is shown in Fig. 5. Some difference in micrometer-scale texture of samples was noted in SEM images of cone biochars. E.g. carbonized cones displayed a special hive-like features with multi-micrometer diameter voids inherited from precursor cones.
Fig. 5 SEM images of precursor biochars of (a) – PAC, (b) – SAC, and (c) – LAC samples (before KOH activation). |
However, the biochars were ground into powder before activation thus minimizing the influence of micrometer-scale architecture of these samples. High-resolution imaging provides more relevant information about the structure of materials on nm-size scale. The SEM images recorded with best resolution reveal highly porous structure with nm-size pores in all AC materials (Fig. 6) as expected from the analysis of pore size distribution provided by N2 isotherm (Fig. 1). The sample with highest SSA value (SAC) was also studied in more details using high resolution STEM imaging not achieved using SEM. High abundance of pores with size of few nm is revealed by these images (Fig. 7).
Fig. 7 High resolution Scanning Transmission Electron Microscopy (STEM) images recorded from SAC sample: (a) lower resolution, (b) higher resolution. |
Second method was to prepare aqueous dispersions using graphene oxide as a main surfactant and to blade coat the electrodes on current collector. Results of testing for supercapacitors prepared using pellet electrodes with 6 M KOH electrolyte are summarized in the Fig. 8 and Table 3.
Samples | Specific capacitance (F g−1, at 0.5 A g−1) | Corrected specific capacitance (F g−1, at 0.5 A g−1) | Energy density (W h kg−1, at 30 A g−1) | Power density (kW kg−1, at 30 A g−1) | Response time (s, in KOH) | Capacitance retention (%, after 10000 cycles at 3 A g−1) |
---|---|---|---|---|---|---|
PAC | 160.79 | 178.60 | 4.09 | 12.80 | 0.08 | 88.31 |
SAC | 174.44 | 180.17 | 4.12 | 12.31 | 0.10 | 93.10 |
LAC | 159.46 | 187.12 | 3.47 | 12.00 | 0.10 | 93.53 |
The CV curve (Fig. 8a) demonstrates the shape close to rectangular as expected for EDLC type of supercapacitors. Analysis of charge–discharge curves (Fig. 8b) was used to evaluate energy storage parameters of electrodes. Highest gravimetric capacitance was found for SAC sample (∼174 F g−1), also the sample with highest surface area and purity.
As expected, AC materials with higher amount of ash impurities showed smaller absolute values of capacitance. The AC samples showed good cycling stability with 88–93% retention after 10000 cycles (Fig. 8e) The Nyquist plots (Fig. 8f) show low-amplitude semicircle reflecting good conductivity of electrodes and low interface resistance. The cells showed similar response in the whole measured frequency range. Impedance data were plotted also as frequency-dependent complex capacitance which allows to evaluate the efficiency of a capacitive rate response55 and as an imaginary capacitance plot, which peaks at a knee frequency (or its reciprocal characteristic capacitor response time τ) separating the domination of capacitance and resistance. The SAC electrodes show the most efficient capacitive response in agreement with highest capacitance values among other samples. The capacitive response (Fig. 8h) shows rather rapid and similar characteristic time of 0.06–0.1 s for all cone-derived electrodes. Energy storage parameters of SC with pellet AC electrodes are summarized in (Fig. 8d).
Once again, the highest energy density is observed for the AC material with highest measured surface area. Using organic electrolytes allows to achieve better energy storage parameters due to broader voltage window.
Therefore, we tested four samples of cone-produced AC in SC with pellet electrodes and TEA-BF4/acetonitrile electrolyte (Fig. 9a–h). Once again, the samples with higher SSA values generally showed better performance and higher energy storage. Excellent gravimetric capacitance of ∼122 F g−1 was achieved for SAC samples. Table 4 shows also values of gravimetric capacitance calculated only for carbon part of samples (excluding ash content).
Samples | Specific capacitance (F g−1, at 1 A g−1) | Corrected specific capacitance (F g−1, at 1 A g−1) | Energy density (W h kg−1, at 30 A g−1) | Power density (kW kg−1, at 30 A g−1) | Response time (s, in KOH) | Capacitance retention (%, after 10000 cycles at 5 A g−1) |
---|---|---|---|---|---|---|
PAC | 118.01 | 131.08 | 20.93 | 22.43 | 0.08 | 91.28 |
SAC | 122.53 | 126.55 | 19.82 | 16.13 | 0.50 | 78.19 |
LAC | 105.61 | 123.93 | 16.71 | 21.11 | 0.10 | 90.10 |
The similarity between performances of electrodes prepared using three different bio precursors become even more obvious in this representation. Notably, the samples with the highest surface area (SAC) showed somewhat worse cycling stability and stronger drop of gravimetric capacitance at higher current density as compared to PAC sample. Ragone plot (Fig. 9d) shows that high energy and power density are achieved using cone-derived AC electrodes, with best sample showing energy density 20.9 W h kg−1 and power density 22.4 kW kg−1.
Energy storage parameters of supercapacitors prepared using blade coated water based dispersions based on AC were also analysed. Very similar electrochemical performance was observed except for somewhat smaller energy storage values as expected due to influence of added components. Fumed silica and graphene oxide are required to make the aqueous dispersions stable but these additions are not contributing to energy storage. The electrodes prepared using blade coating showed also less good cycling stability, most likely due to GO addition. Therefore, this study was mostly focused on pellet electrodes (see ESI File† for the data with performance of blade deposited electrodes).
It is interesting to note that the energy storage parameters recorded here for AC produced from essentially cost free wood waste precursors are very similar to the performance of “activated graphene” in supercapacitor devices.3,28,29 “Activated graphene” is produced starting from graphite oxide which is explosively exfoliated to produce rGO powder. KOH activation of rGO results is a material with BET surface area, pore size distribution and performance in supercapacitors which are very similar to those reported here for ACs. Using essentially cost free wood waste instead of expensive and still rather exotic graphite oxide as a precursor for high surface area carbons is clearly an advantage for practical applications in energy storage.
Hydrogen sorption tests using volumetric method were performed for the LAC sample with SSA (2981 m2 g−1) at ambient and liquid nitrogen temperature. The same samples were tested for hydrogen storage also using Gravimetric Tank Method (GTM). Note that this sample is from different batch compared to the one described in Table 1 since relatively large samples are needed for good precision measurements in the GTM method. The excess H2 sorption of 0.96% (120 bar) was measured for this sample by volumetric method at ambient temperature (296 K) and 4.6 wt% at liquid nitrogen temperature. Volumetric test was also performed at ambient temperature for the sample with SSA of 2656 m2 g−1 (activated carbon from full larch cone) and showed 0.87 wt% (120 bar) excess uptake.
The volumetric method allows to measure H2 uptake as a surface excess value. Direct evaluation of hydrogen storage as amount of hydrogen storage inside of tank is enabled by Gravimetric Tank Method (GTM) developed in our earlier study.
The main goal of this method is to measure hydrogen storage Gain defined as an ability of material to increase (or decrease) amount of hydrogen stored inside of tank compared to hydrogen stored in the same volume as a pure H2 gas. The Gain value of 8.9% was evaluated using GTM method for LAC sample with SSA of 2981 m2 g−1. According to definition of Gain, the tank completely filled with powder of LAC at 120 bar H2 pressure will store 8.9% more hydrogen (by weight) compared to the same tank filled with pure gas (no material added). The GTM also allows to calculate excess wt% providing verification of isotherm recorded using volumetric method (red points in the Fig. 10a). It should be noted that high surface area AC is one of relatively few materials, which demonstrate positive values of Gain at ambient temperature. That is, adding material to the tank help to store more hydrogen due to large surface excess uptake. Many materials capable to physisorb hydrogen (e.g. reduced graphene oxide rGO) were demonstrated to exhibit negative gain, when adding material to the tank results in smaller total hydrogen storage.
The Gain takes into account volumetric properties of material and explicitly includes bulk density into calculation of volume occupied by material. Therefore, Gain value can be improved by compaction of powder into pellets. Low bulk density of high surface area AC (∼0.2 g ml−1 for sample used in this study) is typical obstacle in using this material for practical gas storage application. It is very likely that bulk density of our AC can be increased by 2–3 times using pressurizing into pellets or by using more advanced monolith compaction methods reported in literature. In this case, the Gain value is expected to increase up to 15–20% thus providing sizable increase of hydrogen storage capacity for material-filled H2 tanks.
It is interesting to compare hydrogen uptake measured in our experiments for cone-derived AC with overall trends reported earlier for other nanostructured carbon materials. The Fig. 11 shows data points obtained earlier for hydrogen uptake at 120 bar and 296 K for several types of materials: rGO at lower SSA values, SWNT's, commercial activated carbon with SSA ∼2000 m2 g−1 and variety of “activated graphene” samples produced at different conditions. It should be noted that “activated graphene” samples with highest SSA values shown in this graph were obtained using KOH procedure (applied to rGO precursor) very similar the procedure used in this study for activation of biochars. Fig. 11 shows that hydrogen sorption properties of LAC are following general trends in H2 uptake vs. SSA values and provides values which are among the best reported. However, activated carbons are materials, which are a lot less expensive in production compared to e.g. CNT's or “activated graphene.” Moreover, activated carbons were produced in this study from bio precursors rather abundant in forest industry and currently considered as a waste.
Fig. 11 Excess hydrogen uptake measured at 120 bar H2 pressure and ambient temperature by various carbon materials vs. BET surface area. The measurements performed in this study are shown by red stars. Black symbols reflect general trend reported in our earlier study including rGO (circles) and various arGO samples produced with variations of KOH activation procedure.10 |
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3ra00820g |
This journal is © The Royal Society of Chemistry 2023 |