Achieving Superprotonic Conduction in Metal–Organic Frameworks

Dec 22, 2017 - AC impedance data was collected by cycling between 106 and 1 Hz with 10–200 mV of applied potential using VersaStudio software. .... ...
0 downloads 9 Views 680KB Size
Subscriber access provided by RMIT University Library

Article

Achieving Superprotonic Conduction in Metal-Organic Frame-works through Iterative Design Advances SiRim Kim, Biplab Joarder, Jeff Hurd, Jinfeng Zhang, Karl W Dawson, Benjamin S. Gelfand, Norman E. Wong, and George K. H. Shimizu J. Am. Chem. Soc., Just Accepted Manuscript • DOI: 10.1021/jacs.7b11364 • Publication Date (Web): 22 Dec 2017 Downloaded from http://pubs.acs.org on December 22, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Journal of the American Chemical Society is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 8 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

Achieving Superprotonic Conduction in Metal-Organic Frameworks through Iterative Design Advances SiRim Kim, Biplab Joarder, Jeff A. Hurd, Jinfeng Zhang, Karl W. Dawson, Benjamin S. Gelfand, Norman E. Wong, George K. H. Shimizu* Department of Chemistry, University of Calgary, Calgary, AB, Canada, T2N 1N4 ABSTRACT: Two complementary design strategies – isomorphous ligand replacement and heterocycle doping – have been applied to iteratively enhance the proton conductivity of a metal-organic framework: β-PCMOF2. The resulting materials, PCMOF2½(Pz) and PCMOF2½(Tz) (Pz= 1H-pyrazole, Tz= 1H-1,2,4-triazole), have their proton conduction raised almost two orders of magnitude compared to β-PCMOF2. The bulk conductivities of these materials are over 10-1 Scm-1 at 85 °C and 90% relative humidity (RH), while maintaining the parent MOF structure. A solid state synthetic route for doping 1-D channels is also presented.

INTRODUCTION Metal−organic frameworks (MOFs) have received growing 1-14 interest as proton conducting materials in recent years. This is mainly because of the modular nature of MOF design and synthesis as well as their crystallinity allowing structure 15-21 activity relationships to be developed. Many protonconducting MOFs have focused on incorporating proton transfer agents within the pores, both to improve conductivity and the upper operating temperature, or functionalizing co-ordinatively unsaturated metal sites or tuning the acidity 22-36 of the pores by incorporating specific functional groups. The following Arrhenius equation, which is derived from the Nernst-Einstein relation, provides a fundamental guideline in 37-45 designing high proton conducting MOFs.



  

    

exp

 

→

 

exp

 

 ∙∙∙∙∙∙

(1)

ers can be achieved by introducing acidic moieties such as H2SO4, RSO3H, H3PO4, water/hydronium, RCO2H, or protonated N-heterocycles into the framework or as guests in the 46-56 proton conduction channels. Only about 50% of the charge carriers should be protonated as proton hopping requires non-protonated sites to act as proton acceptors. Finally, a proton carrier will need to reorient itself following a proton transfer and the ease with which this occurs is determined by how the carriers interact with the pore and each other; this is intimately related to the structure of the MOF. Therefore, a desirable proton conducting metal organic framework would have a 3-dimensional pore structure with large number of available protons and proton carriers that have high rotational freedom. Proton carriers would interact weakly with the framework and they should be at a distance to allow for facile proton hopping but not be so close as to hinder reorientation following proton transfer.

σ = Ionic conductivity n = Number of charge carriers (temperature independent mobile site occupancy) e = Charge of the mobile ion Do = Constant which is related to the mechanism of ionic conductivity k = Boltzmann constant T = Temperature ΔSm = Motional entropy Ea = Motional enthalpy (activation energy for ion transport) According to the above equation, better and high conducting MOFs can be achieved via three different strategies: 1) by increasing the number of charge carriers, 2) by decreasing the activation energy of proton transfer, and 3) by increasing the motional entropy. Increasing the number of charge carri-

Figure 1. (Left) The space-filling cross section structure of βPCMOF2 with its 1-dimensional proton conduction channel illustrated with an arrow and the schematic illustration of isomorphous replacement. (Right) A 2-dimensional crystal layout of β-PCMOF2 where the pores are impregnated with heterocycles (teal pentagons). A proton conducting MOF named β-PCMOF2 (Figure 1 left) was reported in our previous study and the conductivity val-3 -1 ue was found to be 1.3 × 10 Scm at 90% relative humidity 57 and 85 °C. β-PCMOF2 consists of a trisodium 2,4,6trihydroxy-1,3,5-trisulfonate benzene(Na3L1) complex ar-

