Open Access Article
Qing Langa,
Jian Zhanga,
Ziqi Tian
a,
Qi Fanb,
Kun Liangb,
Yunan Wanga,
Evgenia Dmitrieva
c,
Liang Chen
a,
Congxue Liu*d,
Mingchao Wang
d,
Maryam AlNahyane,
Panče Naumov
*fghi and
Gang Wang
*a
aZhejiang Key Laboratory of Advanced Fuel Cells and Electrolyzers Technology, Materials Tech Laboratory for Hydrogen & Energy Storage, Ningbo Institute of Materials Technology and Engineering (NIMTE) of the Chinese Academy of Sciences (CAS), Ningbo 315201, P. R. China. E-mail: gang.wang@nimte.ac.cn
bZhejiang Key Laboratory of Data-Driven High-Safety Energy Materials and Applications, Ningbo Institute of Materials Technology and Engineering of the Chinese Academy of Sciences, Ningbo 315201, China
cLeibniz Institute for Solid State and Materials Research Dresden (IFW), Helmholtzstraße 20, 01069 Dresden, Germany
dGuangdong Provincial Key Laboratory of Nano-Micro Materials Research & State Key Laboratory of Advanced Waterproof Materials, School of Advanced Materials, Peking University, Shenzhen Graduate School, Shenzhen 518055, China. E-mail: liucongxue@pku.edu.cn
eDivision of Engineering, New York University Abu Dhabi, POB 129188, Abu Dhabi, United Arab Emirates
fSmart Materials Lab, New York University Abu Dhabi, PO Box 129188, Abu Dhabi, United Arab Emirates
gCenter for Smart Engineering Materials, New York University Abu Dhabi, PO Box 129188, Abu Dhabi, United Arab Emirates
hResearch Center for Environment and Materials, Macedonian Academy of Sciences and Arts, Bul. Krste Misirkov 2, MK-1000 Skopje, Macedonia
iMolecular Design Institute, Department of Chemistry, New York University, 100 Washington Square East, New York, NY 10003, USA. E-mail: pance.naumov@nyu.edu
First published on 10th March 2026
Aqueous proton batteries (APBs) have attracted increasing attention due to their high safety and fast kinetics derived from proton insertion chemistry. Redox-active organic compounds with tunable molecular structures are promising electrode materials for APBs, yet their implementation remains limited by suboptimal cycling stability owing to dissolution of intermediates and structural degradation. Here, we introduce an unexplored concept of exploiting geometric isomerism to tailor structure–property relationships in proton storage. Using highly π-conjugated cis- and trans-perinone isomers as the model system, we reveal that the spatial arrangement of carbonyl groups significantly influences proton insertion processes, overpotential, kinetics, and stability. Proton insertion is confirmed to proceed via a staging mechanism involving both bare and hydrated protons, i.e. H+ and H+·(H2O)n, respectively. Co-insertion of water molecules is comparatively more into the trans-isomer relative to the cis-isomer, resulting in irreversible lattice distortion and inferior potential stability of the former. The cis-isomer delivers exceptional cycling performance with 100% capacity retention after 8000 cycles, exceeding state-of-the-art organic proton-storage materials. This work highlights the crucial role of molecular geometry in dictating electrochemical behavior and offers valuable insights for developing high-performance, practically useful organic rechargeable batteries.
Broader contextAqueous proton batteries (APBs) have emerged as a promising candidate due to the minimal ionic radius and mass and the well-established Grotthuss mechanism of the H+ mobility. Redox-active organic compounds as efficient proton hosts have recently attracted much attention due to their structure flexibility. However, their implementation remains limited by suboptimal cycling stability owing to dissolution of intermediates and structural degradation, although various physical coating or chemical modification strategies have been proposed. There is an urgent need to explore new paradigms for structural regulation from the molecular design level, aiming to fundamentally improve the intrinsic stability. Here, we introduce an unexplored concept of exploiting geometric isomerism to tailor structure–property relationships in proton storage. Using highly π-conjugated cis- and trans-perinone isomers as the model system, it is confirmed that spatial arrangement of carbonyl groups significantly influences proton insertion processes, overpotential, kinetics, and stability. The highly stable organic small molecule cis-perinone stands out with no capacity attenuation during long-term cycling. This work highlights the crucial role of molecular geometry in dictating electrochemical behavior and offers valuable insights for developing high-performance, practically useful organic rechargeable batteries. |
Polycyclic perinone, which is a commonly used industrial pigment, consists of a linearly fused naphthalene and imidazole rings (Fig. 1a).31–34 The extended π-delocalization and expanded planar arylene structure endow perinone with both chemical and thermal stability, but also account for the narrow highest occupied molecular orbital (HOMO)–lowest unoccupied molecular orbital (LUMO) gap, which facilitates charge transport.35 In our previous work, perinone-based poly(benzimidazobenzophenanthroline) (BBL)-type covalent organic frameworks (COFs) were found to have high charge mobilities and fast proton storage capability.36,37 BBL and its phenazine analogue are proven to be stable electrode materials in APBs;38,39 however, the cis-/trans-mixed perinone structure in the reported BBL or COFs poses challenges with exerting control over the spatial carbonyl group arrangement, and the behavior of these isomeric units during the redox reactions on the overall electrochemical performance remains unclear. Small intramolecular changes are expected to affect the electrochemical behavior,40,41 and in that context organic isomers are becoming increasingly relevant.42–45
Herein, we take perinone as the first example of trans/cis-isomerism for proton storage and we report a disparate effect of the disposition of its carbonyl groups on its electrochemical behavior. trans-Isomer 2 suffers from deteriorating potential upon cycling and a large overpotential, while the cis-isomer 1 remains stable both in structure and electrochemical performance. The isomerism effect governs the evolution of distinct microscopic crystal structures and affects the ion intercalation volume, ultimately leading to divergent macroscopic performance. In situ electrochemical quartz crystal microbalance (EQCM), ex situ X-ray diffraction (XRD) and crystallography reveal that the large hydrated proton species (H+·(H2O)n) inserted in 2 leads to irreversible spacing and lattice expansion, causing deteriorating performance. Impressively, 1 demonstrates a capacity retention of 100% after 8000 cycles, which is superior over state-of-the-art organic proton hosts. When coupled with the MnO2 or activated carbon (AC) cathode, both cells exhibited high cycling stability without capacity fading.
:
57. Following a reflux in an alkaline solution, a brilliant orange solid of the trans-isomer (2) precipitated first during cooling, while the blueish-red cis-perinone (1) was obtained from the filtrate via acid treatment. The axisymmetric structure of 1 gives rise to a doublet 1H peak in the nuclear magnetic resonance (NMR) spectra at 9.3 ppm, while 2 has a doublet of doublet proton signal (Fig. 1b and Fig. S1–S3). The powder XRD patterns indicate different crystal structures of the two isomers (Fig. 1c). Modelling confirmed that 1 exhibits C2v molecular symmetry and is arranged in a ‘carbonyl face-to-face’ style along the b axis, while 2 presents a Ci-symmetry in a ‘hunter's fence’ mode (Fig. 1d, e and Table S1).46–48 Solids of the two isomers were also charactered by Raman analysis and scanning electron microscopy (SEM) (Fig. S4 and S5) and both show a cudgel-shaped morphology, with 2 having a higher aspect ratio. The highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) were further evaluated by density functional theory (DFT) calculations to compare the electronic structures of two isomers (Fig. 1f). The narrower HOMO–LUMO gap (2.56 eV) of 1 indicates its superior charge transport efficiency compared to 2 (2.78 eV), 1,4,5,8-naphthaldiimide (NTCDI, 3.52 eV), PTO (3.34 eV), and diquinoxalino [2,3-a:2′,3′-c] phenazine (HATN, 3.60 eV), among other compounds (Table S2).
In a potential window of −0.3 V–0.3 V, galvanostatic measurements of 1 and 2 were carried out at 0.1 A g−1 within the three-electrode cells. The working electrodes were prepared by mixing perinone, carbon black and polytetrafluoroethylene in an optimized mass ratio of 70
:
20
:
10 on a Ti foil (Fig. S8). In the first cycle, 1 showed a discharge plateau at −0.05 V and a charge plateau at 0.05 V (Fig. 2a), affording an overpotential of ∼0.1 V. 2 exhibited a larger overpotential of ∼0.18 V, with a relatively flat discharge plateau at −0.13 V and a charge plateau at 0.05 V (Fig. 2c). The plateau indicates phase transitions of 1 and 2. There exist two pairs of sharp redox peaks for 1 and only one pair of peaks for 2 on the first-cycle dQ/dV curve (Fig. 2b and d), agreeing well with CV results (Fig. 1g) and corresponding to a two-step reaction and a one-step reaction, respectively.
Both 1 and 2 exhibit lower discharge potentials and higher overpotentials in the first cycle than others (Fig. S9) that is related to charge carriers and will be discussed in the following section. After the initial several cycles, 1 showed excellent cycling stability with nearly overlapped charge–discharge profiles within 50 cycles (Fig. 2a and b). In contrast to 1, the charge–discharge curves of 2 varied significantly, along with transition from one charge–discharge platform to three platforms (Fig. 2c and d). Such potential shift of 2 was not limited to sulfuric acid electrolyte and was also observed with 1 M HCl and 1 M H3PO4 (Fig. S10). Moreover, 1 exhibited much lower overpotentials than 2 during the whole cycling process (Fig. 2e and S11). These results indicate that cis–trans isomerism affects the electrochemical behavior of the perinone electrode. The performance of the nascent mixture of isomers as a control sample was inferior to that of 1 in stability even though the potential shift was less than that of 2 (Fig. S12).
Both 1 and 2 exhibited a high reversible capacity of 120 mAh g−1 and over 99.3% coulombic efficiency (CE) at a low current density of 0.1 A g−1, corresponding to 93% active site utilization based on its theoretical capacity of two-electron transfer (130 mAh g−1). Note that some proton hosts suffer from inevitable dissolution; stable cycling was not achievable at low rate (Fig. S13).10 Instead, high current measurements were always adopted to shorten the whole test period. Remarkably, the capacity retention of 1 and 2 was 100% after 100 cycles at 0.1 A g−1 (Fig. S14), ruling out any materials loss due to dissolution (Fig. S15). The potential shift of 2 did not influence the utilization of active sites and thus its redox activity.
In terms of rate performance, 1 delivered a high specific capacity of 90 mAh g−1 at 10 A g−1, maintaining 75% capacity at 0.1 A g−1, while 2 retained a capacity of 70 mAh g−1 at the same rate (Fig. S16 and 17). More importantly, there was only little polarization enhancement for 1 with the augment of current density from 0.1 to 10 A g−1, in sharp contrast to the double overpotential (from 120 mV to 225 mV) observed for 2. The electrochemical kinetics of the perinone isomers were further studied by CV at different scan rates. The calculated b values of dominant peaks are 0.69 and 0.66 for 1 and 0.58 and 0.55 for 2, implying slightly better kinetics of 1 over 2 (Fig. S18). Encouraged by the high electrochemical performance of the cis-isomer 1, galvanostatic charge/discharge (GCD) tests were further carried out under high mass loading. Increasing the areal loading from 1.7 mg cm−2 to 27.6 mg cm−2, 1 electrode could output a stable areal capacity of 3.4 mAh cm−2 at 0.1 A g−1 (Fig. 2f), corresponding to a specific capacity of 123 mAh g−1. Surprisingly, an increase in the areal mass loading of 1 by a factor of over ten times only increased the overpotential from 60 to 91 mV (Fig. 2g). Even at high mass loading, 1 delivered a high capacity of 68 mAh g−1 at 1 A g−1, highlighting the kinetic advantage (Fig. S19). All the results indicated excellent scale up capability of 1 electrode for practical implementation.
The long-term cycling stability of 1 electrode was analysed at 10 A g−1. As shown in Fig. 2h, after 8000 cycles, the capacity retention of 1 electrode was as high as 100%. The Raman, FTIR and SEM-EDX analyses suggested that the chemical structure and morphology of 1 are well preserved after long-term cycling (Fig. S20–22). Notably, the cycling stability of 1 surpassed state-of-the-art organic electrodes for proton storage (Fig. 2h and Table S3).27,49–56 This favorable electrochemical stability of 1 can be attributed to the high planarity, highly π-conjugated structures and strong π–π intermolecular interactions in its structure, which were reflected in the density functional theory (DFT) models (Fig. S23).57,58 The cis-isomer of 1 has a low molecular planarity parameter (MPP) value of 0.001 and balanced electrostatic potential (ESP), which are expected to reduce the molecular dissolution. Furthermore, the π electrons of 1 are well delocalized across the molecule and strong π–π interactions exist between the layers, both of which are stabilizing factors during the protonation/deprotonation process. The model for compound 2 (Fig. S24) is quite similar to that for 1, which explains the absence of capacity attenuation of 2 and active material loss during cycling, despite the unfavorable potential shift.
At mass equilibrium of the electrode, the mass and charge profiles were plotted together to explore the essence of mixed electron-ion-water transfer in detail. The slope of each region (Δm/Q) in profiles reflects mass transmission during the redox process. Results suggest that the profiles for 1 can be divided into two mass transfer stages, corresponding to the Δm/Q values of 0.017 mg C−1, 0.284 mg C−1 for discharge, and −0.228 mg C−1, −0.447 mg C−1 for the charge process, respectively (Fig. 3c). Three stages were observed for 2, in which experimental Δm/Q values with protons as the dopant were 0.02 mg C−1, 0.87 mg C−1 and 0.73 mg C−1 for discharge and −0.01 mg C−1, −0.46 mg C−1, −2.47 mg C−1 for the charge process (Fig. 3d). For the theoretical Δm/Q (0.01 mg C−1) of bare H+, the corresponding charge carriers determined for 1 are H+–H+·(H2O)1.5 during discharge and H+·(H2O)1.2–H+·(H2O)2.4 during charge. For 2, they were H+–H+·(H2O)4.6–H+·(H2O)3.9 and H+–H+·(H2O)2.5–H+·(H2O)13. For better understanding, the apparent molecular weight (Mw) of the transferred species was also calculated from Faraday's law to quantify the average number of transferred water molecules (Fig. 3e and f). Overall, the transferred water molecules (H+·(H2O)n) were 1.24/1.32 and 3.07/3.1 during charge–discharge for 1 and 2 electrodes, indicating that more water molecules were co-inserted into 2 than 1.
Subsequently, in situ FTIR spectroscopy was used to track the structure evolution of 1 during the electrochemical reactions (Fig. 4c and S30; the colorless region represents a decrease in intensity and the purple region represents an increase of band intensity). During discharging, the band from the symmetric C
O stretching at about 1700 cm−1 gradually disappeared and a new broad band from the –OH stretching emerged at 3390 cm−1, indicating the participation of the C
O groups in the redox reaction (Fig. S31). The bands from the C
C stretching of the pyridine rings at 1550 cm−1 and the C
C benzene ring stretching band at 1650 cm−1 were both enhanced after the protonation due to the formation of the quinone-vinyl moiety. Once charging, the peak of C
O groups reemerged along with weakened C
C peaks. In situ Raman and ex situ X-ray photoelectron spectroscopy (XPS) analyses were further carried out to check whether the C
N bonds contributed to the redox reactions. As shown in Fig. 4c, reduction of the band due to the C
O (1700 cm−1) stretching and enhancement of the band due to the C–O stretching (1180 cm−1) are apparent in the Raman spectra (Fig. S32 and S33). The intensity of the strong band at 1595 cm−1 representing the C–C–C stretching of naphthalene63 was reduced, while the C
C (1609 cm−1) stretching was enhanced upon discharging, in line with the FTIR spectra. The characteristic benzimidazole modes γ(C–C–C) and γ(C–C–N)64 at 1380–1400 cm−1 underwent a red-shift, accompanied by an increase in intensity, suggesting that the integrity of the benzimidazole structure was retained during the redox processes. DFT calculations confirmed that the peak shift can be attributed to changes in the conjugation of the benzimidazole fragment in the protonated structures (Fig. S34). In the XPS N 1s spectra, the C
N (398.2 eV) and C–N (400.5 eV) peaks remained unchanged in the pristine, fully discharged and charged states (Fig. 4d), a result which ruled out contribution from the C
N groups. Moreover, the DFT calculations also confirm that C
N sites exhibit much lower reduction potential based on decreased Gibbs free energy (Fig. S35). Consequently, we infer that the C
O groups, rather than the C
N groups, are the redox active site for 1 and two protonated perinone structures are entailed through a two-step one-electron transfer pathway (Fig. 4f).
), (100), (
02), (102), (200) at the respective angles of 5.9°, 7.2°, 10.5°, 11.9°, 12.7° for 1 disappeared, and two new peaks at 5.2° and 9.8° emerged (Fig. 5a).46 For 2, the (100), (002), and (102) peaks at 7.8, 11.0 and 12.9° evolved into three peaks at 5.9°, 11.4°, and 12.3° (Fig. 5b).46 These results indicate significant changes in the crystal structures of 1 and 2 after insertion of hydrated protons, a process that was also accompanied by a visible color change of the electrodes. Upon charging, the structures of both isomers recovered to their initial states.
Based on the XRD patterns of fully discharged electrodes, crystallographic data refinement calculations were then conducted to elucidate the possible protonated structures of 1 and 2. For simplicity, the calculations were performed by associating one C
O with one H3O+ for 1 and one C
O with one H+·(H2O)3 for 2. The simulated models reproduced the experimental results well and returned similar intensities (Fig. S36 and Table S4). Correspondingly, we attribute the XRD peaks at 5.2° and 9.8° of 1 to the (01
) and (02
) planes and the peaks at 5.9°, 11.4°, and 12.3° of 2 to the (100), (002) and (200) planes. Fig. 5c depicts the fully relaxed structure of protonated 1, where a unique hydrogen-bonding network connecting hydroxy–water–hydroxyl groups is formed. The strong hydrogen-bonding interactions make the protonated 1 more compact, which is evident by the shortened O-to-O distance from 7.01 Å to 6.61 Å. The ‘carbonyl face-to face’ stacking mode of 1 facilitates the embedding of the water molecule only into the carbonyl-rich rather than in the carbonyl-free space (Fig. 5d). By contrast, the heteromeral carbonyl structure of 2 made water molecules insert into the perinone layers, where only one hydrogen bond exists between each hydroxyl group and water molecule (Fig. 5e). Moreover, free H2O molecules expand the d spacing of the (100) planes from 11.3 Å to 15.1 Å in protonated 2, resulting in unit cell volume expansion from 861 Å3 to 1182 Å3 (Fig. 5f).
XRD patterns of electrodes from 1 and 2 were collected after 50 cycles. There is no obvious peak displacement for 1, suggesting that the unit cell remains practically unaffected after cycling (Fig. S37a). In contrast to 1, the (100) peak of 2 undergoes a small shift relative to that of the pristine one (Fig. S37b). When combined with the results from theoretical calculations, the irreversible displacement of the (100) peak of 2 can be attributed to the expansion of the d spacing and therefore lattice expansion caused by the insertion of the large hydrated protons, H+·(H2O)3. SEM images revealed an elongated cudgel morphology of 2 after 50 and 5000 cycles, which stands in sharp contrast to the unchanged morphology of 1 (Fig. S37c and d). We conclude that an irreversible lattice expansion caused by the larger carrier insertion is the main cause for the less stable charge–discharge profiles and the potential shift of 2. The results also suggested that the redox reaction would be affected by the local environment of the active sites, including the lattice and the molecular packing.
At a low current density of 0.1 A g−1, the MnO2@CC//1 battery delivered a high reversible capacity of 120 mAh g−1 based on the mass of 1. Impressively, the full cell exhibited high rate with slight polarization enhancement (Fig. 6d and e). Even at large current densities of 4 and 10 A g−1, the MnO2@CC//1 battery still maintained a capacity of 107 and 100 mAh g−1, respectively. The capacity completely recovered once the current was reverted to 0.1 A g−1. Additionally, a long-term cycling test was conducted at 2 A g−1. The MnO2@CC//1 battery performed quite well with a high capacity of 110 mAh g−1 after 1500 cycles (Fig. 6f and S40). It is worth noting that capacity retention approached 99%, indicating excellent long-term cycling stability. Similarly, excellent cycling performance was also achieved in the AC//1 cell except for a low working voltage of 0.5 V, where no capacity fading was observed after 4000 cycles (Fig. S41).
O groups are the active site for reversible proton storage in 1. Moreover, two types of full devices based on 1 anode and MnO2 or AC cathode also delivered high rate and cycling stability. Our work provides a fresh understanding on the isomerism effect in the field of energy storage and advances the development of high-performance organic batteries.
| This journal is © The Royal Society of Chemistry 2026 |