Proton Conduction in a Phosphonate-Based Metal–Organic Framework Mediated by Intrinsic “Free Diffusion inside a Sphere”

Understanding the molecular mechanism of proton conduction is crucial for the design of new materials with improved conductivity. Quasi-elastic neutron scattering (QENS) has been used to probe the mechanism of proton diffusion within a new phosphonate-based metal–organic framework (MOF) material, MFM-500(Ni). QENS suggests that the proton conductivity (4.5 × 10–4 S/cm at 98% relative humidity and 25 °C) of MFM-500(Ni) is mediated by intrinsic “free diffusion inside a sphere”, representing the first example of such a mechanism observed in MOFs.


* S Supporting Information
ABSTRACT: Understanding the molecular mechanism of proton conduction is crucial for the design of new materials with improved conductivity. Quasi-elastic neutron scattering (QENS) has been used to probe the mechanism of proton diffusion within a new phosphonatebased metal−organic framework (MOF) material, MFM-500(Ni). QENS suggests that the proton conductivity (4.5 × 10 −4 S/cm at 98% relative humidity and 25°C) of MFM-500(Ni) is mediated by intrinsic "free diffusion inside a sphere", representing the first example of such a mechanism observed in MOFs. F uel cells represent an appealing option as alternative clean energy systems, 1 and technologies based upon polymer electrolyte membrane fuel cells (PEMFCs) are used widely in portable applications. The design and synthesis of new protonconducting materials are of fundamental importance for the development of PEMFCs; 2 currently, commercially used proton conductors are based upon acidic polymers such as Nafion which exhibit high conductivity of 10 −2 S/cm in the presence of water. 3 However, the amorphous nature of such polymers precludes investigation of the mechanisms and/or pathways for their proton conduction, and thus analysis and feedback in order to improve future materials development are difficult to obtain. 4 Constructed from metal ions and organic linkers, metal−organic framework (MOF) materials often display high surface areas, high porosity, and, more importantly, extended crystalline structures, and significant focus has been placed on their applications in gas storage, separation, and catalysis. 5−7 More recently, MOFs have appeared as promising new candidates as porous materials for proton conduction. 8,9 Functionalization of the organic linker in MOF materials allows the periodic introduction of acidic groups (e.g., SO 3 H, PO 3 H 2 ), 10−12 and their intrinsic porosity enables the loading of additional protonic molecules (e.g., imidazole, histamine) within the pore, 13−15 to yield decorated materials with improved proton conductivity. Most importantly, the crystalline nature of MOFs provides an excellent platform to interrogate possible proton hopping and conduction pathways, thus enabling the construction of structure−activity relationships, which cannot be achieved in polymer-based systems due to their intrinsic lack of long-range order. 1 X-ray crystallographic studies afford average positions of protons within the extended lattice, mostly on O-atoms from hydroxy groups or water molecules. These are often refined by a "riding model" and subject to large uncertainties inherent in these X-ray experiments. Understanding the dynamics of these active protons is of fundamental importance for the design of improved materials. However, gaining such information within porous lattices of MOFs is very challenging, and knowledge on the dynamics of proton diffusion in MOFs is rarely reported. 16 Pulsed field gradient (PFG) NMR is capable of measuring the dynamics of specific species in materials containing intrinsic structural defects. However, PFG-NMR is likely to underestimate the diffusion rate, as it operates on longer time and spatial scales. 17 Here we describe the synthesis and crystal structures of two novel isostructural phosphonate-based MOFs, [M 3 (H 3 L) 2 (H 2 O) 9 (C 2 H 6 SO) 3 ] (M = Ni, Co; H 6 L = benzene-1,3,5-p-phenylphosphonic acid), denoted as MFM-500(Ni) and MFM-500(Co), respectively (MFM = Manchester Framework Material). MFM-500(Ni) and MFM-500(Co) adopt 2D networks, in which the ligand is only 50% deprotonated (H 6 L → H 3 L 3− ) and bound to the metal ions. This affords free acidic protons from the partially protonated ligand within the lattices of these coordination complexes. The materials show proton conductivities of 4.5 × 10 −4 [MFM-500(Ni)] and 4.4 × 10 −5 S/cm [MFM-500(Co)] at 98% relative humidity (RH) and 25°C. Single-crystal X-ray structural analyses reveal a potential proton hopping pathway constructed from the free phosphonic acid groups and coordinated water molecules on the metal centers, subject to uncertainties as discussed above.
More importantly, quasi-elastic neutron scattering (QENS) has been used to investigate the intrinsic mechanism of proton diffusion within MFM-500(Ni), and this study suggests that the proton conduction in MFM-500(Ni) is mediated by the model of "free diffusion inside a sphere" rather than "jump diffusion between sites".
Benzene-1,3,5-p-phenylphosphonic acid, H 6 L, was synthesized following a literature procedure to afford a white powder in 67% yield. 18 Upon mixing M(NO 3 ) 2 (M = Ni, Co) with H 6 L in a 2:1 molar ratio in a solution of H 2 O/DMSO/DMF, hexagonal column-shaped single crystals of MFM-500(Ni) and MFM-500(Co) grew within 2 days at 40°C. Synchrotron single-crystal X-ray diffraction (XRD) reveals that the two materials are isostructural, and both crystallize in the hexagonal space group P6 3 /m (Table S1) with a 2D layered structure. The tris-phosphonate ligand, H 3 L 3− , acts as one type of threeconnected node in overlapping pairs of 2D hexagonal (6,3) networks lying in the ab plane, in which the other type of threeconnected node is a disordered mixture of Figure  1a). Running down the c-axis are stacks of pairs of overlapping ligands in which each pair is rotated by 60°with respect to the pairs above and below (ligand separations along the c-axis are 3.85 Å between overlapping pairs). Columns of disordered metal centers running along the c-axis all have pseudooctahedral coordination geometries made up of O donors from deprotonated phosphonate hydroxyl groups, neutral water molecules, and partially resolved DMSO molecules. These O donors act as either monodentate or μ 2 -bridging ligands, depending on the occupancy of the disordered adjacent metal cation sites. The columns of M 1 and M 2 nodes along the c-axis reside in the hexagonal holes within the pairs of offset hexagonal networks above and below them ( Figure S1).
In MFM-500(Co), residual electron density peaks were observed for H-atoms on the water molecules and one of the phosphonate hydroxy groups; this has allowed their positions to be refined using suitable geometric restraints, revealing an extensive H-bonding network within and between the columns of metal cations (Figure 1b; Table 1). The H-bond donor− acceptor (D−A) distances range from 2.45 to 3.01 Å and form crossover chains that run between adjacent columns of metal cations via reciprocally disordered H-bonds between symmetryrelated phosphonate−hydroxy oxygen atoms O33 (D−A = 2.45 Å). O33 is also involved in disordered H-bonding with a bridging water molecule O2W (D−A = 2.73 Å), while the bridging water molecule O1W donates a bifurcated H-bond to two symmetrically equivalent unbound deprotonated phosphonate−hydroxy oxygen atoms O34 (D−A = 3.01 Å). The passage of H-bond chains along the c-axis is interrupted by the presence of a partially occupied disordered DMSO solvent molecule; it is likely that the chains continue through this region conveyed by an unmodeled disorder solvent component.  (Table 1).
Thermal gravimetric analysis (TGA) of MFM-500(Ni) and MFM-500(Co) shows similar weight loss steps, with a slightly lower stability observed for MFM-500(Co) ( Figure S3). In situ variable-temperature PXRD data confirmed a reversible phase transition occurring between 75 and 100°C for both materials (Tables S3 and S4). This result is consistent with the first step in the TGA plots and is associated with a color change from green to brown for MFM-500(Ni) and from pink to purple for MFM-500(Co). The original phase can be recovered by exposing the dehydrated samples to air (or water vapor), suggesting a reversible change in the coordinated water molecules upon dehydration/rehydration. This is further evidenced by the dehydrated samples remaining intact and stable under a flow of dry O 2 . Due to the quality of the data and the complexity of these structures, attempts to determine the crystal structures of the dehydrated phases have thus far been unsuccessful. Samples of desolvated MFM-500(M) (M = Ni, Co) showed very low BET surface areas (<10 m 2 /g), as determined from N 2 adsorption isotherms at 77 K.