ACS Paragon Plus Environment

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ranged in a 1D columnar structure where the pores are lined 57 with sulfonate oxygen atoms. The pores are 5.6 Å in diameter and are hydrophilic in nature. Further, the proton conductivity value of β-PCMOF2 was improved 1.5 orders of magnitude by isomorphous ligand replacement with a more protic ligand of comparable size and shape (1,3,5benzenetriphosphonic acid, Na3H3L2) and the resultant -2 -1 MOF named β-PCMOF2½ showed 2.1 × 10 Scm at 90% 57 relative humidity and 85 °C. This enhancement in conductivity is due to increase the number of freely available acidic 57,58 protons from the backbone. In another study, we successfully enhanced the anhydrous proton conductivity of βPCMOF2 by loading the pores with less volatile, amphiprotic guest molecules (Figure 1 right). The report describes a βPCMOF2(1,2,4-Triazole) system where β-PCMOF2 was loaded with 0.3 equivalent of triazoles in the pores and showed a five orders higher anhydrous conductivity of magnitude -4 -1 relative to unloaded β-PCMOF2 to reach 5 × 10 Scm at 150 49 °C. Thus combination of nonvolatile heterocyclic compound loading and isomorphous ligand replacement in βPCMOF2 would be an excellent strategy to further enhance both the water assisted and anhydrous proton conductivity for MOF materials.

Scheme 1. Chemical structure of L1, H3L2, 1H-pyrazole and 1H-1,2,4-triazole. Herein, we report a study assessing the effect of heterocycle loading, isomorphous protic ligand replacement, and the effect of merging these approaches on proton conduction in a single family of isostructural MOFs. Three new compounds, PCMOF2½(Pz), PCMOF2½(Tz) and PCMOF2(Pz) are reported and compared under common conditions to PCMOF2(Tz) and PCMOF2½. The ligand H3L2 had previously been re59 ported in a structure named PCMOF3 and so the mixed ligand systems were named PCMOF2½. The systematic variation of components provides insights to the extent different factors can affect total proton conduction and activation energy. To obtain good heterocycle loadings in these systems, solid state syntheses are employed, a method that is potentially applicable to other MOF syntheses where competition from solvent is problematic. Ultimately, the application of the two design strategies results in the synergistic enhancement of the β-PCMOF2’s proton conductivity to -1 -1 yield MOFs that conduct protons over 10 Scm . EXPERIMENTAL PROCEDURE

General synthesis of mixed ligand-heterocycle loaded MOFs. As-synthesized PCMOF2(Heterocycle) and Na3H3L2 were individually ground into fine powders using a mortar and pestle and mixed in a 2:1 ratio with respect to Na3L1 and Na3H3L2. The ground powders were placed in a vial and mechanically shaken. The resulting mixture was placed into a silicon tube, sealed and pressed under a hydrostatic pressure 2 of 10,000 pounds/inch for 2 minutes. The resulting pellet was removed from the silicone tube and placed into a 23 mL Teflon autoclave, along with a vial of water (1.7 mL). The Teflon autoclave was then sealed in a stainless steel jacket and heated to 80 °C for 48 hours then cooled back to room temperature over 12 hours. The pellet was then placed in a desiccator for minimum of 12 hours to remove excess moisture, and ground into a fine powder using a mortar and pestle. The resultant powder was analyzed using powder X-ray diffraction (PXRD), thermogravimetric analysis (TGA), ele1 31 mental analysis, and H and P NMR to verify the completion of the isomorphous ligand replacement as well as the chemical composition of the mixed ligand-heterocycle loaded PCMOFs. Detailed compositional analysis on all PCMOFs is presented in the supplementary information. We have observed that higher phosphonate ligand loadings than 33% give new phases in the PXRD. AC Impedance Measurements. Powdered samples of PCMOF2½(Pz), PCMOF2½(Tz), PCMOF2(Pz), PCMOF2(Tz), (25 to 30 mg each) were placed in a glass cell and compressed between two solid titanium electrodes (0.3175 cm diameter). The sample length was measured by the difference between the empty cell and the filled cell, which was typically to 1 to 2 mm in length. The sample cells were placed inside a humidity and temperature controlled chamber (ESPEC BTL-433) and connected to a Princeton Applied Research VersaSTAT 3 impedance analyzer using a 2 probe setup. AC impedance 6 data was collected by cycling between 10 and 1 Hz with 10200 mV of applied potential using VersaStudio software. Samples were equilibrated for between 20 to 24 hours after each step in temperature (20 to 85 °C) and for 48 to 120 hours for each step in relative humidity (30% to 90% RH). These lengthy times are not arbitrary as short equilibration times, especially for RH, yield non-Arrhenius behavior and much lower conduction values. Exposure of the samples to humid environment was performed for conditions at and below 90% RH using the ESPEC BTL-433 humidity control oven. Proton conductivity of all samples was measured under 90% RH condition at various temperatures. The temperature was varied from 20 °C to 85 °C for minimum of two heating and cooling cycles with sufficient time for sample equilibration between each step. Following the heating and cooling cycles, the temperature was held at 25 °C and the humidity was decreased from 90% to the minimum of 30% to measure the humidity dependent proton conductivity of the samples. RESULT AND DISCUSSIONS Three new crystalline compounds, PCMOF2½(Pz), PCMOF2½(Tz) and PCMOF2(Pz) were synthesized and characterised through PXRD, elemental analysis, thermogravimetry and NMR analyses (See supporting information).In synthesizing PCMOF2½(Pz) and PCMOF2½(Tz), a three component mixed ligand MOF, it was critical to employ appropriate synthetic parameters to enable mutual compatibility between the two design strategies. At first the syntheses of PCMOF2(heterocycle) (heterocycle= pyrazole and tria-

ACS Paragon Plus Environment

Page 2 of 8

Page 3 of 8 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society zole) were carried out. This synthesis processes of PCMOF2(heterocycle) were tailored to kinetically capture pyrazole within the β-PCMOF2 framework via a solution precipitation. In the second design strategy PCMOF2(heterocycle) were then mixed with Na3H3L2 and placed under solid state reaction conditions to yield corresponding PCMOF2½(heterocycle) (Supporting information). Pelletization accelerates the process because it increases the contact between the multi-phases. Additionally, the high temperature and humidity condition accelerates the mobility of the molecules and drives the reaction because it increases the positional and vibrational entropy of the solid constituents. This synthetic approach is conceptually similar to that of Accelerated Aging which describes the natural mineral 60 weathering process. Figure 2 shows the powder XRD analysis of the various PCMOFs featured in this work with the focus on PCMOF2½(Pz) and PCMOF2½(Tz). We performed powder XRD analysis throughout the synthesis to confirm the preservation of the original β-PCMOF2 structure (Figure 2F). In synthesizing PCMOF2½(Pz) and PCMOF2½(Tz), finely ground powders of PCMOF2(Pz)/PCMOF2(Tz) and Na3H3L2 were mechanically mixed. The PXRD pattern of the mechanical mixture exhibited a simple super positioning of the patterns of its two components (Figure 2I, J & M and 2K, L & M). The mechanical mixture was then pressed into a pellet and placed under high temperature and humidity conditions

