Open Access Article
Samala Venkateswarlua,
Ravinder Singha,
Songbo Cuia,
Andrew Stellaa,
Jimmy Papazotosa,
Naixin Zhaoa,
Yi Yuana,
Haitao Liub,
Xu Lib,
Jinliang Wangb and
Yuning Li
*a
aDepartment of Chemical Engineering and Waterloo Institute for Nanotechnology (WIN), University of Waterloo, Waterloo, Ontario N2L 3G1, Canada. E-mail: yuning.li@uwaterloo.ca
bInstitute of Chemistry, Henan Academy of Sciences, 56 Hongzhuan Road, Jinshui District, Zhengzhou, Henan 450002, China
First published on 20th April 2026
Conjugated polymers with strong absorption in the near-infrared-II (NIR-II, 1000–1700 nm) region are highly attractive for optoelectronic and biomedical applications. However, their availability remains limited, and their synthesis often relies on costly transition-metal catalysts and multistep monomer functionalization. Here, we report an economical and operationally simple FeCl3-mediated oxidative polymerization strategy for the synthesis of two alkoxythiophene-flanked diketopyrrolopyrrole (DPP) homopolymers, PDPPC20OT-C8 and PDPPC12OT-C12. Both polymers exhibit intense NIR-II absorption, with maximum absorption wavelengths reaching ∼1190 nm, and ultralow optical band gaps as low as 0.87 eV. Side-chain engineering enables good solubility, suppressed crosslinking, and more coplanar backbone conformation, leading to ambipolar charge transport in organic thin-film transistors. Optimized devices based on PDPPC12OT-C12 show balanced hole and electron mobilities of up to 1.74 × 10−3 and 2.24 × 10−3 cm2 V−1 s−1, respectively. Density functional theory calculations reveal enhanced backbone planarity arising from intramolecular S⋯O interactions, consistent with the experimentally observed red-shifted absorption and reduced band gaps. This work establishes FeCl3-mediated oxidative polymerization as a viable route to high-performance ultralow-band-gap DPP homopolymers and highlights their potential for ambipolar thin film transistors, NIR-II photodetection, and photothermal applications.
Thiophene-flanked diketopyrrolopyrrole (DPPT) is one of the most widely used electron-acceptor building blocks for constructing D–A polymers.26–31 Coupling DPPT with strong donor comonomers, such as alkoxythiophenes, effectively narrows the band gap and redshifts the λmax to approximately 900 nm.25 However, this absorption remains within the NIR-I region. An alternative and more effective strategy is to pair DPPT with strong electron-acceptor comonomers. Benzobisthiadiazole (BBT),30 thienoisoindigo,32 and dioxothienopyrrolebenzodifurandione (BTPBF)33 have each been demonstrated to significantly extend polymer absorption into the NIR-II region. Stille and Suzuki C–C coupling polymerization methods using Pd catalysts have been used to synthesize these polymers.
More recently, efforts have made to replace the flanking thiophenes of DPP with the more electron-rich 4-methoxythiophene to further reduce the optical band gap. Domokos et al.34 synthesized a series of homo- and copolymers based on 2-hexyldecyl-substituted DPP flanked by 4-methoxythiophene (DPPC1OT) units and observed pronounced redshifts relative to their thiophene-flanked counterparts. The homopolymer PDPP-2T-OMe exhibited λmax values of 1019 nm in solution and 998 nm in thin film. Giri et al.35 reported a structurally related homopolymer (P3) derived from DPPC1OT bearing longer 2-octyldodecyl side chains on the DPP unit, which displayed λmax values of ∼1000 nm in both solution and thin films and showed promise as an energy-efficient infrared electrochromic optical attenuator. Collectively, these studies indicate that homopolymers derived from DPPC1OT are promising semiconductor materials for NIR-II photodetection. However, their charge transport properties have not yet been investigated. In addition, both studies relied on Yamamoto coupling of dibrominated monomers using bis(1,5-cyclooctadiene)nickel(0) (Ni(COD)2), which is costly and suffers from limited thermal and chemical stability. Consequently, the development of a more economical and operationally simple polymerization strategy remains highly desirable.
In this study, we first attempted to synthesize the homopolymer PDPPC1OT-C20 (Scheme 1), which corresponds to P3 reported by Giri et al.,35 via FeCl3-mediated oxidative polymerization of the non-brominated monomer M2. However, the resulting polymer was completely insoluble in common organic solvents, which is primarily attributed to the crosslinking side reactions during polymerization. To improve solubility, the long, branched 2-octyldodecyl side chains were relocated from the DPP core to the flanking thiophenes to replace the methyl substituents, while the DPP unit was substituted with 1-octyl side chains to afford a redesigned monomer, M3. The corresponding polymer, PDPPC20OT-C8, prepared by FeCl3 oxidative polymerization, exhibited good solubility in chloroform, chlorobenzene, and 1,2-dichlorobenzene. In addition, a second redesigned monomer, M4, a 4-dodecyloxythiophene-flanked DPP bearing 2-butyloctyl side chains, was synthesized. Polymerization of M4 under identical conditions yielded another soluble homopolymer, PDPPC12OT-C12.
PDPPC20OT-C8 exhibits λmax values of 1035 nm in solution and 1190 nm in thin films, whereas PDPPC12OT-C12 shows λmax values of 1176 nm in solution and 1186 nm in thin films. All of these absorption maxima fall within the NIR-II regime. Notably, these λmax values are substantially redshifted relative to previously reported DPPC1OT homopolymers,34,35 indicating more extended π-conjugation. This enhancement is likely associated with improved backbone planarity induced by long alkoxy side chains on the thiophene units, which reduce backbone rotational freedom. When used as the semiconducting channel layers in OTFTs, both polymers display ambipolar charge transport with well-balanced hole and electron mobilities reaching approximately 10−3 cm2 V−1 s−1.
In contrast, introduction of electron-donating alkoxy substituents on the thiophene units in M2–M4 lowers the oxidation potential and stabilizes the cation radical during oxidative polymerization.42,43 Indeed, M2 underwent rapid polymerization with FeCl3 at room temperature, yielding a dark bluish-green solid in 95% yield after purification (Fig. S4). However, the resulting PDPPC1OT-C20 is completely insoluble in common organic solvents. Even after heating in 1,2,4-trichlorobenzene (TCB) under stirring at temperatures up to 200 °C for 8 h, the polymer remained largely undissolved (Fig. S5). Since the same polymer prepared by Yamamoto coupling was reported to be readily soluble in chloroform at room temperature and TCB at 120 °C, the poor solubility of PDPPC1OT-C20 is likely due to the formation of crosslinked structures via coupling at the 3-position of the thiophene units.44 This interpretation is further supported by the highly disordered structure of this polymer observed by powder X-ray diffraction (Fig. S6).
To improve solubility and suppress potential crosslinking, monomer M3 was designed with bulky branched 2-octyldodecyl side chains on the thiophene units and linear octyl side chains on the DPP core. These sterically demanding substituents were expected to enhance solubility and reduce undesired side reactions. Compared with M2, M3 exhibited reduced reactivity, yielding only oligomeric products after 24 h at room temperature. Therefore, the polymerization temperature was increased to 50 °C and the reaction time was extended to 48 h. Under these conditions, the target polymer PDPPC20OT-C8 was obtained as a dark bluish-green solid (Fig. S4) in a yield of 70% after purification. This polymer is readily soluble in chloroform, chlorobenzene, and 1,2-dichlorobenzene. High-temperature gel permeation chromatography (HT-GPC) revealed a weight-average molecular weight (Mw) of 54.3 kDa and a dispersity (Đ) of 3.17 (Fig. S34 and Table 1). The good solubility and relatively narrow molecular weight distribution suggest minimal crosslinking during polymerization.
| Homo polymer | Mna [kDa] | Mwa [kDa] | Đ | Tdb [°C] | λmax [nm] | λonset [nm] | Eoptg c [eV] |
EHOMOd [eV] | ELUMOd [eV] | |
|---|---|---|---|---|---|---|---|---|---|---|
| Solution | Film | Film | Film | Film | Film | |||||
| a Measured from HT-GPC using 1,2,4-trichlorobenzene as eluent at 150 °C.b Thermal decomposition temperatures (Td) at 5% weight loss.c Optical band gaps (Eopt'sg) calculated from the onset absorption wavelengths of the thin films.d EHOMO of polymers were calculated using their onset oxidation potentials, with ferrocene as a reference (EHOMO = −4.8 eV), while ELUMO's of polymers were calculated using the equation: ELUMO = Eoptg + EHOMO. | ||||||||||
| PDPPC20OT-C8 | 17.1 | 54.3 | 3.17 | 333 | 1035 | 1190 | 1375 | 0.90 | −5.06 | −4.16 |
| PDPPC12OT-C12 | 8.7 | 25.3 | 2.91 | 350 | 1176 | 1186 | 1425 | 0.87 | −5.02 | −4.15 |
Monomer M4, bearing less bulky linear dodecyl side chains on the thiophene units and 2-butyloctyl substituents on the DPP core, was also polymerized under identical conditions at 50 °C for 48 h. The corresponding polymer, PDPPC12OT-C12, was obtained in a comparable yield of 71% after purification. This polymer exhibits a lower Mw of 25.3 kDa with a dispersity of 2.91 and shows good solubility in chloroform, chlorobenzene, and 1,2-dichlorobenzene. These results indicate that long alkoxy-substituted thiophene units are effective in promoting oxidative polymerization while simultaneously suppressing crosslinking side reactions.
Thermo-gravimetric analysis (TGA) revealed decomposition temperatures of 333 °C for PDPPC20OT-C8 and 350 °C for PDPPC12OT-C12 (Fig. S1), confirming their good thermal stability.
The large red shift in λmax and the reversal in the relative intensities of the 0–0 and 0–1 bands indicate a substantial change in chain conformation and intermolecular organization from solution to the solid state. In solution, the relatively strong 0–1 transition suggests that the polymer chains are more conformationally relaxed and less planar due to good solvation by chloroform and the steric effect of the long branched C20 side chains. Upon film formation, solvent removal promotes closer interchain packing and enhanced π–π interactions, while also restricting torsional motion and increasing backbone planarization. As a result, the effective conjugation length increases, leading to the marked bathochromic shift. In addition, the dominance of the 0–0 transition in the film and the 0–1/0–0 intensity ratio below unity are consistent with J-aggregate-like character,45,46 further supporting stronger interchain donor–acceptor interactions and a more ordered solid-state packing structure.
In contrast, PDPPC12OT-C12 exhibits λmax values of 1176 nm in solution and 1186 nm in the solid state, both of which are attributed to the 0–0 transition, with only a small redshift of 10 nm upon film formation. Unlike PDPPC20OT-C8, the solution spectrum of PDPPC12OT-C12 is dominated by the 0–0 transition, indicating that significant aggregation is already present in solution. This behavior is consistent with the stronger aggregation tendency of PDPPC12OT-C12, likely arising from its more coplanar backbone and the reduced steric hindrance provided by the linear dodecyloxy side chains on the thiophene units together with the shorter branched side chains on the DPP core. These structural features are less effective at preventing interchain association than the longer branched side chains in PDPPC20OT-C8. As a result, PDPPC12OT-C12 likely exists as pre-aggregated species in solution, as commonly reported for DPP-based polymers with strong intermolecular interactions.47 Therefore, its solution absorption spectrum already resembles that of the solid film, and only a minimal additional red shift is observed upon film formation.
The optical band gaps of PDPPC20OT-C8 and PDPPC12OT-C12 were estimated to be 0.90 eV and 0.87 eV, respectively, based on the onset absorption wavelenghts of their thin films. Compared with previously reported homopolymers of 4-methoxythiophene-flanked DPP (DPPC1OT),34,35 both polymers exhibit substantially redshifted absorption maxima and narrower optical band gaps. The shorter λmax (998 nm) reported by Domokos et al.34 may be attributed to the low molecular weight of the polymer (Mw = 2.9 kDa). However, even the higher-molecular-weight polymer (Mw = 87.7 kDa) reported by Giri et al.35 displayed a λmax of only ∼1000 nm. These observations suggest that replacing the methyl substituents on the methoxythiophene units with longer linear or branched alkyl chains significantly extends the effective conjugation length in DPP homopolymers. Such long alkyl side chains can pack more tightly in the solid state, restricting the rotational freedom of thiophene units along the polymer backbone and promoting a more coplanar conformation, whereas smaller methyl substituents are more mobile and contribute less effectively to backbone planarity.48
Most CPs exhibited λmax either in visible (400–780 nm) or NIR I (780–1000 nm) regions, whereas polymers with λmax in NIR-II region remain rare.25 Notably, in this work, both newly prepared polymers display absorption maxima in the NIR-II window. Even more importantly, these polymers exhibit very weak π–π* absorption in the visible region compared with their strong absorption in the NIR-I and NIR-II regions, resulting in high optical transparency in the visible range. Such selective infrared absorption combined with high visible transparency is highly desirable for a range of applications, including NIR photodetectors, infrared camouflage, transparent laser welding, thermal-shielding coatings on glass, and transparent electronic devices.25,49–52
Cyclic voltammetry (CV) measurements were conducted on thin films of PDPPC20OT-C8 and PDPPC12OT-C12 deposited on platinum electrodes (Fig. 1b). The HOMO energy levels of PDPPC20OT-C8 and PDPPC12OT-C12 were calculated to be −5.06 eV and −5.02 eV, respectively, based on their onset oxidation potentials. Despite the different side chain lengths, the two polymers exhibit very similar HOMO energy levels because these are governed mainly by the electronic structure and effective conjugation length of the same conjugated backbone rather than by the alkyl side chains. The slightly deeper HOMO level of PDPPC20OT-C8 may arise from side-chain-induced steric effects that increase backbone twisting and less ordered chain packing, consistent with its slightly larger optical bandgap (0.90 eV vs. 0.87 eV).
Owing to the absence of discernible reduction peaks, their LUMO energy levels were derived by subtracting the optical band gaps (Eoptg) from the HOMO energy levels, yielding values of −4.03 and −4.05 eV, respectively. It should be noted that the LUMO energy levels estimated in this way may be significantly lower than the actual values because organic semiconductors typically exhibit large exciton binding energies, often around 0.3 eV or higher.53,54
The dihedral angle between the thiophene and DPP units (C–C–C–S), defined as θ1, and the dihedral angle between adjacent thiophene units in the bithiophene, defined as θ2, were used to evaluate the backbone conformation. The monomer models, DPPT and DPPC1OT, exhibit θ1 values of 0.14° and 0.88°, respectively, indicating that both structures are highly coplanar, with DPPT being slightly more planar. Compared with DPPT, DPPC1OT shows elevated frontier energy levels due to the strong electron-donating effect of the methoxy substituent. In particular, the HOMO level increases by 0.14 eV, which facilitates the significantly enhanced oxidative polymerization of M2–M4 relative to M1. In addition, the HOMO–LUMO gap of DPPC1OT decreases by 0.10 eV compared to DPPT primarily due to the significantly elevated HOMO energy level.
The dimer [DPPT]2 exhibits θ1 and θ2 values of 7.2° and 6.5°, respectively, which are substantially larger than those of [DPPC1OT]2 (θ1 = 1.3° and θ2 = 0.17°). The markedly reduced θ2 in [DPPC1OT]2 is attributed to conformational locking of the bithiophene segment induced by intramolecular S⋯O interactions55 at a distance of 2.81 Å. In contrast, such interactions are absent in [DPPT]2, resulting in a more twisted backbone. Notably, θ1 in [DPPC1OT]2 increases only slightly relative to the monomer, whereas [DPPT]2 shows a pronounced increase in θ1 relative to monomer, indicating a strong disruption of backbone coplanarity upon chain extension. The enhanced coplanarity of [DPPC1OT]2 promotes more effective π-conjugation and contributes to band-gap narrowing in the alkoxy-substituted polymers.
The calculated HOMO and LUMO energy levels are −4.83 and −2.87 eV for [DPPT]2 and −4.49 and −2.73 eV for [DPPC1OT]2, respectively, corresponding to band gaps of 1.96 and 1.76 eV. These results indicate that the incorporation of alkoxy substituents significantly raises the HOMO energy level and narrows the band gap. The computational trends are consistent with the experimentally observed electrochemical and optical properties of PDPPC20OT-C8, PDPPC12OT-C12, and PDPPT-C20, as well as related analogues bearing different alkyl side chains on the DPP units.34,35,56
For PDPPC20OT-C8, both the out-of-plane and in-plane patterns show no discernible diffraction peaks in the as-cast state or after annealing up to 200 °C, indicating a predominantly amorphous microstructure. Upon annealing at 250 °C, a weak lamellar (100) reflection appears in the out-of-plane direction at 2θ = 6.73°, corresponding to a d-spacing of 13.1 Å, along with a broad π–π stacking (010) peak at 2θ = 23.8° (d = 3.7 Å). These features suggest the onset of limited molecular ordering only at elevated temperatures.
In contrast, PDPPC12OT-C12 exhibits a distinct lamellar (100) peak at 2θ = 8.0°, corresponding to a d-spacing of 11.1 Å, together with a broad (010) peak at 2θ = 23.8° (d = 3.7 Å) in the out-of-plane pattern of the as-cast film. After annealing at 100 °C, the lamellar peak shifts slightly to 2θ = 7.6° (d = 11.6 Å) and increases in intensity, indicating enhanced crystallinity. Further annealing at 150, 200, and 250 °C induces no significant additional changes in the diffraction patterns, suggesting that the film microstructure is largely stabilized after mild thermal treatment.
Compared with DPPT-based homopolymers bearing unsubstituted thiophene units, which typically exhibit high crystallinity,56 the present results indicate that alkoxy substitution on the thiophene units hinders efficient chain packing. In particular, the long branched 2-octyldodecyloxy side chains in PDPPC20OT-C8 substantially suppress long-range order, whereas the linear octyloxy side chains in PDPPC12OT-C12 partially preserve lamellar ordering.
To evaluate the surface morphology of the polymer thin films, tapping-mode atomic force microscopy (AFM) measurements were performed, and height images were collected for both polymers in the as-cast state and after thermal annealing (Fig. 4). PDPPC20OT-C8 exhibits root mean-square (RMS) roughness values of 1.01 nm for the as-cast film and 0.91 nm for the film after annealing at 200 °C. In comparison, PDPPC12OT-C12 shows RMS roughness values of 0.65 nm in the as-cast state and 0.71 nm after annealing at 250 °C. The lower RMS roughness of PDPPC12OT-C12 indicates that it forms a smoother film surface than PDPPC20OT-C8. Nonetheless, both homopolymers exhibit RMS roughness values of 1 nm or less, indicating the formation of uniform and smooth films, which is favorable for improving interfacial contact between the dielectric and semiconducting layers in organic thin-film transistor devices.
![]() | ||
| Fig. 4 AFM height images (1 µm × 1 µm) of (a) PDPPC20OT-C8 (as-cast) (b) PDPPC20OT-C8 (annealed at 200 °C) (c) PDPPC12OT-C12 (as-cast) (d) PDPPC12OT-C12 (annealed at 250 °C). | ||
Representative output and transfer curves are shown in Fig. 5, and the extracted charge carrier mobilities in the saturation regime are summarized in Table 2. Both PDPPC20OT-C8 and PDPPC12OT-C12 exhibit typical p-type and n-type transistor behavior at their respective optimal annealing temperatures, indicating ambipolar charge transport.
![]() | ||
| Fig. 5 Representative transfer and output characteristics of OTFT devices fabricated from (a) and (c) PDPPC20OT-C8 after annealing at 200 °C and (b) and (d) PDPPC12OT-C12 after annealing at 250 °C. | ||
| Homo polymer | Annealing temperature [°C] | Max. µha [cm2 V−1 s−1] | Avg. µhb ± Std [cm2 V−1 s−1] | Max. µea [cm2 V−1 s−1] | Avg. µeb ± Std [cm2 V−1 s−1] | Ion/Ioffc hole/electron | Vthd [V] hole/electron |
|---|---|---|---|---|---|---|---|
| a The maximum hole/electron mobility in the saturated region (Max. µh/µe).b The average hole mobility in the saturated region and standard deviation (std) from more than 5 devices.c The on-current (Ion) to off-current (Ioff) ratio.d The average threshold voltage (Vth). | |||||||
| PDPPC20OT-C8 | As-cast | 1.93 × 10−4 | 1.69 × 10−4 ± 2.4 × 10−5 | 3.32 × 10−4 | 2.63 × 10−4 ± 7.9 × 10−5 | 102/101 | 24.0/72.1 |
| 100 | 1.82 × 10−4 | 1.71 × 10−4 ± 9.3 × 10−6 | 3.56 × 10−4 | 3.19 × 10−4 ± 3.2 × 10−5 | 102/101 | 9.7/69.7 | |
| 150 | 2.95 × 10−4 | 2.61 × 10−4 ± 3.0 × 10−5 | 3.82 × 10−4 | 3.35 × 10−4 ± 4.1 × 10−5 | 102/101 | −1.8/64.9 | |
| 200 | 6.83 × 10−4 | 6.05 × 10−4 ± 1.1 × 10−4 | 5.74 × 10−4 | 5.11 × 10−4 ± 9.8 × 10−5 | 102/101 | −5.2/62.8 | |
| 250 | 5.86 × 10−4 | 5.74 × 10−4 ± 1.7 × 10−5 | 3.44 × 10−4 | 3.39 × 10−4 ± 5.4 × 10−6 | 102/101 | −4.2/61.1 | |
| PDPPC12OT-C12 | As-cast | 1.27 × 10−3 | 0.85 × 10−3 ± 3.8 × 10−4 | — | — | 102/— | 51.4/— |
| 100 | 0.86 × 10−3 | 0.82 × 10−3 ± 4.5 × 10−4 | — | — | 102/— | 23.7/— | |
| 150 | 1.23 × 10−3 | 1.17 × 10−3 ± 7.4 × 10−5 | — | — | 102/— | 16.9/— | |
| 200 | 1.57 × 10−3 | 1.43 × 10−3 ± 1.6 × 10−4 | 1.57 × 10−3 | 1.31 × 10−3 ± 3.0 × 10−4 | 101/100 | 3.2/65.6 | |
| 250 | 1.74 × 10−3 | 1.66 × 10−3 ± 1.3 × 10−4 | 2.24 × 10−3 | 2.14 × 10−3 ± 1.3 × 10−4 | 101/101 | −4.9/59.9 | |
For PDPPC20OT-C8, the device annealed at 200 °C shows the best overall performance, with a highest hole mobility of 6.83 × 10−4 cm2 V−1 s−1 (average: 6.05 × 10−4 ± 1.1 × 10−4) and a highest electron mobility of 5.74 × 10−4 cm2 V−1 s−1 (average: 5.11 × 10−4 ± 9.8 × 10−5) (Table 2). In comparison, PDPPC12OT-C12 achieves highest hole and electron mobilities of 1.74 × 10−3 cm2 V−1 s−1 (average: 1.66 × 10−3 ± 1.3 × 10−4) and 2.24 × 10−3 cm2 V−1 s−1 (average: 2.14 × 10−3 ± 1.3 × 10−4), respectively, at annealing temperature of 250 °C. No measurable electron transport is observed for PDPPC12OT-C12 in the as-cast state or after annealing at temperatures up to 150 °C, which may be attributed to the presence of a high density of electron traps in the thin films (Table 2).
Overall, PDPPC12OT-C12 exhibits higher charge carrier mobilities compared to PDPPC20OT-C8, which is is likely associated with its improved chain packing and smoother film morphology, as evidenced by XRD and AFM analyses (Fig. 3 and 4).
Gold contacts were used in this study. Their work function of approximately −5.0 eV is well aligned with the HOMO energy levels of PDPPC20OT-C8 and PDPPC12OT-C12, but poorly aligned with their LUMO energy levels, resulting in large barriers for electron injection. This may account for the large threshold voltages (Vth) observed for n-channel operation (Table 2). Using an electrode material with a higher Fermi energy level or modifying the surface of the Au electrodes58–61 may facilitate electron injection and thereby reduce the threshold voltage.
To assess device operational stability in ambient air, OTFTs based on PDPPC20OT-C8 and PDPPC12OT-C12 annealed at their optimal temperatures 200 °C and 250 °C, respectively, in the glove box were characterized under ambient conditions (22 °C, 50% relative humidity) (Fig. S7 and Table S1). While their hole-transport performance remained largely stable (µ = 1.27 × 10−3 cm2 V−1 s−1 for PDPPC20OT-C8 and 1.58 × 10−3 cm2 V−1 s−1 for PDPPC12OT-C12), their electron transport performance vanished when characterized in ambient air. As noted earlier, the actual LUMO energy levels of these polymers may be significantly higher than the values estimated by combining the HOMO energy level measured by CV with the optical bandgap, owing to their large exciton binding energies. Consequently, their actual LUMO energy levels may exceed the threshold required for air-stable electron transport, which is typically around −4.0 eV.62,63
| This journal is © The Royal Society of Chemistry 2026 |