Journal of the American Chemical Society
Communication plots contain an incomplete semicircle in the high frequency region and a pronounced tail at low frequencies, consistent with blocking of protons at the electrodes (Figure 2). At room temperature and 98% RH, the proton conductivities for MFM-500(Ni) and MFM-500(Co) were measured as 4.5 × 10 −4 and 4.4 × 10 −5 S/cm, respectively. The difference in the conductivities between these isostructural materials probably correlates to the bond strength between coordinated water molecules and the metal cations, as suggested by the TGA plots and solid-state UV/vis absorption spectra (see SI respectively. 20 Impedance data were also measured at 98%, 75%, 45%, and 0% RH at 25°C for MFM-500(Ni) and MFM-500(Co); both materials showed a steady decrease of conductivity with decreasing %RH (Figures S6 and S7; Table  S5). At 0% RH both materials showed no apparent conductivity (<10 −9 S/cm). This is typically observed for water-mediated proton conductors. It is worth noting that the dehydrated materials show no apparent proton conductivity, and their conductivities (4.5 × 10 −4 and 4.4 × 10 −5 S/cm for the Ni and Co compounds, respectively) can be recovered by exposing the dehydrated samples to water vapor, consistent with the reversible phase change shown by the in situ PXRD data. The activation energy (E a ) of the proton conduction in MFM-500(Ni) was estimated from the impedance spectra recorded at 98% RH between 15 and 32°C to be 0.43 eV ( Figure S8). Two main mechanisms for proton diffusion are the Vehicle mechanism (typically E a > 0.4 eV) and the Grotthuss mechanism (typically E a < 0.4 eV). 8 For MFM-500(Ni) E a lies at the boundary of the two mechanisms, indicating that it is likely that proton conduction in MFM-500(Ni) is governed by an intermediate process between the Grotthuss and Vehicle mechanisms. This behavior has been reported in a few cases. 2,21,22 For example, (NH 4 ) 2 (adp)[Zn 2 (ox) 3 ]·3H 2 O (adp = adipic acid; ox 2− = oxalate) showed a high proton conductivity of 8 × 10 −3 S/cm at 25°C under 98% RH with E a = 0.63 eV, and the mechanism of proton conduction was assigned as mixed Grotthuss and Vehicle types. 21,22 We sought to gain further understanding of the mechanism of proton conduction in MFM-500(Ni) using QENS. Data for MFM-500(Ni) were collected between −23 and 150°C under both anhydrous and 98% RH conditions to study the dynamics of protons through the framework lattice. The elastic incoherent structure factors (EISF) were extracted from the QENS spectra to gain the geometrical information on the molecular motions of active protons in MFM-500(Ni) (see SI). The EISF plots showed clear Q-dependence and were carefully fitted with all well-known theoretical models for proton diffusion (i.e., jumping between n sites and various free diffusion models as shown in Figure S10). It has been found that the proton diffusion in MFM-500(Ni) is best described by the model of "free diffusion inside a sphere" (eq 1) rather than the model of jumping between n sites (Figures 3a and  S10). 23 where j 1 is the first-order spherical Bessel function, r is the radius of the sphere, and p and (1 − p) are the immobile and mobile fractions of the protons involved in this process, respectively. The best fitting to EISF plot was observed for r = 2.25 Å ( Figure S11), entirely consistent with the observed H··· A distances from the single-crystal X-ray structure (taking an O−H bond distance of ∼0.84 Å). The dependency of the spherical free diffusion on the distance r is confirmed by the extraction of the half-width of the half-maximum (HWHM, Γ) of the QENS spectra as a function of Q 2 at different temperatures ( Figure 3b) Figure S12).
In conclusion, we have combined the single-crystal XRD and QENS spectroscopy to study the molecular mechanism for proton conduction in a new phosphonate-based MOF material. These complementary static and dynamic approaches yield highly consistent results and afford direct visualization of the pathway and mechanisms of proton transfer in the framework

Journal of the American Chemical Society
Communication lattice. For the first time, the model of "free diffusion inside a sphere" has been experimentally confirmed in protonconducting MOFs. Further efforts to optimize the hydrogenbonding network in these materials via ligand modification and post-synthetic approaches are currently underway.

* S Supporting Information
The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/jacs.6b02194. CCDC-1450010 and CCDC-1450011 contain the supplementary crystallographic data for this paper.