crease of ~0.1 & 0.1 Å and 0.03 & 0.04 Å in the d-spacing, respectively (Figure S11 and Table S1). The presence of hetero1 cycles were demonstrated by H-NMR studies, where the digested samples of PCMOF2(heterocycle) showed peaks corresponding to the heterocycle molecules in the framework (supporting information). The molecular formula of all the compounds were calculated merging elemental analysis and thermogravimetric analysis (supporting information). The maximum ratio of H3L2 that can be incorporated with retention of structure is 33%. While we do not have direct structural support, a possible arrangement based on this ratio would be a phosphonate ligand between two sulfonate ligands. To characterize the proton conductivity, AC impedance measurements were performed on samples of PCMOF2, PCMOF2(heterocycle) and PCMOF2½(heterocycle) at 90% RH in air. Table 1 summarises the conductivity values and the corresponding activation energies of the various PCMOFs under their respective conditions. The Nyquist plots obtained from the second heating cycle, to ensure humidity equilibration, are shown in Figure 4 and it shows that PCMOF2½(heterocycle) exhibited very low resistances, resulting in only the tail end of a semi-circle being observed at high frequencies due to instrument limitations. At lower frequencies, the capacitive tail is also observed, as expected for blocking effects of the mobile charge at the electrode interface. The conductivity was calculated from the real-axis intercept of the Nyquist plot (Figure 3). Table 1. List of PCMOF samples featured in this work and their conductivity profiles. Name

Conductivity Scm

Figure 2. Powder XRD patterns of (A) PCMOF2½(Pz), postimpedance; (B) PCMOF2½(Tz), post-impedance; (C) PCMOF2½, post-impedance; (D) PCMOF2(Pz), postimpedance; (E) PCMOF2(Tz), post-impedance; (F) βPCMOF2, post-impedance; (G) PCMOF2½(Pz), preimpedance; (H) PCMOF2½(Tz), pre-impedance; (I) Mechanical mixture of PCMOF2(Pz) and Na3H3L2 (Resembles J + M); (J) PCMOF2(Pz); (K) Mechanical mixture of PCMOF2(Tz) and Na3H3L2 (Resembles L + M); (L) PCMOF2(Tz) and (M) Na3H3L2. to achieve isomorphous ligand replacement, a thermody57 namically driven solid state reaction. The resultant powder XRD of PCMOF2½(Pz) and PCMOF2½(Tz) resembled that of β-PCMOF2 with a slight 2θ shift for the peaks that correspond to the (2,-1,0) plane and the (0,0,1) plane with a de-

1.1 × 10

PCMOF2½(Pz)

[i]

1.23 × 10

PCMOF2½(Pz)

[i]

7.2 × 10

PCMOF2½(Tz)

[i]

1.17 × 10

