Byoung Gwan
Lee
a,
Dongwook
Kim
b,
Jin Young
Bae
a,
Ji Woo
Jeong
a and
Dae-Woon
Lim
*a
aDepartment of Chemistry, Yonsei University, 1 Yonseidae-gil, Wonju, Gangwon-do 26493, Republic of Korea. E-mail: limdaewoon@yonsei.ac.kr
bCenter for Catalytic Hydrocarbon Functionalization and Institute of Basic Science, Daejeon 34141, Republic of Korea
First published on 20th August 2025
Understanding the correlation between structural variations and proton transport in coordination polymers (CPs) is essential for developing efficient solid-state proton conductors (SSPCs). In this study, we demonstrate that post-synthetic transmetalation via alkali metal exchange enhances proton conductivity in a pseudo-three-dimensional yttrium-based CP, {[Y(H2O)4(HDSNDC)]·H2O}n (H4DSNDC = 4,8-disulfonaphthalene-2,6-dicarboxylic acid). Immersion in 1 M KCl solution induces transmetalation, yielding a potassium-substituted analogue, {K3(H2O)2(HDSNDC)}n. Structural analysis reveals that transmetalation induces structural reorganization, in which terminal H2O ligands are converted into μ-bridging ones. These bridging H2O ligands exhibit increased acidity due to cooperative polarization by adjacent K+ ions, facilitating proton dissociation and significantly enhancing conductivity from 7.23 × 10−5 S cm−1 for Y-DSNDC to 2.50 × 10−3 S cm−1 for K-DSNDC at 363 K under 95% RH. This work highlights transmetalation-induced coordination tuning as an effective strategy to enhance proton transport in CPs.
Recently, coordination polymers(CPs), constructed by coordination bonds between metal ions and organic linkers, have garnered increasing attention as solid-state proton conductors owing to their structural designability, tunability, thermal stability, and superior proton transport properties under mild temperature conditions.9–16 To enhance their proton conductivity, various strategies have been explored, including the introduction of predesigned functionality in the frameworks via de novo synthesis17,18 (e.g., employing the ligands including extra acidic functional groups) and post-synthetic metal and ligand modifications, such as controlling the coordination environment of metal centres,19,20 incorporating hydrophilic acid functional groups in the ligands (e.g., –SO3H, –COOH, –PO3H2),21,22 and introducing guest molecules such as H2O, NH3, imidazole, and inorganic acids.23–27 These approaches have not only facilitated hydrogen bonding networks but also increased the mobile proton concentration.28,29 While these strategies have shown promise, they do not inherently guarantee efficient proton conduction. The effectiveness often depends on the spatial arrangement and chemical composition of the resulting conduction pathways. The decisive factor for high conductivity lies not merely in the presence of conducting media (e.g., H2O, NH3, or acids) and functional groups, but in how these components cooperate to create geometrically favourable and dynamically interactive environments for proton diffusion. Therefore, structural designs that promote the formation of conjugated acid–base systems, such as H3O+ + H2O ⇌ H2O + H3O+ or OH− + H2O ⇌ H2O + OH− between guest molecules and nearby acidic or basic sites, can be more critical for facilitating efficient proton transport.30–33 Among emerging strategies to engineer such environments, post-synthetic metal ion exchange (transmetalation) has attracted considerable attention for its ability to modulate local coordination environments while preserving crystallinity via solution-mediated transmetalation.34–36 In this approach, the coordination mode of H2O molecules, whether terminal or bridging, can often be modulated, leading to variations in their pKa. Notably, the replacement of high valence metal centres with alkali metal cations, which possess lower valency and large ionic radii, enables high metal ion densities, thereby favouring the formation of μ-bridging water ligands between adjacent centres.37–39 These bridging waters exhibit enhanced acidity due to cooperative polarization, which strengthens hydrogen-bond networks and increases their proton-donating ability, thereby increasing the concentration of mobile protons.40–42 Although structural modulation via transmetalation is a promising approach to promote efficient proton conduction, to the best of our knowledge, there are scarce reports, and the detailed mechanistic understanding of how exchanged alkali ions govern proton transport remains incomplete. Herein, we demonstrate an enhancement in proton conductivity of Y-DSNDC {[Y(H2O)4(HDSNDC)]·H2O}n (H4DSNDC = 4,8-disulfonaphthalene-2,6-dicarboxylic acid) via transmetalation with K+ under ambient conditions, resulting in {K3(H2O)2(HDSNDC)}n (Scheme 1). During the process, four terminally coordinated H2O molecules in {[Y(H2O)4(HDSNDC)]·H2O}n were changed to two bridging H2O (μ-H2O) ligands protruding into the pore direction and exhibiting an improved proton conductivity from 7.23 × 10−5 to 2.50 × 10−3 S cm−1 at 363 K.
![]() | ||
Scheme 1 (a) Schematic illustration of the transmetalation process. (b) Crystal images of Y-DSNDC (left) and K-DSNDC (right). |
Two structurally defined coordination polymers, Y-DSNDC {[Y(H2O)4(HDSNDC)]·H2O}n and its metal-exchanged analogue K-DSNDC {K3(H2O)2(HDSNDC)}n, were synthesized using H4DSNDC, a sulfonic acid-functionalized naphthalene dicarboxylic acid ligand. The dual –SO3H/–COOH functionality of the ligand can enhance its hydrophilicity, thereby promoting dense hydrogen-bonding networks essential for proton conduction. Y-DSNDC was obtained as single crystals via a hydrothermal reaction using Y(NO3)3·6H2O and H4DSNDC at 120 °C. Single-crystal X-ray structural analysis reveals that Y-DSNDC has a monoclinic crystal system with the P21/m space group (Table S1). The coordination of yttrium ions and H4DSNDC extends the structure along the bc plane, forming a 2D framework, incorporating one lattice H2O molecule. The two carboxyl groups of HDSNDC each coordinate to yttrium ions in the ac plane in a monodentate fashion, with the coordinating oxygen atoms originating from the carbonyl (CO) moieties. One carboxyl group is deprotonated (–COO−), while the other remains protonated (–COOH) and forms a hydrogen bond with the adjacent carboxylic acid and carboxylate groups (O⋯O = 2.552 Å) (Fig. S2). One of the two sulfonate functional groups remains uncoordinated, while the other bridges yttrium ions along the b-axis (Fig. 1c). Specifically, Y1 exhibits eight-coordinate geometry, coordinated by four terminal H2O molecules (O24, O25, O26, O26*), two oxygen atoms from two sulfonate groups (O16, O16*), and two oxygen atoms from the carboxylate group (O19) and carboxylic acid group (O23) of the HDSNDC ligands (Fig. 1a). The Y1–O bond lengths range from 2.237 to 2.439 Å. According to continuous shape measurement (CShM) calculations of metal coordination geometry using the SHAPE software (Table S2),43 the coordination geometry of Y1 approximates a triangular dodecahedron. When viewed along the c-axis, the yttrium centres are arranged in a zigzag pattern along the b-axis (Fig. 1c). The bc planes are interconnected through hydrogen bonding between the terminally coordinated H2O molecules and the uncoordinated sulfonate functional group along the a-axis. Overall, a two-dimensional coordination network of metal ions and ligands forms a pseudo-three-dimensional structure through interlayer hydrogen bonding (Fig. S3). Hydrogen bonding between the ligand and coordinated H2O molecules along the b-axis exhibits distances of 2.675 and 2.824 Å (OH2O⋯OSO3−) (Fig. S6b). The ligand molecules are planar and oriented parallel to the ac-plane, generating rectangular pores extending along the c-axis (6.8 × 3.4 Å, Fig. 1c). As shown in Scheme 1 and Fig. S4, immersing the rod-shaped Y-DSNDC crystals in a 1 M aqueous KCl solution for 72 hours induces time-dependent transmetalation, gradually transforming them into plate-shaped K-DSNDC. K-DSNDC crystals could not be directly synthesized under the same conditions. Single-crystal X-ray diffraction analysis reveals that K-DSNDC retains the monoclinic crystal system but adopts a different space group, C2/c (Table S1). For charge balance, one trivalent Y3+ was substituted by three monovalent K+ ions, resulting in a reorganization of the coordination environment. A metal cluster is formed with two distinct K+ coordination geometries, and terminally coordinated H2O molecules in Y-DSNDC are converted to μ-bridging H2O ligands. Furthermore, all sulfonate functional groups in the ligands coordinate to metal centres, resulting in 3D networks. Within the metal cluster, K1 adopts an eight-coordinate biaugmented trigonal prismatic geometry. It is coordinated by two μ-H2O molecules (O15, O15*) that bridge K1 with K2 and K2*, four oxygen atoms from two sulfonate ligands acting in a bidentate mode (O4, O5, O4*, O5*), and two oxygen atoms (O6*) from two monodentate sulfonate groups, as shown in Fig. 1b. In contrast, K2 adopts a seven-coordinate capped octahedral geometry, bonded to one μ-H2O (O15), one monodentate sulfonate oxygen (O5*), two oxygens from bidentate sulfonate (O4*, O6*), one μ-bridging sulfonate oxygen (O5), and two monodentate oxygens from carboxylic acids (O13, O14*). The orientation of the ligands remains similar to that of Y-DSNDC, with adjacent carboxylate groups of ligands in K-DSNDC connected along the c-axis through O–H⋯O hydrogen bonding (O⋯O = 2.470 Å) (Fig. S6). Although the ligands lie parallel to the ac-planes, they exhibit a torsion angle of 28.6° between the naphthalene rings and carboxylic acid groups due to coordination of all oxygen atoms in –SO3− functional groups, resulting in a tilted stacking configuration. This conformational change enlarges the one-dimensional rectangular pores along the c-axis, yielding channels approximately 7.3 Å wide and 3.6 Å high (Fig. 1d).
![]() | ||
Fig. 1 The coordination environment of the (a) Y3+ ion and (b) K+ ions. The cell packing along the c-axis of (c) Y-DSNDC and (d) K-DSNDC. |
The powder X-ray diffraction (PXRD) patterns of the as-synthesized Y-DSNDC and after immersion of Y-DSNDC crystals in 1 M aqueous KCl solution for 72 hours are equivalent to the simulated patterns obtained from single-crystal X-ray data, confirming the successful formation of a phase-pure crystalline material. This is consistent with recrystallization into a new phase via transmetalation, without any detectable remnants of the original Y-DSNDC phase (Fig. 2). Both coordination polymers exhibit high chemical stability in water and common organic solvents (Fig. S9). In the FT-IR spectrum of H4DSNDC, the CO stretching band appears at 1695 cm−1. Upon coordination to Y3+ and K+, this band shifts to 1683 and 1682 cm−1, respectively, confirming metal-carboxylate bond formation. In addition, the free ligand shows characteristic sulfonic acid vibrations at 1192, 1112, and 1046 cm−1, corresponding to S
O and S–O stretching modes. In Y-DSNDC, these bands are split and shifted into distinct pairs at 1195/1164, 1115/1089, and 1054/1025 cm−1. This splitting arises from the coexistence of coordinated and non-coordinated sulfonate groups within the structure (Fig. S1a), resulting in asymmetric vibrational environments in line with previous reports.44 In K-DSNDC, vibration characteristics of sulfonate groups are observed at 1191, 1116, and 1044 cm−1, confirming the presence of sulfonate groups in the framework (Fig. S11).
The structural features of Y-DSNDC and K-DSNDC are further supported by X-ray photoelectron spectroscopy (XPS) analysis (Fig. S12 and S13). In Y-DSNDC, the S 2p region exhibits two distinct doublets, indicative of asymmetric coordination of the sulfonate groups, in agreement with the IR spectral splitting and the single-crystal X-ray structure. In contrast, K-DSNDC exhibits a single S 2p doublet, reflecting a uniform coordination environment. Both compounds exhibit excellent thermal stability in thermogravimetric analysis (TGA), as evidenced by structural decomposition occurring only above 450 °C. The initial weight loss below 200 °C arises from the release of lattice and coordinated water molecules, which is fully reversible. Y-DSNDC shows a three-step weight loss of approximately 15% below 200 °C, consistent with five water molecules in the structure (calculated value: 16.3%). The first step, approximately 6% below 90 °C, corresponds to the release of two molecules, followed by another 6% up to 150 °C, and a final 3% by 200 °C. In contrast, K-DSNDC exhibits a single weight loss of 6.89% between 110 and 150 °C, consistent with the loss of two bridging water molecules (calculated: 6.8%) (Fig. S14). To confirm the porosity, N2 gas sorption measurements were carried out at 77 K. Neither sample exhibits significant N2 uptake in the low-pressure region (P/P0 < 0.2), consistent with the single-crystal X-ray structures indicating pore apertures smaller than the kinetic diameter of N2 (3.64 Å). However, K-DSNDC shows a pronounced increase in uptake near saturation (P/P0 ≈ 1.0), which can be attributed to multilayer adsorption on the external surface or condensation within interparticle voids (Fig. 3a).45 To further investigate the hydration behaviour, H2O vapor sorption measurements were carried out after thermal activation at 373 K for 24 h under vacuum. During this process, the removal of lattice and coordinated water induced structural changes, as confirmed by PXRD (Fig. S15). H2O vapor sorption isotherms recorded at 298 K revealed stepwise uptake in both materials with increasing relative humidity (RH). Y-DSNDC adsorbed up to 5.6 H2O molecules per formula unit until P/P0 = 0.95, of which approximately 5 H2O molecules remained after desorption, indicating strong retention within the framework. Similarly, K-DSNDC uptakes 2.3 H2O molecules per formula unit and retains approximately 2 H2O molecules upon desorption. In both cases, PXRD patterns collected after H2O vapor sorption revealed a full structural recovery to the original structure, suggesting reversible framework rehydration driven by water incorporation (Fig. 3b and S15). To further investigate the uptake behaviour and host–guest interactions of water molecules after structural recovery, second-cycle H2O vapor sorption measurements were conducted. Second-cycle H2O vapor sorption measurements revealed that K-DSNDC exhibited slightly higher uptake of guest water molecules than Y-DSNDC in the range of P/P0 = 0.1–0.8, which is attributed to its expanded pore dimensions as shown in structural analysis. In addition, the steep uptake observed for K-DSNDC at P/P0 ≈ 0.9 can be ascribed to multilayer adsorption in the pore. A hysteresis loop between the adsorption and desorption branches was observed, consistent with the thermodynamic irreversibility of capillary condensation, which means the presence of strong interaction sites with guest water molecules in the pores.46 This behaviour may also be influenced by the slow desorption of water molecules due to strong hydrogen bonding interactions within the framework (Fig. S16). Based on the structural analysis, this phenomenon can be attributed to the μ-bridging H2O molecules protruding into the pore direction. The structural modifications induced by metal exchange, including expanded pore dimensions and changes in metal-aqua coordination environments, provided a compelling rationale to investigate their impact on proton conductivity. For proton conductivity measurements, pelletized samples of Y-DSNDC and K-DSNDC were analyzed using electrochemical impedance spectroscopy (EIS) under variable relative humidity (RH) and temperature conditions. As shown in Fig. 4a, proton conductivity was measured at 298 K while increasing the relative humidity from 30% to 95%. At 30% RH, Y-DSNDC exhibited a conductivity of 9.19 × 10−11 S cm−1, while K-DSNDC showed a markedly higher value of 2.44 × 10−8 S cm−1. Under low-humidity conditions, the higher conductivity of K-DSNDC compared to Y-DSNDC can be attributed to the lower pKa of μ-bridging H2O molecules, despite Y-DSNDC containing more terminally coordinated H2O. As humidity increased, the proton conductivities of both materials showed strong moisture dependence, reaching 1.56 × 10−4 S cm−1 for Y-DSNDC and 2.89 × 10−4 S cm−1 for K-DSNDC, respectively, at 95% RH. The temperature-dependent proton conductivities at 95% RH are presented in Fig. 4b. Y-DSNDC exhibits a gradual increase in conductivity with rising temperature, reaching a maximum of 2.55 × 10−4 S cm−1 at 343 K. Above this point, the conductivity progressively declined, reaching 7.23 × 10−5 S cm−1 at 363 K. In contrast, K-DSNDC shows a continuous increase, attaining 2.5 × 10−3 S cm−1 at 363 K (Fig. 4c). Activation energies (Ea), derived from Arrhenius plots, were 0.13 eV for Y-DSNDC and 0.30 eV for K-DSNDC, consistent with the Grotthuss-type proton diffusion, where protons migrate via hydrogen-bond-assisted hopping (Fig. 4b). Notably, the maximum proton conductivity of K-DSNDC (2.5 × 10−3 S cm−1 at 363 K and 95% RH) surpasses the highest values reported for K-based CPs or MOFs, including the 3D frameworks [K2(OOCCH(OH)PO3H)(H2O)2] (K-HPAA) (1.3 × 10−3 S cm−1 at 297 K under 98% RH)47 and {[K8(PTC)2(H2O)1.5]·4H2O}n (1.0 × 10−3 S cm−1 at 298 K under 98% RH),48 highlighting its outstanding proton transport performance among K-based materials. In Y-DSNDC, the decrease in conductivity from 343 K to 363 K coincides with the release of two coordinated water molecules from the structure, as confirmed by TGA (Fig. S14). A comparison of the hydrogen-bonding environments around the Y3+ centre indicates that the two coordinated water molecules (O26, O26*) each form only two hydrogen bonds, the fewest among all structural water molecules, and are thus weakly bound and preferentially lost upon heating (Fig. S17). Their loss disrupts the hydrogen-bonding network, as shown in Fig. S18, leading to a decline in the proton transport pathway. Although PXRD analysis after EIS measurements at 363 K (Fig. S21) confirms that the overall crystalline framework remains intact, the removal of these water molecules significantly weakens the hydrogen-bonding network and diminishes proton transport efficiency. In contrast, K-DSNDC exhibited a continuous increase in proton conductivity with rising temperature, which reflects the thermal retention of tightly bound water molecules that preserve the hydrogen-bonding network. The superior proton conductivity of K-DSNDC across the entire humidity and temperature ranges, relative to Y-DSNDC, can be ascribed to multiple structural advantages. First, the expanded pore apertures in K-DSNDC allow for greater uptake of guest water molecules, which enhances the formation of extended hydrogen-bonding networks. Second, μ-bridging water molecules are likely to exhibit higher acidity, facilitating proton dissociation and increasing the concentration of mobile protons. Finally, the spatial proximity between acidic bridging water molecules and guest water enables efficient hydrogen-bonded interactions that facilitate proton hopping (Fig. S22). Notably, no degradation in performance was observed over several days, suggesting that the proton conductivity of K-DSNDC remains stable in the long term (Fig. 4d).
CCDC 2463061 and 2463062 contains the supplementary crystallographic data for this paper.49
The authors confirm that the data supporting this article have been included as part of the SI. Crystallographic data for Y-DSNDC (CCDC2463061) and K-DSNDC (CCDC2463062) have been deposited at the Cambridge Crystallographic Data Center (CCDC).
This journal is © The Royal Society of Chemistry 2025 |