Open Access Article
Xiuru
Yang
a,
Mansour
Alhabradi
bc,
Anurag
Roy
b,
Manal
Alruwaili
bd,
David
Benson
b,
Hong
Chang
a,
Xiaohong
Li
b,
Asif Ali
Tahir
*b and
Yanqiu
Zhu
*a
aFaculty of Environment, Science and Economy, University of Exeter, Exeter EX4 4QF, UK
bFaculty of Environment, Science and Economy, University of Exeter, Penryn TR10 9FE, UK. E-mail: Y.Zhu@exeter.ac.uk
cDepartment of Physics, Faculty of Science, Majmaah University, Majmaah, 11952, Saudi Arabia
dDepartment of Physics, Faculty of Science, Jouf University, 2014, Sakaka 42421, Saudi Arabia
First published on 23rd December 2025
A CuBi2O4 photocathode with interconnected nanoparticle textured morphology has achieved a photocurrent density of −0.94 mA cm−2 at 0.52 V vs. RHE. It was successfully fabricated via electrodeposition using ethylene glycol (EG) containing a specific concentration of Bi(NO3)3·5H2O and CuCl2·2H2O as the electrolyte, followed by 2 h of calcination at 550 °C. Using urea as a complexing agent in the EG electrolyte enhanced the photocurrent density of the CuBi2O4 photocathode. Adding 0.15 g of urea to the electrodeposition solution improved film uniformity, enhanced PEC water splitting efficiency, and achieved a photocurrent density of −1.44 mA cm−2 at 0.52 V vs. RHE. This value is higher than those of previously reported CuBi2O4 photocathodes, which typically exhibit photocurrent densities below −1.0 mA cm−2. To understand the factors contributing to this enhanced PEC performance, this study investigated the effects of varying urea concentrations (0, 0.1, 0.15, and 0.2 g per 100 mL EG) on the crystallite domain size, morphology, surface roughness, light absorption, band gap, electronic band structure, and PEC performance. A mechanism was proposed to account for the long-term stability based on its inadequate valence band potential and irreversible degradation behaviour. This work provides insights for optimizing CuBi2O4 thin films to enhance their stability and efficiency in PEC water splitting applications.
Various techniques, including spin coating,14,15 pulsed laser deposition (PLD),16,17 physical vapor deposition (PVD) techniques (sputtering and electron beam evaporation),18–20 doctor blade coating,21,22 spray pyrolysis,23,24 chemical bath deposition,25,26 and electrodeposition,27,28 have been used to fabricate CuBi2O4 thin films. Among these, electrodeposition is particularly attractive due to its low cost, mild operating temperature, and ability to produce adherent interfaces, making it highly suitable for the creation of durable thin films. However, non-uniform electrodeposition can reduce film quality and photocathode performance. To address these issues, various strategies have been explored, including adjusting the metal-ion concentration, implementing pH control through buffering, and adding organic additives during the deposition process. Nakabayashi et al. improved the uniformity of CuBi2O4 thin films by using tartrate ions to stabilize the Cu2+ and Bi3+ in solution. This stabilization enabled controlled co-deposition of CuO and Bi2O3 at a high anodic potential.29 After annealing at 500 °C, the resulting CuBi2O4 thin films exhibited enhanced coverage, stronger adhesion, and better uniformity. Similarly, Hahn et al. prepared Cu–Bi electrodeposition baths by dissolving Bi3+ and Cu2+ in 10% nitric acid. This step was essential for fully dissolving bismuth nitrate and producing a homogeneous electrolyte, which is critical for consistent electrodeposition.27 The homogeneous bath enabled effective cathodic co-deposition and subsequent annealing yield well-crystallized CuBi2O4 thin films. Citric acid has also been employed as a complexing agent in CuBi2O4 synthesis. By forming stable metal–ligand complexes, it regulates the release of metal ions during electrodeposition, thereby promoting uniform film growth.30
Urea is a small, inexpensive, and environmentally friendly molecule with strong hydrogen-bonding and coordination abilities.31–33 Its –NH2 group can donate hydrogen bonds while its carbonyl oxygen can accept them, enabling interactions with solvent molecules and metal ions.31–33 These interactions tune the solvation environment, improve ion dissolution, and influence key solution properties such as ionic conductivity and viscosity, ultimately promoting uniform metal deposition and controlled thin-film growth.31,34,35 Motivated by these advantages, this work investigates the effect of urea at varying concentrations (0, 0.1, 0.15, and 0.2 g per 100 mL of ethylene glycol, EG) on the synthesis of CuBi2O4 thin films on FTO glass. The resulting thin films are systematically characterized in terms of crystallite size, morphology features, surface roughness, electronic band structures, and PEC performance. Based on the observed enhancements, a potential mechanism involving charge separation and transfer is proposed to account for the enhancements in PEC performance.
902, 300 W xenon lamp with an air mass (AM) of 1.5 was used to approximate a light intensity of 100 mW cm−2 (1 sun condition). The photocurrent density vs. EAg/AgCl curves (J–V plots) were obtained by using the linear sweeping voltammetry (LSV) mode at a constant scan rate of 10 mV s−1 from −0.50 to 0.2 V vs. Ag/AgCl. At a bias potential of −0.45 V vs. Ag/AgCl, the electrochemical impedance spectroscopy (EIS) curves were recorded with a frequency range of 100 kHz to 0.1 Hz under a modulation amplitude of 10 mV. To compare with literature results, the experimental potential measured versus the Ag/AgCl reference electrode (EAg/AgCl) was converted to the potential versus the reversible hydrogen electrode (ERHE) via the following Nernst eqn (1):36,37| ERHE = EAg/AgCl + E0Ag/AgCl + 0.0591 V × pH (E0Ag/AgCl = 0.1976 V vs. NHE at 25 °C) | (1) |
![]() | ||
| Fig. 1 (a) XRD patterns and (b) corresponding crystallite size of different CuBi2O4 thin films deposited on an FTO substrate. | ||
The crystallite size of the CuBi2O4 thin films was estimated based on the (211) facets using the Scherrer equation,41 as presented in Fig. 1b, having values of 1137, 998, 730, and 997 Å for 0 g-urea-6C-9mM, 0.1 g-urea-6C-9mM, 0.15 g-urea-6C-9mM, and 0.2 g-urea-6C-9mM, respectively. The variation in crystallite size can be attributed to the supersaturation-controlled nucleation and growth kinetics.42,43 Although all samples contain identical Bi3+ and Cu2+ in 100 mL EG, the effective instantaneous supersaturation at the electrode is modulated by urea content.42,43 Urea coordinates with Cu2+/Bi3+ and interacts with the solvent through hydrogen bonding, controlling the rate of ion release during electrodeposition.31,32 With 0.15 g urea, this produces a sustained supersaturation, promoting high nucleation density and restricting individual crystallite growth, resulting in the smallest crystallite size (730 Å).42,43 Without urea, rapid ion release generates a brief supersaturation spike; although the nucleation barrier is low during the spike, its short duration limits the number of nuclei, leading to large crystallites (1137 Å).42,43 At higher urea content (0.2 g), increased viscosity and slower ion diffusion again reduce nucleation, allowing the crystallite to grow larger (997 Å).42,43
The surface morphologies of the CuBi2O4 thin films were examined using SEM, with representative low-magnification images shown in Fig. S3 and high-magnification images in Fig. 2. All CuBi2O4 thin films exhibit uneven and rough surface features, with some regions containing small nanoparticles, while in others the nanoparticles coalesce into larger branched structures, forming interconnected networks. These networks display irregular, coral-like morphologies reminiscent of finger coral, characterized by pronounced variations in particle size, spatial distribution, and height. The size differences of nanoparticles among the four CuBi2O4 thin films were estimated from the particle diameter distribution histograms shown in Fig. 2. Specifically, the average diameter of the sample without urea (0 g-urea-6C-9mM) is 0.15 µm, whilst being 0.21 and 0.17 µm, respectively, for 0.1 g-urea-6C-9mM and 0.2 g-urea-6C-9mM, the largest average diameter of 0.24 µm being achieved for sample 0.15 g-urea-6C-9mM. This behavior can be understood by combining classical nucleation and growth with non-classical aggregation mechanisms:42,43 small crystallites with high surface energy tend to aggregate, forming larger particles. With 0.15 g urea, moderated ion release sustains supersaturation, generating many small crystallites that aggregate in a controlled manner, producing well-packed particles and uniform films. Without or with excess urea, nucleation is limited and aggregation is slowed, resulting in less uniform particles. The low-resolution SEM images in Fig. S3 further highlight the differences in coral-like interconnected networks distribution, with the urea-containing thin films, especially 0.1 g-urea-6C-9mM and 0.15 g-urea-6C-9mM thin films, showing the most pronounced branched structures.
The cross-sectional SEM images in Fig. 3 reveal that all CuBi2O4 thin films are composed of nanoparticle aggregates coalesced into porous, irregular networks. Their vertical height of these branched regions varies, and the mean thickness was determined from the thickness histogram in Fig. 3. The thin film prepared under urea-free conditions (0 g-urea-6C-9mM) exhibits a moderately compact structure with a mean thickness of ∼1.1 µm. In contrast, films deposited from urea-containing electrolytes display increased thicknesses of ∼1.5 µm (0.1 g-urea-6C-9mM), ∼1.4 µm (0.15 g-urea-6C-9mM), and ∼1.5 µm (0.2 g-urea-6C-9mM). These results show that urea promotes vertical growth by moderating ion release during electrodeposition, which enhances nucleation and aggregation along the vertical direction, leading to thicker, more pronounced branched structures. EDS spectra displayed in Fig. S4 show that Cu, Bi, O, Sn, Si, and the coated C are present in all thin films. The Cu/Bi atomic ratios of all samples, derived from EDS analysis, are 3.09
:
6.39 for 0 g-urea-6 C-9 mM, 3.49
:
7.62 for 0.1 g-urea-6C-9mM, 9.86
:
19.48 for 0.15 g-urea-6C-9mM, and 0.97
:
1.94 for 0.2 g-urea-6C-9mM. These values correspond closely to a Cu
:
Bi ratio of approximately 1
:
2, consistent with the stoichiometric composition of CuBi2O4. This result demonstrates the successful synthesis of CuBi2O4 thin films, which is consistent with the XRD analysis.
Three-dimensional (3D) perspective topographical pictures displayed in Fig. 4 present the vertical height variations caused by the bumps, pits, or other surface irregularities. The topologies of all the CuBi2O4 thin films demonstrate non-uniform growth, characterized by regions that develop laterally large and vertically tall features, while other areas exhibit minimal growth. Specifically, thin films of 0.1 g-urea-6C-9mM and 0.15 g-urea-6C-9mM show more uniform distributions of laterally large and vertically tall features compared with thin films 0 g-urea-6C-9mM and 0.2 g-urea-6C-9mM. The two-dimensional (2D) views, showing the height signals along the given section lines, the corresponding mean roughness (Ra), the maximum height (Rmax) values, and other related parameters, are presented in Fig. S5. Here Rmax is defined as the height difference between the highest and lowest points on the cross-sectional profile relative to the centre line over the length of the profile. The variations in the Rmax values indicate the non-uniform topology features of the thin films. The surface average Ra values of the CuBi2O4 thin films are 181.1 nm for sample 0 g-urea-6C-9mM, 174.9 nm for sample 0.1 g-urea-6C-9mM, sample 0.15 g-urea-6C-9mM for 154.2 nm, and 185.1 nm for 0.2 g-urea-6C-9mM (Fig. S5), indicating a more uniform surface for samples 0.1 g-urea-6C-9mM and 0.15 g-urea-6C-9mM.
![]() | ||
| Fig. 4 AFM 3D perspectives of (a) 0 g-urea-6C-9mM, (b) 0.1 g-urea-6C-9mM, (c) 0.15 g-urea-6C-9mM, and (d) 0.2 g-urea-6C-9mM. | ||
According to the XPS survey scan results shown in Fig. S6, both samples have similar entire spectra with Cu 2p, Bi 4f, and O 1s peaks. Fig. 5a reveals the high-resolution Bi 4f spectra of samples 0 g-urea-6C-9mM and 0.15 g-urea-6C-9mM. The binding energies at 158.6 and 163.9 eV are responsible for the typical spin–orbit split Bi 4f7/2 and Bi 4f5/2 peaks, respectively, which are brought on by the Bi3+ component in both samples.44,45 The small peaks fitted at lower binding energies of 162.2 and 156.9 eV are attributed to residual metallic Bi on the surface of the CuBi2O4 lattice.46,47 In Fig. 5b, the raw data exhibit four peaks: a doublet for Cu 2p and two additional shake-up peaks. These shake-up peaks are attributed to a strong interaction of the final states, involving charge transfer from O 2p band to the Cu 3d band.48,49 Cu 2p spectra are typically featured by focusing on the Cu 2p3/2 peak and the lower binding energy shake-up peak, owing to their distinctive spectral features that allow for reliable fitting. Peaks with binding energies of 933.9 and 932.3 eV, fitted from the Cu 2p3/2 peak, are attributed to the Cu2+ in the CuBi2O4 lattice and Cu+ species, respectively.44,49 Peaks located at 943.7 and 941.2 eV, obtained from the lower binding energy shake-up peak, further suggest the coexistence of these Cu2+ and Cu+ species within the CuBi2O4 thin films. The peak area ratios of Bi3+ and Bi (99
:
1), as well as those of Cu2+ and Cu+ (8
:
2), in the 0.15 g-urea-6C-9mM thin film and 0 g-urea-6C-9mM thin film, are maintained without significant change. This indicates that the addition of urea into the electrolyte during the electrodeposition process does not significantly affect the chemical composition and the chemical states of the CuBi2O4 thin film. To further confirm the chemical states and quantify the surface concentration of Cu species in the CuBi2O4 thin film, the Cu LM2 Auger spectra were analysed. As depicted in Fig. 5c, the Cu LM2 Auger electron spectra show peaks corresponding to the lattice Cu2+ (blue peaks) and the low-valence Cu+ species (red peaks).50,51 This finding is consistent with the results obtained from the Cu 2p spectra. The O 1s spectra (Fig. 5d) for both CuBi2O4 samples were fitted into three distinct peaks at binding energies of approximately 529.4, 530.5, and 531.4 eV. These peaks are attributed to lattice oxygen (OL),50–52 surface hydroxyl group,50,51 and absorbed water species, respectively. The presence of surface hydroxyl groups in CuBi2O4 is commonly associated with surface oxygen environments influenced by the coexistence of Cu2+/Cu+ states, as observed in the Cu LM2 Auger spectra. The OL peak area ratios for CuBi2O4 thin films of 0 g-urea-6C-9mM and 0.15 g-urea-6C-9mM are 49.5% and 47.2%, respectively, while the surface hydroxyl species for these two CuBi2O4 thin films are 39.3% and 42.1%, respectively. This variation may be attributed to the different distribution of aggregated, finger coral-like structures on the surface of the thin films.
![]() | ||
| Fig. 5 XPS spectra of prepared CuBi2O4 thin films: 0 g-urea-6C-9mM and 0.15 g-urea-6C-9mM. (a) Bi 4f, (b) Cu 2p, (c), Cu LM2, and (d) O 1s. | ||
The UV-visible absorbance spectra (Fig. 6a) exhibit strong visible-light absorptions in all thin films. The incident light absorbances of the CuBi2O4 thin films, 0.15 g-urea-6C-9mM and 0.2 g-urea-6C-9mM, are higher compared with those of 0 g-urea-6C-9mM and 0.1 g-urea-6C-9mM CuBi2O4. Slight red shifts can be observed upon introducing 0.15 g and 0.2 g urea into the electrodeposition process. This increase in optical absorption and reduction in reflectance (Fig. 6b, inset) are attributed to the varying surface topographical features induced by the urea. Fig. 6b presents the direct transition band gap values for the samples—1.88 eV for sample 0 g-urea-6C-9mM, 1.89 eV for sample 0.1 g-urea-6C-9mM, 1.87 eV for sample 0.15 g-urea-6C-9mM, and 1.86 eV for sample 0.2 g-urea-6C-9mM. These values are consistent with those reported in previous studies.13
![]() | ||
| Fig. 6 (a) UV-visible absorbance spectra of the CuBi2O4 thin films and (b) the corresponding predicted band gap values from the Kubelka–Munk function,53 with the UV-visible diffuse reflectance spectra (DRS) shown in the inset. | ||
Mott–Schottky measurements are utilized at 1 kHz in 0.1 M NaOH solution in dark conditions, to estimate the flat band potential (Efb) and the hole density (NA) of the CuBi2O4 thin films. The Mott–Schottky plots depicted in Fig. 7a were obtained from the following eqn (2):54
![]() | (2) |
| Ec = Ev − Eg | (3) |
| Thin film | Conduction band potential (V vs. RHE) | Valence band potential (V vs. RHE) | Value of the slope | N A (cm−3) |
|---|---|---|---|---|
| 0 g-urea-6C-9mM | −0.47 | 1.41 | −3.18 × 1010 | 5.54 × 1019 |
| 0.1 g-urea-6C-9mM | −0.48 | 1.41 | −5.94 × 1010 | 2.97 × 1019 |
| 0.15 g-urea-6C-9mM | −0.47 | 1.40 | −3.57 × 1010 | 4.94 × 1019 |
| 0.2 g-urea-6C-9mM | −0.46 | 1.40 | −3.46 × 1010 | 5.09 × 1019 |
The potentials of both the conduction and valence band positions relative to the water redox potential were constructed and illustrated in Fig. 7b. The valence band potentials of all the thin films at pH 12.8 are above 1.98 V vs. RHE (corresponding to 1.23 V vs. RHE at pH 0), while the conduction band potentials are above 0.75 V vs. RHE (corresponding to 0 V vs. RHE at pH 0).58 Therefore, the CuBi2O4 thin films demonstrate the capability to perform the hydrogen evolution reaction under a pH of 12.8. Additionally, the band edge positions differ by only 10 mV, which is within the range of experimental error and would not significantly affect the overall band structure. The ultraviolet photoelectron spectroscopic (UPS) analysis using a He(I) laser (hv = 21.22 eV) as the incident light was utilized to determine the work functions of the CuBi2O4 thin films. As illustrated in Fig. S7, the work functions for samples 0 g-urea-6C-9mM and 0.15 g-urea-6C-9mM are 5.03 eV and 5.21 eV, respectively. Therefore, the energy required to extract electrons from the materials differs by a maximum of 0.2 eV, which is quite similar and might not be enough to affect the PEC water splitting performance.
The J–V and Nyquist plots are displayed in Fig. 8, which were obtained from linear sweeping voltammetry (LSV) and frequency response analyser (FRA) impedance potentiostatic mode, respectively, for the investigation of the PEC performance of the CuBi2O4 thin films. In Fig. 8a, sample 0 g-urea-6C-9mM exhibits the highest photocurrent density of −0.94 mA cm−2 at 0.52 V vs. RHE, whilst the photocurrent density was further improved after the introduction of urea during the electrodeposition process, reaching −1.04 and −1.44 mA cm−2 for the 0.1 and 0.15 g urea samples, which is higher than previously published values (Table S2). Further increasing the amount of urea to 0.2 g led to a decreased photocurrent density of −0.99 mA cm−2, which is still greater than that of sample 0 g-urea-6C-9mM. Therefore, the charge separation and transfer efficiencies have been improved. The Nyquist spectra (Fig. 8b) fitted via the following equivalent circuit inserted in the Fig. 8b are applied to investigate the charge-transfer resistance of the CuBi2O4 thin films. Rs represents the intercept of the semicircle with the real axis (Re(Z)) at low frequency, indicating the resistance of the electrolyte solution and any inherent resistance of the setup. The diameter of the semicircle in the Nyquist plot is utilized to determine the charge transfer resistance (Rct) at the interface of the CuBi2O4 thin film and FTO glass. The double-layer capacitance (Cdl) represents the capacitive behaviour at the interface, while the Warburg impedance (W) indicates the influence of the diffusion processes on the impedance response in combination with the double-layer capacitance. According to Fig. 8b and Table 2, sample 0.15 g-urea-6C-9mM, which exhibits the smallest semicircle diameter, has the lowest charge transfer resistances (Rct) with a value of 2430 Ω. This is followed by sample 0.2 g-urea-6C-9mM (Rct = 3193 Ω), and sample 0.1 g-urea-6C-9mM (Rct = 4169 Ω), whilst the urea-free reference CuBi2O4 thin film exhibits the highest charge transfer resistance of 5441 Ω. Therefore, the charge transfer efficiencies of the CuBi2O4 thin films are ranked as follows: 0.15 g-urea-6C-9mM > 0.2 g-urea-6C-9mM > 0.1 g-urea-6C-9mM > 0 g-urea-6C-9mM under the bias of 0.52 V vs. RHE. The Cdl values, however, exhibit an opposite trend, with the 0.15 g-urea-6C-9mM thin film showing the highest Cdl value. This implies a larger area for charge separation, potentially enhancing the PEC performance. The observed decrease in charge transfer resistance, coupled with the increased Cdl, indicates an improved efficiency in utilizing photogenerated charges for these thin films.
| 0 g-urea-6C-9mM | 0.1 g-urea-6C-9mM | 0.15 g-urea-6C-9mM | 0.2 g-urea-6C-9mM | |
|---|---|---|---|---|
| R s | 26.22 Ω | 23.06 Ω | 20.75 Ω | 25.31 Ω |
| R ct | 5441 Ω | 4169 Ω | 2430 Ω | 3193 Ω |
| C dl | 15.26 µF | 16.58 µF | 18.29 µF | 15.25 µF |
| W | 14 593 Ω s−0.5 |
9099 Ω s−0.5 | 3925 Ω s−0.5 | 6932 Ω s−0.5 |
The stability of the CuBi2O4 thin films has been investigated using I–t curves measured through chronoamperometry under continuous solar light illumination for 2 h at 0.52 V vs. RHE. Thin film 0.15 g-urea-6C-9mM, which exhibits the highest photocurrent density, was used to investigate the stability performance in comparison to thin film 0 g-urea-6C-9mM. As shown in Fig. S8, both thin films showed poor PEC stabilities. Specifically, initial sharp drops in photocurrent density were observed within approximately 10 seconds. After 100 seconds of continuous illumination, both films reached their steady-state photocurrents, with sample 0.15 g-urea-6C-9mM exhibiting approximately −0.27 mA cm−2, which is about 42% higher than the reference (−0.19 mA cm−2).
The variations in PEC performance (photocurrent density and photostability) of the CuBi2O4 thin films can be primarily attributed to the urea-induced morphological differences. Examination of the charge–time curves (Fig. S9) recorded during electrodeposition reveals that the presence of urea slows the deposition process under the same applied charges. This effect arises from the coordination of urea with Cu2+ and Bi3+ ions, which stabilizes them in solution and maintains a steady concentration of electroactive species.59,60 This behavior further supports the sustained-supersaturation condition proposed earlier, which governs both the nucleation rate and subsequent particle aggregation. Under the influence of urea, the resulting more continuous particle networks and slightly increased film thickness provide more efficient charge-transport pathways and reduce recombination at interparticle boundaries, accounting for the higher photocurrent densities of the urea-assisted thin films. Moreover, the improved structural coherence and uniformity help to moderate photocorrosion, explaining the modest enhancement in stability compared to the referent thin film (0 g-urea-6C-9mM).
Based on above analyses, a possible charge separation and transfer mechanism has been proposed to explain the improved PEC water splitting performance. As illustrated in Fig. 9, under solar light illumination, the CuBi2O4 thin film generates photocarriers after absorbing light with energy greater than its band gap. A built-in electric field separates the photo-generated electrons and holes. Ideally, the electrons will transfer to the surface of the CuBi2O4 thin film to participate in the reduction reaction eqn (4):61
| 4H2O + 4e− → 2H2 + 4−OH | (4) |
Meanwhile, the Pt counter electrode serves as the site for the oxidation reaction eqn (5):61
| 4−OH → O2 + 2H2O + 4e− | (5) |
The generated electrons at the counter electrode flow through the external circuit back to the CuBi2O4 thin film, completing the electrical circuit. However, an undesirable photocorrosion reaction occurs within the CuBi2O4 thin film due to the insufficiently anodic valence band potential of the synthesized CuBi2O4 photocathode (1.4 V vs. RHE at pH 12.8). This thermodynamic limitation prevents the efficient oxygen evolution reaction (OER), leading to the photogenerated holes oxidizing the CuBi2O4 lattice and degrading the Cu–O–Bi framework into CuO and Bi2O3 rather than driving water oxidation (1.98 V vs. RHE at pH 12.8).62–64 Furthermore, in an alkaline electrolyte, Bi3+ within the CuBi2O4 lattice reacts with hydroxide ions to form soluble hydroxo complexes,62 leading to the bismuth leaching effect and structural destabilization. Simultaneously, Cu+ species undergo the following disproportionation eqn (6):63
| 2Cu+ → Cu2+ + Cu0 | (6) |
| This journal is © The Royal Society of Chemistry 2026 |