PCMOF2½(Tz)

[i]

-1

90

0.16

40

n/a

-7

150

0

0.20 [vii] 0.98

-1

85

90

0.22

-2

25

40

n/a

-2

85

90

0.10

1.08 × 10

PCMOF2(Pz)

3.28 × 10

PCMOF2(Pz)

[ii]

1.3 × 10

PCMOF2(Tz)

[ii]

1.9 × 10

PCMOF2(Tz)

[ii]

[iii]

4.6 × 10

-4

25

40

150

0

-2

85

90

0.10

3.18 × 10

25

40

n/a

-4

150

0

-2

85

90

0.21

-5

-5

5 × 10

2.1 × 10

2.4 × 10

85

90

0.28

1.8 × 10

-6

20

50

n/a

-9

100

0

n/a

β-PCMOF2

ACS Paragon Plus Environment

[ix]

n/a

[iii]

PCMOF2½(Im)

0.34

50

1.3 × 10

[i]

[viii]

20

[iv]

PCMOF2(Im)

n/a 0.42

-3

β-PCMOF2

[i]

[vi]

-3

[ii,iv]

β-PCMOF2

eV

25

[ii]

PCMOF2½

[v]

85

PCMOF2(Pz)

[iv]

Ea

-3

[ii]

[iv]

R H %

[i]

PCMOF2½

°C

-1

PCMOF2½(Pz)

PCMOF2(Tz)

T

1 × 10

These PCMOFs exhibited proton -6 -1 conductivity below 10 S cm at 85 °C and 90% RH.

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

i. New compounds and new conductivity data. ii. Compounds previously synthesized by Hurd et al; its conductivity 49 was measured under new conditions. iii. Compound and conductivity previously reported by Shimizu et al (the au49 thor). iv. Compounds and conductivity previously reported by 57 Shimizu et al (the author). v. Activation energy measured between the temperature of 85 °C and 25 °C. vi. Activation energy measured between the temperature of 150 °C and 110 °C. vii. Activation energy measured between the temperature of 110 °C and 90 °C. viii. Activation energy measured between the temperature of 150 °C and 60 °C.ix. Activation energy measured between the temperature of 150 °C and 90 °C. As a control, PCMOF2(Pz) and PCMOF2(Tz) were measured and their proton conductivities reproduced previously measured values on samples prepared by a different co-worker: -2 -1 -2 PCMOF2(Pz) (4.6 × 10 Scm ) and PCMOF2(Tz) (1.9 × 10 -1 Scm ) at 85 °C and 90% RH (Table 1). The mixed ligand, loaded samples were then measured and gave data almost an order of magnitude higher under identical conditions: 1.1 × -1 -1 -1 -1 10 Scm and 1.17× 10 Scm for PCMOF2½(Pz) and PCMOF2½(Tz), respectively, at 85 °C and 90% RH.

ty of the sample was maintained, as suggested by the reproducible conduction data, and corroborated by the pre- and post-impedance PXRD patterns (Figure 2). Notably, all PCMOFs investigated showed the same PXRD pattern as βPCMOF2 both before and after impedance measurements, confirming the preservation of the structure. The conductivity was found to be highly dependent on the humidity, which suggests that the water plays a critical role in proton conduction. Interestingly, both PCMOF2½(Pz) and PCMOF2½(Tz) showed a steady decrease in conductivity when the relative humidity was lowered from 90% to 40%, where as the conductivity of PCMOF2(Pz) and PCMOF(Tz) experienced a much more abrupt decrease (Table 1 & figure S8). The activation energy for proton transfer in the PCMOF2 58 series under humid condition is between 0.10-0.28 eV, which is consistent with a Grotthuss transfer mechanism (Figure 3). β-PCMOF2 had the greatest activation energy, 0.28 eV among all the PCMOFs investigated in this study. PCMOF2½ showed activation energy of 0.21 eV. PCMOF2(Tz) and PCMOF2(Pz) showed very low activation energies (0.10 eV). PCMOF2½(Pz) and PCMOF2½(Tz) showed activation energies of 0.16 and 0.22 eV respectively (Table 1 and figure S6-7).

Figure 4. Proton conductivity data (90% RH) for various samples. Lines are present only as visual guides

Figure 3. Nyquist plot for a) PCMOF2½(Pz) and b) PCMOF2½(Tz) at 90% RH. Since both PCMOF2½(Pz) and PCMOF2½(Tz) showed pro-1 -1 ton conducting properties over 10 Scm (at the time of initial measurement, the highest in any proton conducting MOF), new samples were remade (by a second researcher) and conductivity remeasured for four heating and cooling cycles to confirm stability and reproducibility. Each cycle was measured at 90% RH and from 20 °C to 85 °C with at least an 8 hour equilibration time between temperatures. The integri-

The activation energy (Ea) of proton transfer is dependent on + the facility of H transfer and the reorientation of proton 58 carriers post-transfer. Subsequently, proton mobility is a function of activation energy, and proton conductivity is a function of mobility. It may seem that a PCMOF with the lowest activation energy would yield the highest proton conductivity, but based on the results of this work (Figure 4 and Table 1), lowest activation energy does not necessarily equal to highest conductivity. In order to address this, multiple factors that affect conductivity – amphiprotic proton donor/acceptor sites, pKas, host/guest interactions, and the size of the heterocycles – need to be addressed since they all affect the resultant conductivity and the activation energy. Empirical comparison of conductivity can be made using the relative pKa values (Table 2) of the groups lining the pores and the conjugate acid form of the heterocycles.

ACS Paragon Plus Environment

Page 4 of 8

Page 5 of 8

Journal of the American Chemical Society -7

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Table 2. A list of pKas relevant to this work. Name

pKa

Benzenesulfonic acid

-0.60

Phenylphosphonic acid

1.88

1H-pyrazol-2-ium

2.50

1H-1,2,4-triazol-2-ium

2.39

1H-imidazol-3-ium

6.90

Hydroxonium

-1.74

-1

10 S cm at 150 °C in air, a much lower value than the anhydrous conductivity of PCMOF2(Pz) and PCMOF2(Tz), yet a -8 -1 higher value than that of pure β-PCMOF2 (10 S cm ) (Figure S11). It is speculated that the loaded heterocycle in PCMOF2½(Pz) still imparts some degree of enhanced anhydrous enhancement but nowhere as drastic as the 5 orders of magnitude increase as PCMOF2(Pz) versus β-PCMOF2. We hypothesized that the conductivity enhancement by Na3H3L2 requires sufficient amount of pore water molecules to activate/dissociate its free acidic protons. This is supported by the rapid decrease in PCMOF2½(Pz)’s and PCMOF2½(Tz)’s conductivity when the humidity was decreased from 90% to 40% RH.

PCMOF2(Pz) and PCMOF2(Tz) has an activation energy for proton transfer of 0.10 eV, a very low value indicative of highly facile proton transport pathway. PCMOF2½(Pz) and PCMOF2½(Tz), on the other hand, have activation energies of 0.16 eV and 0.22 eV, respectively, double that of PCMOF2(Pz) and PCMOF2(Tz). PCMOF2 itself has a higher 57 activation energy than PCMOF2½. This result is in accord with the pKa values in Table 2 as the sulfonate group has lower affinity for protons. However, despite the higher activation energy, in all three cases, the PCMOF2½(heterocycle) complexes show higher conductivity than the pure sulfonate counterparts (PCMOF2(heterocycle)) . Revisiting equation [1], the conductivity in the PCMOF2½(heterocycle) series is greater because the higher activation energy is offset by the increase in the number of freely available acidic protons resulting in enhanced conductivity. Concurrently, the greater difference in pKa will also lead to greater host-guest interaction between Na3H3L2 and the heterocycles, inhibiting reorientation and increasing the activation energy. pKa is an empirical solution measurement and may be a good predictor for the resultant conductivity but in the PCMOF framework, other factors such as their size, hydrogen bonding proficiency and host/guest interactions (which in turn affects the motional entropy and the rotational freedom of the heterocycles) also impact the resultant mobility and conductivity. Being similar in terms of pKa and molecular size, both triazole & pyrazole showed similar conductivity. Perhaps this is because of both the pKa of pyrazolium and triazolium cation is much closer to the pKa of hydronium cation (-1.74), enabling a more facile reversible proton transfer between the two conjugate acids. The pKa values indicates that the triazole-triazolium or pyrazole-pyrazolium pair would form a more balanced mix of protonated and unprotonated molecules – it is just as important to provide sufficient unprotonated sites that can act as facile proton acceptors in addition to maximizing the number of proton donors.

The field of proton conducting MOF materials has seen steady progress in the levels of proton conduction being re-1 -1 45 ported now with accounts of values exceeding 10 Scm . These advances are enabled by design advances in parallel with increasingly more robust MOF structures being reported. In this case, the MOF materials are not highly robust and would not survive a long-term fuel cell test; the value of this work is the design principles put forth. The fundamental chemistry governing proton transfer is identical regardless of the specific material or even class of solid. There must be a continuous proton transfer pathway with minimal potential wells to trap the protons. MOFs, with their crystalline structures offer the ability to correlate crystallographic structure and function. This coupled with the well-established routes for systematic structural variation, as we have shown here, make iterative cycles of design and assessment a rational rather than stochastic process.

As a contrast, we also investigated imidazole (Im, pKa1 = 7.18, pKa2 = 14.52) as a dopant. We synthesized PCMOF2(Im) and PCMOF2½(Im) and measured their respective conductivities. Neither sample showed significant proton conductivity at 85 °C and 90% RH. We hypothesized the minimal conductivity to be caused by the much higher pKa of imidazolium meaning much more localized protons as well as the resulting stronger host-guest interaction between the imidazole molecule and the MOF framework limiting reorientation.

Details of synthesis, powder XRD, impedance analysis. This material is available free of charge via the Internet at http://pubs.acs.org.

Comparisons of the proton conductivity in these systems was also carried out under anhydrous conditions. PCMOF2½(Pz) was dehydrated and the conductivity was found to be 7.2 x

The authors declare no competing financial interest.

CONCLUSION We have successfully applied two design strategies - isomorphous ligand replacement and heterocycle doping - to enhance the proton conductivity of a proton conducting metalorganic framework: β-PCMOF2. The rational and sequential variation of structures gives proton conductivity data that correlates with expected performance based on molecular properties. We have also shown a solid state route to controllably loading 1-D pores in the MOFs with proton carriers. Of six PCMOFs investigated in this study, two materials, PCMOF2½(Pz) and PCMOF2½(Tz), show exceptionally high -1 -1 proton conduction values (>1 × 10 Scm ) at 85 °C and 90% RH, while maintaining the parent MOF structure. The principles outlined are applicable to other proton conducting systems beyond MOF materials.

ASSOCIATED CONTENT

AUTHOR INFORMATION Corresponding Author [email protected]

Notes

ACKNOWLEDGMENT

ACS Paragon Plus Environment

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

We thank the Natural Science and Engineering Research Council (NSERC) of Canada for support of this research and for graduate scholarships to S. K. and B. S. G.. We thank Alberta Innovates Technology Futures for a doctoral scholarship to K.W.D. and a postdoctoral fellowship to B. J.

REFERENCES (1) Yamada, T.; Otsubo, K.; Makiura, R.; Kitagawa, H. Chem. Soc. Rev. 2013, 42, 6655. (2) Ramaswamy, P.; Wong, N. E.; Shimizu, G. K. H. Chem. Soc. Rev. 2014, 43, 5913. (3) Meng, X.; Wang, H.-N.; Song, S.-Y.; Zhang, H.-J. Chem. Soc. Rev. 2017, 46, 464. (4) Ramaswamy, P.; Wong, N. E.; Gelfand, B. S.; Shimizu, G. K. H. J. Am. Chem. Soc. 2015, 137, 7640. (5) Sadakiyo, M.; Yamada, T.; Honda, K.; Matsui, H.; Kitagawa, H. J. Am. Chem. Soc. 2014, 136, 7701. (6) Yoon, M.; Suh, K.; Natarajan, S.; Kim, K. Angew. Chem., Int. Ed. 2013, 52, 2688. (7) Tominaka, S.; Coudert, F.-X.; Dao, T. D.; Nagao, T.; Cheetham, A. K. J. Am. Chem. Soc. 2015, 137, 6428. (8) Shimizu, G. K. H.; Taylor, J. M.; Kim, S. Science 2013, 341, 354. (9) Nguyen, N. T. T.; Furukawa, H.; Gándara, F.; Trickett, C. A.; Jeong, H. M.; Cordova, K. E.; Yaghi, O. M. J. Am. Chem. Soc. 2015, 137, 15394. (10) Jeong, N. C.; Samanta, B.; Lee, C. Y.; Farha, O. K.; Hupp, J. T. J. Am. Chem. Soc. 2012, 134, 51. (11) Cai, K.; Sun, F.; Liang, X.; Liu, C.; Zhao, N.; Zou, X.; Zhu, G. J. Mater. Chem. A 2017, 5, 12943. (12) Bazaga-García, M.; Colodrero, R. M. P.; Papadaki, M.; Garczarek, P.; Zoń, J.; Olivera-Pastor, P.; Losilla, E. R.; LeónReina, L.; Aranda, M. A. G.; Choquesillo-Lazarte, D.; Demadis, K. D.; Cabeza, A. J. Am. Chem. Soc. 2014, 136, 5731. (13) Pili, S.; Argent, S. P.; Morris, C. G.; Rought, P.; GarcíaSakai, V.; Silverwood, I. P.; Easun, T. L.; Li, M.; Warren, M. R.; Murray, C. A.; Tang, C. C.; Yang, S.; Schröder, M. J. Am. Chem. Soc. 2016, 138, 6352. (14) Taylor, J. M.; Komatsu, T.; Dekura, S.; Otsubo, K.; Takata, M.; Kitagawa, H. J. Am. Chem. Soc. 2015, 137, 11498. (15) Furukawa, H.; Cordova, K. E.; O’Keeffe, M.; Yaghi, O. M. Science 2013, 341. (16) Choi, K. M.; Jeong, H. M.; Park, J. H.; Zhang, Y.-B.; Kang, J. K.; Yaghi, O. M. ACS Nano 2014, 8, 7451. (17) Meng, X.; Song, X.-Z.; Song, S.-Y.; Yang, G.-C.; Zhu, M.; Hao, Z.-M.; Zhao, S.-N.; Zhang, H.-J. Chem. Commun. 2013, 49, 8483. (18) Borges, D. D.; Devautour-Vinot, S.; Jobic, H.; Ollivier, J.; Nouar, F.; Semino, R.; Devic, T.; Serre, C.; Paesani, F.; Maurin, G. Angew. Chem., Int. Ed. 2016, 55, 3919. (19) Li, D.-S.; Wu, Y.-P.; Zhao, J.; Zhang, J.; Lu, J. Y. Coord. Chem. Rev. 2014, 261, 1. (20) Wu, Y.-P.; Zhou, W.; Zhao, J.; Dong, W.-W.; Lan, Y.-Q.; Li, D.-S.; Sun, C.; Bu, X. Angew. Chem., Int. Ed. 2017, 56, 13001. (21) Zhai, Q.-G.; Bu, X.; Zhao, X.; Li, D.-S.; Feng, P. Acc. Chem. Res. 2017, 50, 407. (22) Pili, S.; Argent, S. P.; Morris, C. G.; Rought, P.; GarciaSakai, V.; Silverwood, I. P.; Easun, T. L.; Li, M.; Warren, M.

R.; Murray, C. A.; Tang, C. C.; Yang, S. H.; Schroder, M. J. Am. Chem. Soc. 2016, 138, 6352. (23) Zhai, Q.-G.; Mao, C.; Zhao, X.; Lin, Q.; Bu, F.; Chen, X.; Bu, X.; Feng, P. Angew. Chem., Int. Ed. 2015, 54, 7886. (24) Stankiewicz, J.; Tomas, M.; Dobrinovitch, I. T.; ForcenVazquez, E.; Falvello, L. R. Chem. Mat. 2014, 26, 5282. (25) Phang, W. J.; Lee, W. R.; Yoo, K.; Ryu, D. W.; Kim, B.; Hong, C. S. Angew. Chem., Int. Ed. 2014, 53, 8383. (26) Nagarkar, S. S.; Unni, S. M.; Sharma, A.; Kurungot, S.; Ghosh, S. K. Angew. Chem., Int. Ed. 2014, 53, 2638. (27) Chandra, S.; Kundu, T.; Kandambeth, S.; BabaRao, R.; Marathe, Y.; Kunjir, S. M.; Banerjee, R. J. Am. Chem. Soc. 2014, 136, 6570. (28) Okawa, H.; Sadakiyo, M.; Yamada, T.; Maesato, M.; Ohba, M.; Kitagawa, H. J. Am. Chem. Soc. 2013, 135, 2256. (29) Umeyama, D.; Horike, S.; Inukai, M.; Itakura, T.; Kitagawa, S. J. Am. Chem. Soc. 2012, 134, 12780. (30) Mallick, A.; Kundu, T.; Banerjee, R. Chemical communications (Cambridge, England) 2012, 48, 8829. (31) Wei, Y.-S.; Hu, X.-P.; Han, Z.; Dong, X.-Y.; Zang, S.-Q.; Mak, T. C. W. J. Am. Chem. Soc. 2017, 139, 3505. (32) Li, R.; Wang, S.-H.; Chen, X.-X.; Lu, J.; Fu, Z.-H.; Li, Y.; Xu, G.; Zheng, F.-K.; Guo, G.-C. Chem. Mat. 2017, 29, 2321. (33) Joarder, B.; Lin, J.-B.; Romero, Z.; Shimizu, G. K. H. J. Am. Chem. Soc. 2017, 139, 7176. (34) Zhang, G.; Fei, H. Chem. Commun. 2017, 53, 4156. (35) Nagarkar, S. S.; Horike, S.; Itakura, T.; Le Ouay, B.; Demessence, A.; Tsujimoto, M.; Kitagawa, S. Angew. Chem., Int. Ed. 2017, 56, 4976. (36) Ye, Y.; Zhang, L.; Peng, Q.; Wang, G.-E.; Shen, Y.; Li, Z.; Wang, L.; Ma, X.; Chen, Q.-H.; Zhang, Z.; Xiang, S. J. Am. Chem. Soc. 2015, 137, 913. (37) Nagao, Y.; Kubo, T.; Nakasuji, K.; Ikeda, R.; Kojima, T.; Kitagawa, H. Synth. Met. 2005, 154, 89. (38) Yamada, T.; Sadakiyo, M.; Kitagawa, H. J. Am. Chem. Soc. 2009, 131, 3144. (39) O̅kawa, H.; Shigematsu, A.; Sadakiyo, M.; Miyagawa, T.; Yoneda, K.; Ohba, M.; Kitagawa, H. J. Am. Chem. Soc. 2009, 131, 13516. (40) Sahoo, S. C.; Kundu, T.; Banerjee, R. J. Am. Chem. Soc. 2011, 133, 17950. (41) Pardo, E.; Train, C.; Gontard, G.; Boubekeur, K.; Fabelo, O.; Liu, H.; Dkhil, B.; Lloret, F.; Nakagawa, K.; Tokoro, H.; Ohkoshi, S.-i.; Verdaguer, M. J. Am. Chem. Soc. 2011, 133, 15328. (42) Sadakiyo, M.; O̅kawa, H.; Shigematsu, A.; Ohba, M.; Yamada, T.; Kitagawa, H. J. Am. Chem. Soc. 2012, 134, 5472. (43) Colodrero, R. M. P.; Olivera-Pastor, P.; Losilla, E. R.; Aranda, M. A. G.; Leon-Reina, L.; Papadaki, M.; McKinlay, A. C.; Morris, R. E.; Demadis, K. D.; Cabeza, A. Dalton Trans. 2012, 41, 4045. (44) Colodrero, R. M. P.; Papathanasiou, K. E.; Stavgianoudaki, N.; Olivera-Pastor, P.; Losilla, E. R.; Aranda, M. A. G.; León-Reina, L.; Sanz, J.; Sobrados, I.; ChoquesilloLazarte, D.; García-Ruiz, J. M.; Atienzar, P.; Rey, F.; Demadis, K. D.; Cabeza, A. Chem. Mater. 2012, 24, 3780. (45) Yang, F.; Xu, G.; Dou, Y.; Wang, B.; Zhang, H.; Wu, H.; Zhou, W.; Li, J.-R.; Chen, B. Nature Energy 2017. (46) Zhang, F.-M.; Dong, L.-Z.; Qin, J.-S.; Guan, W.; Liu, J.; Li, S.-L.; Lu, M.; Lan, Y.-Q.; Su, Z.-M.; Zhou, H.-C. J. Am. Chem. Soc. 2017, 139, 6183. (47) Wu, H.; Yang, F.; Lv, X.-L.; Wang, B.; Zhang, Y.-Z.; Zhao, M.-J.; Li, J.-R. J. Mater. Chem. A 2017, 5, 14525.

ACS Paragon Plus Environment

Page 6 of 8

Page 7 of 8 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society (48) Inukai, M.; Horike, S.; Itakura, T.; Shinozaki, R.; Ogiwara, N.; Umeyama, D.; Nagarkar, S.; Nishiyama, Y.; Malon, M.; Hayashi, A.; Ohhara, T.; Kiyanagi, R.; Kitagawa, S. J. Am. Chem. Soc. 2016, 138, 8505. (49) Hurd, J. A.; Vaidhyanathan, R.; Thangadurai, V.; Ratcliffe, C. I.; Moudrakovski, I. L.; Shimizu, G. K. H. Nat. Chem. 2009, 1, 705. (50) Sadakiyo, M.; Yamada, T.; Kitagawa, H. J. Am. Chem. Soc. 2009, 131, 9906. (51) Ramaswamy, P.; Matsuda, R.; Kosaka, W.; Akiyama, G.; Jeon, H. J.; Kitagawa, S. Chem. Commun. 2014, 50, 1144. (52) Horike, S.; Umeyama, D.; Inukai, M.; Itakura, T.; Kitagawa, S. J. Am. Chem. Soc. 2012, 134, 7612. (53) Bureekaew, S.; Horike, S.; Higuchi, M.; Mizuno, M.; Kawamura, T.; Tanaka, D.; Yanai, N.; Kitagawa, S. Nat. Mater. 2009, 8, 831. (54) Wang, K.; Jin, Y.; Jiang, L.; Wang, Z.; Zhang, Q. CrystEngComm 2017, 19, 3997.

(55) Homburg, T.; Hartwig, C.; Reinsch, H.; Wark, M.; Stock, N. Dalton Trans. 2016, 45, 15041. (56) Kang, D. W.; Lim, K. S.; Lee, K. J.; Lee, J. H.; Lee, W. R.; Song, J. H.; Yeom, K. H.; Kim, J. Y.; Hong, C. S. Angew. Chem., Int. Ed. 2016, 55, 16123. (57) Kim, S.; Dawson, K. W.; Gelfand, B. S.; Taylor, J. M.; Shimizu, G. K. H. J. Am. Chem. Soc. 2013, 135, 963. (58) Shigematsu, A.; Yamada, T.; Kitagawa, H. J. Am. Chem. Soc. 2011, 133, 2034. (59) Taylor, J. M.; Mah, R. K.; Moudrakovski, I. L.; Ratcliffe, C. I.; Vaidhyanathan, R.; Shimizu, G. K. H. J. Am. Chem. Soc. 2010, 132, 14055-14057. (60) Mottillo, C.; Lu, Y.; Pham, M.-H.; Cliffe, M. J.; Do, T.O.; Friscic, T. Green Chem. 2013, 15, 2121.

ACS Paragon Plus Environment

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 8

TOC

ACS Paragon Plus Environment

8