ACS Publications. Most Trusted. Most Cited. Most Read
Effects of Ligand Substituents on the Character of Zn-Coordination in Zeolitic Imidazolate Frameworks
My Activity

Figure 1Loading Img
  • Open Access
C: Physical Properties of Materials and Interfaces

Effects of Ligand Substituents on the Character of Zn-Coordination in Zeolitic Imidazolate Frameworks
Click to copy article linkArticle link copied!

Open PDFSupporting Information (1)

The Journal of Physical Chemistry C

Cite this: J. Phys. Chem. C 2023, 127, 43, 21456–21464
Click to copy citationCitation copied!
https://doi.org/10.1021/acs.jpcc.3c06054
Published October 25, 2023

Copyright © 2023 The Authors. Published by American Chemical Society. This publication is licensed under

CC-BY 4.0 .

Abstract

Click to copy section linkSection link copied!

Due to their favorable properties and high porosity, zeolitic imidazolate frameworks (ZIFs) have recently received much limelight for key technologies such as energy storage, optoelectronics, sensorics, and catalysis. Despite widespread interest in these materials, fundamental questions regarding the zinc coordination environment remain poorly understood. By focusing on zinc(II)2-methylimidazolate (ZIF-8) and its tetrahedrally coordinated analogues with Br-, Cl-, and H-substitution in the 2-ring position, we aim to clarify how variations in the local environment of Zn impact the charge distribution and the electronic properties of these materials. Our results from density-functional theory confirm the presence of a Zn coordinative bond with a large polarization that is quantitatively affected by different substituents on the organic ligand. Moreover, our findings suggest that the variations in the Zn coordination induced by the functionalization have a negligible effect on the electronic structure of the considered compounds. On the other hand, halogen terminations of the ligands lead to distinct electronic contributions in the vicinity of the frontier region which ultimately reduce the band gap size by a few hundred millielectron volts. Experimental results obtained from X-ray absorption spectroscopy (Zn K-edge) confirm the trends predicted by theory and, together with them, contribute to a better understanding of the structure–property relationships that are needed to tailor ZIFs for target applications.

This publication is licensed under

CC-BY 4.0 .
  • cc licence
  • by licence
Copyright © 2023 The Authors. Published by American Chemical Society

Introduction

Click to copy section linkSection link copied!

Zeolitic imidazolate frameworks (ZIFs) are a popular family of metal–organic frameworks (MOFs) and ZIF-8, a polymorph of zinc(II)2-methylimidazolate, (1,2) is their most prominent member. Although transition-metal imidazolates have been known for decades, (3) these materials have entered the spotlight after the discovery of their high porosity (4,5) and exceptional thermal and chemical stability. (1,2) Thanks to their straightforward synthesis (6,7) and characterization (8) as well as their intriguing properties (9−11) even under pressure (12−17) and in response to external stimuli, (18) ZIFs have become relevant for a number of technological areas ranging from gas storage to photocatalysis, and from biomedicine to optoelectronics. (19−26)
From a fundamental perspective, these characteristics are ruled by the electronic charge distribution in the material and the nature of the chemical bonds therein. Despite the general consensus regarding the ionic nature of the Zn–N bond in ZIFs, (27−30) the actual strength of this bond and the influence exerted on it by the chemical environment are still a matter of debate. (15,31,32) In particular, recent studies have unraveled nontrivial characteristics in the electronic structure of ZIF-8, (33−35) which point to the existence of some degree of covalence in the bond between Zn and N. This hypothesis is supported by ab initio calculations on the transport properties (33) and by the vibrational activity of ZIFs. (34,35) Experimental evidence of the (partial) covalent character of Zn–N coordination is, to date, only indirect. Observations of chemical bond stability of ZIF-8 and related compounds especially in acidic environments point to their susceptibility to hydrolytic cleavage, which is only known for covalent bonds. (36) Likewise, the photocatalytic activity of ZIFs is sensitive to structural and chemical changes within the framework, (37) while the polarizability of metallic ions should be hindered in a structure with purely ionic coordination. In order to clarify the nature of Zn–N coordination in ZIF-8, a dedicated study combining both theoretical and experimental analyses is urgently needed.
Here, we present a joint investigation of ZIF-8 and related compounds with ligand substitution carried out with first-principles calculations and X-ray spectroscopy. In this analysis, we consider halogen ligands terminations such as Cl and Br [zinc(II)2-chloroimidazole and zinc(II)2-bromoimidazole, abbreviated hereafter as ZIF-8-Cl and ZIF-8-Br], which are more electronegative than the CH3 group in zinc(II)2-methylimidazolate (ZIF-8); additionally, we include H termination in zinc(II)imidazolate (ZIF-8-H), which is expected to be slightly more electropositive than ZIF-8. By means of density functional theory (DFT), we investigate the charge distribution of these compounds both in their (simplified) crystal structure as well as in the isolated metal–organic complex building up the framework and in the isolated imidazole unit. This way, we assess the role of the local chemical environment of Zn coordination in determining the charge distribution in these systems. Furthermore, through analysis of the density of states, we clarify the influence of the functional groups on the character and composition of the electronic states close to the frontier. We complement this study with the analysis of the near-edge X-ray absorption spectrum of ZIF-8 and its Br-terminated analogue measured from the Zn K-edge.

Methods

Click to copy section linkSection link copied!

Modeling ZIF-8 and Its Functionalized Analogues

In this study, we focus on ZIF-8, an established MOF structure with methyl-functionalized linker molecules. To investigate the effects induced in the Zn coordination by different terminations of the organic ligands, we consider the ZIF-8 analogues with H-, Cl-, and Br-substituents of the methyl groups in the 2-ring position (see Figure 1). While the H-passivated ZIF-8 has been only predicted computationally, (38) the halogen-functionalized systems have been synthesized by several groups (see, among others, refs (39)(42)) and adopted in theoretical studies as representative ZIF-8 variants. (40,43,44) To model ZIF-8-H, we have adopted the predicted crystal structure published in ref (38). On the other hand, we have taken the structural information on ZIF-8 and its halogenated counterparts directly from the experimental references.

Figure 1

Figure 1. Ball-and-stick representation of the ZIF-8 compounds considered in this work including the (a) crystal structure with R3 space-group symmetry, (b) isolated metal–organic complex replicated in the crystal, and (c) imidazole ligand. Carbon atoms are depicted in black, hydrogen atoms in white, nitrogen atoms in blue, zinc atoms in gray, and the terminating group in green (CH3 in ZIF-8, H, Cl, and Br in the substituted siblings). The atoms in the ligand are marked in panel c: identical labels indicate equivalent atoms; Ha is the proton added to the imidazole unit and to the complex to passivate one of the two N atoms.

For ZIF-8, we have considered the system with the R3 space group. (45) With a script developed in-house, we have adjusted the structure provided in ref (45) by removing redundant H atoms. The adopted procedure includes finding the nearest hydrogen for each substituted C atom of the methyl group and the next H atoms with angles comprised between 100 and 130° to obtain the CH3 functional group. For both ZIF-8-Cl and ZIF-8-Br, we took the structural data from ref (41). In these cases, the linker molecule assumes two different orientations in the scaffold. In our bulk models, we took the molecular orientation giving rise to a bulk with 43m symmetry. All of these crystals have the same network topology (SOD). (46) The rhombohedral unit cell adopted in the calculations of all the considered ZIF-8-X (with X = CH3, H, Br, Cl), see Figure 1a, is retrieved from the initial cubic cell using the code seek-path. (47) To preserve the symmetry of the initial structures, the input lattice parameters of all of the considered ZIF-8 variants (see Table S1 of the Supporting Information) are kept constant during optimization, which thus entails only the minimization of the interatomic forces.
To model the isolated complexes (Figure 1b), the atomic positions are extracted from the ZIF-8-H crystal. (38) Based on this structure, the other functional groups are placed to substitute for the terminating hydrogen. The same procedure is adopted for the imidazole molecule (Figure 1c). In these calculations, the nonperiodic structures are embedded in a supercell including 7.5 Å of vacuum in each direction to prevent spurious interactions among periodic replicas.

Computational Details

All calculations presented in this work are performed in the framework of DFT (48) solving the Kohn–Sham equations (49) with the Perdew–Burke–Ernzerhof (PBE) exchange–correlation functional. (50) The code Quantum ESPRESSO (51) is adopted with the projector augmented wave method (52) and pseudopotentials from the ps-library. (53) A 2 × 2 × 2 k-mesh is used to sample the Brillouin zone of all of the considered periodic structures. The plane-wave and electron density cutoffs are set to 150 and 1000 Ry, respectively. Structure optimization is performed by minimizing interatomic forces until the convergence threshold of 10–4 Ry/bohr in the crystal and of 10–3 Ry/bohr is achieved in the metal–organic complexes and in the imidazole units, respectively. The Bader partial charges (54) and the planar deformation density plots are calculated using the postprocessing software package critic2. (55) The method derived by Yu-Trinkle (56) is used to integrate the partial charges.

Experimental Methods

X-ray absorption near-edge structure spectroscopy (XANES) measurements were performed at the BAMline (BESSY-II, Helmholtz Centre Berlin for Materials and Energy, Berlin, Germany). (57) The beam was monochromatized using a double-crystal monochromator (DCM) with a Si crystal with [111] orientation. The size of the beam was 3 mm (h) × 1 mm (v). The measurements were performed at Zn K-edge (9659 eV) in transmission, with two ionization chambers as detectors. For the pre-edge region, the energy was varied in 10 eV steps; for the region around the edge, the energy was tuned first in 0.5 eV steps and from then on in 1 eV steps. For the measurement, the samples were mixed with boron nitride, placed in polycarbonate hole plates with a thickness of 1 mm, and sealed with polyimide tape (Kapton) on both sides. Before collecting the sample spectra, a zinc foil was used as a reference for the respective K edge. The relative energies of the spectra were calibrated against the first inflection point from the first derivative of the zinc metal absorption edge. The XAS data were processed by using ATHENA which belongs to the main package IFEFFIT (v. 1.2.11). (58)

Results and Discussion

Click to copy section linkSection link copied!

Structural Properties

We start our study by inspecting the structural properties of the considered systems. The gained information is important to support the charge-density analysis reported in the next section. The results obtained for the Zn–N bond in all of the considered ZIFs are reported in Table 1, where the initial values in the unrelaxed input structures are included for comparison. The errors associated with the values of the optimized structures are due to the averaging. The separation between Zn and N is close to 2 Å in all crystal structures with variations induced by the ligand substituents. In ZIF-8, the calculated Zn–N distance is equal to 1.989 Å, slightly larger than in the initial experimental structure (45) and in excellent agreement with other measurements. (1) Ligand substitution leads to an increase in the Zn–N bond length by less than 0.01 Å upon halogen functionalization and by 0.015 Å in the presence of H. Notice the different trend of the initial structures, where the Zn–N separation is smaller in ZIF-8-Cl than in ZIF-8. ZIF-8-H is the only structure in which the Zn–N distance decreases upon relaxation (see Table 1). The other bond lengths undergo variations of the same magnitude as Zn–N in all of the considered systems (see Table S3). In particular, it is worth mentioning that in ZIF-8 and its halogen-functionalized counterparts, the distance between the carbon atoms C2 (see Figure 1c) increases upon relaxation, while it remains essentially equal to the experimental value in ZIF-8-H. Likewise, the Cn-N distances (n = 1 and 2, see Figure 1c), which are almost identical in the initial structures, rearrange themselves upon relaxation.
Table 1. Average Zn–N Bond Lengths of the Initial Structures and of the Optimized Crystals and Complexesa
 ZIF-8 (Å)ZIF-8-H (Å)ZIF-8-Cl (Å)ZIF-8-Br (Å)
initial1.9652.0251.9421.987
crystal1.989(0.003)2.004(0.001)1.997(0.001)1.998(0.001)
complex2.031(0.089)2.019(0.086)2.022(0.087)2.029(0.085)
a

The uncertainties due to the averaging reported in parenthesis.

The same analysis performed on the isolated metal–organic complexes reveals a different trend (see Table 1). ZIF-8 is characterized by the longest Zn–N bond (2.031 Å), followed by those in the halogenated compounds (2.029 Å with Br and 2.022 Å with Cl). The complex terminated with H exhibits the shortest Zn–N distance of 2.019 Å. These findings, in excellent agreement with earlier results obtained on such clusters on an analogous level of theory, (27) indicate that in the absence of any structural constraints given by the periodic boundary conditions, the ligand termination in the metal–organic complexes exerts a more pronounced influence on the local environment of the Zn coordination. In particular, the halogen species, characterized by a significantly larger atomic radius and electronegativity than hydrogen, give rise to an increase in the Zn–N bond length, while in the bulk, the opposite trend is found. In elaborating on these values, two important points should be considered. First, the Zn–N bond lengths reported for the complexes carry an uncertainty that is almost 2 orders of magnitude higher than their counterparts in the crystals. This is due to the larger variations in the nonperiodic system among the four Zn–N bonds that are present therein. Second, the atoms in the complexes are completely free to move in the DFT optimization, in contrast with those in the crystals, which are constrained by the periodicity of the lattice. In this regard, it is worth recalling that no volume optimization was performed for the crystal structures: the atoms were allowed to relax within the experimental unit cells taken as input, with the symmetry of the synthesized crystals being preserved.
Overall, the analysis of all bond lengths in the crystals and complexes reveals that only the Zn–N distances undergo significant variations (Table S2), thus confirming that the influence of ligand functionalization is partly counterbalanced in the periodic environment. Our choice to study isolated metal–organic complexes as they are extracted from the crystalline scaffold, i.e., without any additional “capping” of the N bonds to mimic the dative nature of the Zn–N bond in the crystal, (59) is motivated by our interest in understanding how the Zn coordination is affected by ligand terminations in isolated and periodic arrangements.
The results obtained in this study for the considered ZIF-8 structures are in general agreement with the experimental Zn–N bond lengths reported in the literature for various organometallic compounds. (60,61) Specifically, in the crystal structure of hexaimidazole-zinc(II) dichloride tetrahydrate, where Zn is coordinated with N, the Zn–N distances range between 2.153 and 2.264 Å. (62) Shorter values are reported for dichlorobis(N-n-propylsalicylaldimine)zinc(II), where Zn is twofold coordinated with N, and Zn–N separations are slightly below 2 Å (1.998 Å). (63) This variability of the Zn–N bond length in different chemical environments reflects, on the one hand, the current discussions regarding the nature of such a bond and, on the other hand, highlights the need for a deeper investigation.

Partial Charge Analysis

Equipped with the knowledge of the structural properties presented above, we proceed with the partial charge analysis. This way, we aim to gain a better understanding of the electron density distribution around the Zn–N bond and the effects of ligand functionalization. The adopted Bader scheme, (54) in particular, offers reliable results for the periodic systems and the isolated clusters on equal footing. With this method, we evaluate the relative amount of electronic charge on the atomic species of ZIF-8 and its related compounds with ligand substitution. The results reported in Figure 2 correspond to the partial charges averaged on all atoms of a given species within the unit cell adopted in the simulations (for further details, see Supporting Information, Tables S4–S8). In ZIF-8, an excess of positive charge (1.17 e) is found on the metal ion while the (negative) electronic charge on nitrogen increases by −1.21 e compared to its nominal value in the free atom. The difference between the absolute values of these quantities is compensated by the functional group (CH3 in ZIF-8) which takes up a positive charge of 0.1 e. The remaining charge density has a positive sign and is distributed among the C and H atoms. Carbon atoms bound to the functional group are labeled as C1 in Figure 1 and host a sizable amount of charge equal to 0.93 e. Conversely, the other two equivalent C atoms in the imidazole ring (C2) accommodate only 0.352 e each. The remaining ∼0.1 e is found on the H atoms added for passivation (see Figure 2). These trends can be understood by considering that C1 shares two bonds with the electronegative N atoms in contrast with C2, which only has one.

Figure 2

Figure 2. Partial charges per atom in the considered ZIF structures. Results reported for the “additional H” concern only the complexes and isolated imidazole units, where protonation of a dangling bond was needed to neutralize the system.

Moving now to the substituted compounds, we notice the effect of the electronegativity of the functional group in redistributing the charge density and hence in modifying the character of the Zn–N coordination. When the CH3 group is replaced by H, the positive charge on Zn increases by 0.01 e, while the negative charge on each N atom remains almost unaltered with respect to ZIF-8 (see Figure 2). The difference is covered by the H atom, which donates 70% of its electronic charge to the network, and by the C2 atoms, which increment their positive charge up to 0.358 e. Upon halogen functionalization, the picture changes significantly. In the presence of Cl termination, where Cl is the most electrophilic element among the considered functional groups, the excess of positive charge on Zn becomes almost as high as 1.20 e, while the negative charge on N decreases by 0.04 e with respect to the reference value in ZIF-8 (see Figure 2). This result shows that this behavior is not driven by an increased ionicity of the Zn–N coordination but is rather an effect of the substituents that modify the chemical environment of this bond. The Cl atom takes up −0.2 e and induces a charge imbalance also on the C and H atoms, which all become more positively charged. Notice that the partial charge in C1, the carbon atom sharing a bond with Cl, is equal to 1.08 e, while on C2, it is as high as 0.368 e. The increase of the positive charge on H is a consequence of this redistribution, too (see Figure 2). Finally, Br termination induces similar results as Cl as far as the Zn–N bond is concerned: In this system, the metal ion hosts a positive charge of about 1.2 e and nitrogen a negative charge of a very similar magnitude (−1.18 e). However, the lower electronegativity of Br compared to Cl leads to an almost vanishing partial charge on Br and to a charge distribution on C1 that is very close to the one in ZIF-8 and ZIF-8-H (Figure 2). For C2 and especially for the H atoms, the values of partial charges are close to those obtained with Cl termination but the small absolute values should be noted.
For a deeper understanding of these findings, it is interesting to evaluate the impact of crystalline coordination in the charge distribution around the Zn–N bond. For this purpose, we examine the results of the partial charge analysis performed on the isolated complex inside the unit cell of the considered structures (see Figure 1b). The absolute values of the partial charges δ on the functional groups are almost identical to those obtained for the periodic structures (see Figure 2). The substantial similarities among the computed values of δ on C1 in the crystals and complexes confirm that their bonds formed with the functional groups are determined only by the charge distribution among the involved species. In contrast, the partial charges on Zn and N are subject to substantial differences in the isolated clusters compared with the crystalline arrangements, in agreement with the variations of their distances discussed above (see Table 1). In all complexes, regardless of the functional groups, the values of positive and negative charges on Zn and N, respectively, decrease systematically with respect to their counterparts in the crystals (Figure 2). This variation can be understood by the passivation of the N dangling bonds. The additional H atoms host a positive partial charge on the order of 0.4 e, which is maximized by the presence of the halogen terminations and minimized in ZIF-8.
While this analysis suggests a reduction of Zn–N bond polarization in the isolated clusters compared to the crystals, the inspection of the Zn–N bond lengths computed for the complexes points in the opposite direction (see Table 1). The Zn–N distances increase in the nonperiodic systems exceeding in all cases 2 Å. In the complexes, the shorter Zn–N distance is found in the presence of H termination and the longest one is with methyl. Despite the larger uncertainties obtained for these values compared with their counterparts in the crystals (see Table 1), these results contrast with the reduced polarization of the Zn–N coordination inferred from the partial charge distribution (Figure 2). In an earlier DFT study on ZIF-8 complexes, similar values of the Zn–N distances were reported. (27) This agreement should exclude artifacts from our calculations and, concomitantly, suggests that in periodic systems the long-range potential may contribute to enhancing the polar character of the Zn–N bond.
Finally, we inspect the partial charge distribution in the imidazole ligands with their respective terminations. The model systems adopted for this analysis are taken from the crystal structures and optimized in vacuo. While this analysis evidently does not give us any information about the Zn–N bond, since the metal ion is not a part of the molecules, it enables us to better assess the charge distribution within the conjugated network and the functional group. In the absence of Zn, the excess of negative charge on N decreases in magnitude by a few hundredths of e in all systems (see Figure 2). The most significant reduction affects the Cl-functionalized compound, where the electron-withdrawing ability of the halogen atom enhances the ionicity of the whole system. The charge distribution within the conjugated network, consisting in this case only of three C and three H atoms (see Figure 1c), is almost identical to the one obtained in the complexes. From this analysis, we conclude that the presence of Zn in the complex impacts merely the local charge distribution on N.
For a better overview of the charge distribution in the considered systems, we complement the partial charge analysis presented above with the visualization of the electron deformation density (EDD) plotted along the molecular plane of the linker molecule (Figure 3). The EDD is defined as the difference between the electron density of the considered structure and the so-called promolecular density, which is a superposition of the densities of the free atoms positioned at the molecular/crystalline sites within the structure: (55) positive values (red domains in Figure 3) represent an accumulation of electrons while negative values (blue domains in Figure 3) represent an electron depletion with respect to the free atoms. The extension of the former over the carbon-based network can be associated with the presence of covalent bonds therein, which seems to be independent of the state of matter (crystals or complexes). In the presence of CH3 and H terminations, the functional groups are covalently bound to the backbone (Figure 3, left panels). In contrast, the halogen atoms are connected to the inner network by a highly polar bond, as demonstrated by the blueish regions between them (Figure 3, right panels). Likewise, the ionicity of the Zn–N bond is visible in all plots in the middle and bottom panels.

Figure 3

Figure 3. Electron deformation density plotted along the molecular plane of a single molecule terminated by CH3, H, Cl, and Br, in the three considered configurations: isolated imidazole (top), complex (middle), and crystal (bottom). Electron depletion is indicated in blue and electron accumulation in red.

X-ray Absorption Spectroscopy

We complement the analysis of the charge-density distribution by inspecting the results of X-ray absorption measurements performed from the Zn K-edge. These experiments provide an indication, at least on a qualitative level, of the Zn coordination in the ZIF-8 environment in the presence of different ligand terminations. To this end, measurements were performed on a sample of zinc(II)2-bromoimidazolate (ZIF-8-Br) and zinc(II)2-methylimidazolate (ZIF-8) as a reference (Figure 4).

Figure 4

Figure 4. (a) X-ray absorption from the Zn K-edge and (b) its first derivative measured in zinc(II)2-bromoimidazolate and zinc(II)2-methylimidazolate.

Figure 4a displays XANES spectra at the Zn K-edge of zinc(II)2-methylimidazolate (blue) and zinc(II)2-bromoimidazolate (orange) measured in transmission mode. The spectra reveal only subtle differences in the position of both absorption maxima, which are characteristic of (functionalized) ZIF-8 structures. (64) The first maximum at 9664 eV is due to transitions from Zn-1s → Zn-4sp + ligand-2p states, while the second maximum at 9672 eV stems from Zn-1s → Zn-4sp + ligand-2p orbitals. This characteristic becomes more evident when looking at the first derivative of the signal plotted in Figure 4b, where the spectrum of the former compound exhibits a minimum at about 0.5 eV below in energy compared to its counterpart in the latter. This result can be interpreted considering that a higher-energy resonance at the onset of absorption from the Zn K-edge is compatible with a lower electron density in this metal atom. Hence, according to this line of reasoning, the partial charge on Zn is more positive in zinc(II)2-bromoimidazolate than in zinc(II)2-methylimidazolate. This is consistent with the picture provided by the DFT results.

Electronic Properties

We conclude our study by turning to the electronic structures of the considered systems. To this end, we focus only on the periodic crystals and assess the effects of the different ligand terminations on the character and energy of the electronic states in the vicinity of the frontier. We inspect the projected density of states (PDOS) of the considered crystal structures including the atom-resolved orbital components (see Figure 5; in Figure S1, the corresponding band structures are reported). Before proceeding with this discussion, it is worth mentioning that these results are obtained from DFT using the PBE functional, which notoriously underestimates band gaps. The corresponding values in Figure 5 are therefore affected by this methodological shortcoming and should not be interpreted quantitatively. However, this issue does not undermine the following analysis of the qualitative features, which are instead well reproduced within the adopted approximations.

Figure 5

Figure 5. PDOS of (a) ZIF-8 and its counterparts with (b) H, (c) Cl, and (d) Br functionalization. The shaded gray area indicates the total density of states (TDOS) summing up all atomic (color code in the legend) contributions. The Fermi energy is set to 0 eV at the top of the valence band and is marked by a vertical dashed bar.

In all systems, the top of the valence band is dominated by a sharp peak, corresponding to a localized band manifold (see Figure S1 and Table S9 for the related band gap analysis), mainly including contributions from the C atoms. Deeper valence states form instead a broader band with prevailing N contributions. In ZIF-8 with halogen terminations, Cl and Br p-states are present in this manifold, too (see Figure 5c,d). This result is unsurprising considering the electronic structure of these atomic species with an almost filled p-shell at the top of their occupied region. At even lower energies, below −4 eV, we find the contributions of Zn d-orbitals hybridized with C and N states within the inner network. By comparing the energy distribution of the two valence-band manifolds between −3 and −1 eV, and below −4 eV, we can assess the influence of the various ligand terminations. Halogen atoms enhance the energy separation of the N-dominated region shifting it down by about 0.5 eV compared to ZIF-8 and ZIF-8-H (see Figure 5). Moreover, the contributions of the Cl and Br atoms in the same range extend energetically this band down to about −3 eV (ZIF-8-Br, Figure 5d) to even below that threshold (ZIF-8-Cl, Figure 5c). Notice, for comparison, that the N-dominated valence-band manifold extends only down to −2.5 eV in both ZIF-8 and ZIF-8-H (Figure 5a,b).
Moving now to the conduction region, we notice at a glance that the Br-termination leads to a reduction of the band gap by more than 0.5 eV compared to ZIF-8 (Figure 5d,a). The bottom of the unoccupied region in ZIF-8-Br is indeed formed by a Br–C hybridized state. In ZIF-8-Cl, halogen contributions to the conduction bands are visible above the band onset (Figure 5c). Likewise, in both ZIF-8 and ZIF-8-H, the relevant contributions of the functional groups to the conduction region are above the energy ranges visualized in Figure 5a,b. In these systems, the lowest-unoccupied band manifold is dominated by hybridized C and N states which are present also in the halogen-terminated structures except that in ZIF-8-Br, these states appear above two unoccupied Br-based states at the conduction-band bottom and are thus up-shifted by approximately 1 eV compared to the other systems.
The analysis of the PDOS suggests that the main role of the different ligand terminations is to alter the size of the band gap and the relative energies of the band manifold in both the valence and conduction regions. These effects are mostly pronounced with halogen functionalizations due to their electronegativity, which is known to affect the electronic properties of the carbon-based backbone hosting them through the insertion of states in the gap region (as for ZIF-8-Br here) or by downshifting the electronic states as a whole. A similar trend was found and rationalized in the context of carbon-based nanostructures. (65)

Conclusions

Click to copy section linkSection link copied!

In summary, we presented a joint first-principles and X-ray spectroscopic analysis of the influence of ligand substitution on the character of Zn-coordination in ZIFs. We inspected the established ZIF-8 compound and its derivatives with halogen (Cl and Br) as well as H-substitution in the organic ligand. By performing DFT calculations both on the crystalline compounds and their molecular building blocks as extracted from the bulk, we assessed the influence of long-range interactions in the local environment of the Zn–N coordination. Our results indicate that the electronegativity of the ligand substituents impacts the polarization of the Zn–N bond. The trends suggested by the DFT calculations are confirmed by the outcomes of X-ray absorption measurements from the Zn K-edge on ZIF-8 and its Br-functionalized counterpart. We completed our study by inspecting the electronic properties of the periodic MOFs considered in this work. We found that the influence of the ligand substituents touches the size of the band gap, with a reduction on the order of 0.5 eV upon halogen functionalization. Furthermore, due to their high electronegativity, Cl and Br atoms attached to the linkers lead to an overall downshift of the electronic bands in comparison to ZIF-8 or its hydrogen-terminated variant.
To conclude, this study offers new insights into the nature of Zn coordination in ZIF-8 and related compounds and its effects on the charge-density distribution, including the polarization of the Zn–N bond, as well as on the electronic structure of these systems. As such, our findings confirm that chemical manipulation, for example, through ligand functionalization, represents a viable route to tune their physical properties in a flexible and controlled way.

Supporting Information

Click to copy section linkSection link copied!

The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.jpcc.3c06054.

  • Lattice parameters of all the investigated structures; average bond lengths of all the considered structures before and after optimization; partial charges and band gaps, including the energies of the frontier levels; band structure plots (PDF)

Terms & Conditions

Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.

Author Information

Click to copy section linkSection link copied!

  • Corresponding Authors
    • Sebastian Beyer - Federal Institute for Materials Research and Testing (BAM), Berlin 12485, GermanyDepartment of Biomedical Engineering, Hong Kong Special Administrative Region of China, The Chinese University of Hong Kong, Hong Kong 999077, ChinaOrcidhttps://orcid.org/0000-0001-7945-1477 Email: [email protected]
    • Caterina Cocchi - Carl von Ossietzky Universität Oldenburg, Institute of Physics, Oldenburg 26129, GermanyCenter for Nanoscale Dynamics, Carl von Ossietzky Universität Oldenburg, Oldenburg 26129, GermanyPhysics Department and IRIS Adlershof, Humboldt-Universität zu Berlin, Berlin 12489, GermanyOrcidhttps://orcid.org/0000-0002-9243-9461 Email: [email protected]
  • Authors
    • Joshua Edzards - Carl von Ossietzky Universität Oldenburg, Institute of Physics, Oldenburg 26129, Germany
    • Holger-Dietrich Saßnick - Carl von Ossietzky Universität Oldenburg, Institute of Physics, Oldenburg 26129, GermanyOrcidhttps://orcid.org/0000-0001-8887-9248
    • Ana Guilherme Buzanich - Federal Institute for Materials Research and Testing (BAM), Berlin 12485, GermanyOrcidhttps://orcid.org/0000-0001-5543-9924
    • Ana M. Valencia - Carl von Ossietzky Universität Oldenburg, Institute of Physics, Oldenburg 26129, GermanyPhysics Department and IRIS Adlershof, Humboldt-Universität zu Berlin, Berlin 12489, GermanyOrcidhttps://orcid.org/0000-0003-0095-3680
    • Franziska Emmerling - Federal Institute for Materials Research and Testing (BAM), Berlin 12485, GermanyOrcidhttps://orcid.org/0000-0001-8528-0301
  • Notes
    The authors declare no competing financial interest.

Acknowledgments

Click to copy section linkSection link copied!

J.E., H.-D.S., A.M.V., and C.C. acknowledge funding from the German Federal Ministry of Education and Research (Professorinnenprogramm III) as well as from the State of Lower Saxony (Professorinnen für Niedersachsen and DyNano). S.B. is supported by the Faculty of Engineering of The Chinese University of Hong Kong (CUHK), grant number 4055120. Computational resources were provided by the North-German Supercomputing Alliance (HLRN), project nic00069, and by the local high-performance computing cluster CARL at the University of Oldenburg, financed by the German Research Foundation (project no. INST 184/157-1 FUGG) and by the Ministry of Science and Culture of the State of Lower Saxony. Experiments were performed at the BAMline at the BESSY-II storage ring (Helmholtz Center Berlin). We thank the Helmholtz-Zentrum Berlin für Materialien und Energie for the allocation of synchrotron radiation beamtime.

References

Click to copy section linkSection link copied!

This article references 65 other publications.

  1. 1
    Park, K. S.; Ni, Z.; Côté, A. P.; Choi, J. Y.; Huang, R.; Uribe-Romo, F. J.; Chae, H. K.; O’Keeffe, M.; Yaghi, O. M. Exceptional chemical and thermal stability of zeolitic imidazolate frameworks. Proc. Natl. Acad. Sci. U.S.A. 2006, 103, 1018610191,  DOI: 10.1073/pnas.0602439103
  2. 2
    Huang, X.-C.; Lin, Y.-Y.; Zhang, J.-P.; Chen, X.-M. Ligand-directed strategy for zeolite-type metal–organic frameworks: zinc (II) imidazolates with unusual zeolitic topologies. Angew. Chem., Int. Ed. 2006, 45, 15571559,  DOI: 10.1002/anie.200503778
  3. 3
    Bauman, J. E., Jr; Wang, J. C. Imidazole complexes of nickel (II), copper (II), zinc (II), and silver (I). Inorg. Chem. 1964, 3, 368373,  DOI: 10.1021/ic50013a014
  4. 4
    Yaghi, O. M.; Li, G.; Li, H. Selective binding and removal of guests in a microporous metal–organic framework. Nature 1995, 378, 703706,  DOI: 10.1038/378703a0
  5. 5
    Li, H.; Eddaoudi, M.; O’Keeffe, M.; Yaghi, O. M. Design and synthesis of an exceptionally stable and highly porous metal-organic framework. Nature 1999, 402, 276279,  DOI: 10.1038/46248
  6. 6
    Beyer, S.; Prinz, C.; Schürmann, R.; Feldmann, I.; Zimathies, A.; Blocki, A. M.; Bald, I.; Schneider, R. J.; Emmerling, F. Ultra-sonication of ZIF-67 crystals results in ZIF-67 nano-flakes. ChemistrySelect 2016, 1, 59055908,  DOI: 10.1002/slct.201601513
  7. 7
    Beyer, S.; Schürmann, R.; Feldmann, I.; Blocki, A.; Bald, I.; Schneider, R. J.; Emmerling, F. Maintaining stable Zeolitic Imidazolate Framework (ZIF) templates during polyelectrolyte multilayer coating. Colloids Interface Sci. Ser. 2018, 22, 1417,  DOI: 10.1016/j.colcom.2017.11.004
  8. 8
    Buzanich, A. G.; Kulow, A.; Kabelitz, A.; Grunewald, C.; Seidel, R.; Chapartegui-Arias, A.; Radtke, M.; Reinholz, U.; Emmerling, F.; Beyer, S. Observation of early ZIF-8 crystallization stages with X-ray absorption spectroscopy. Soft Matter 2021, 17, 331334,  DOI: 10.1039/D0SM01356K
  9. 9
    Phan, A.; Doonan, C. J.; Uribe-Romo, F. J.; Knobler, C. B.; O’Keeffe, M.; Yaghi, O. M. Synthesis, Structure, and Carbon Dioxide Capture Properties of Zeolitic Imidazolate Frameworks. Acc. Chem. Res. 2010, 43, 5867,  DOI: 10.1021/ar900116g
  10. 10
    Tan, J. C.; Bennett, T. D.; Cheetham, A. K. Chemical structure, network topology, and porosity effects on the mechanical properties of Zeolitic Imidazolate Frameworks. Proc. Natl. Acad. Sci. U.S.A. 2010, 107, 99389943,  DOI: 10.1073/pnas.1003205107
  11. 11
    Chen, B.; Yang, Z.; Zhu, Y.; Xia, Y. Zeolitic imidazolate framework materials: recent progress in synthesis and applications. J. Mater. Chem. A 2014, 2, 1681116831,  DOI: 10.1039/C4TA02984D
  12. 12
    Spencer, E. C.; Angel, R. J.; Ross, N. L.; Hanson, B. E.; Howard, J. A. K. Pressure-Induced Cooperative Bond Rearrangement in a Zinc Imidazolate Framework: A High-Pressure Single-Crystal X-Ray Diffraction Study. J. Am. Chem. Soc. 2009, 131, 40224026,  DOI: 10.1021/ja808531m
  13. 13
    Hu, Y.; Kazemian, H.; Rohani, S.; Huang, Y.; Song, Y. In situ high pressure study of ZIF-8 by FTIR spectroscopy. Chem. Commun. 2011, 47, 1269412696,  DOI: 10.1039/c1cc15525c
  14. 14
    Hu, Y.; Liu, Z.; Xu, J.; Huang, Y.; Song, Y. Evidence of pressure enhanced CO2 storage in ZIF-8 probed by FTIR spectroscopy. J. Am. Chem. Soc. 2013, 135, 92879290,  DOI: 10.1021/ja403635b
  15. 15
    Widmer, R. N.; Lampronti, G. I.; Anzellini, S.; Gaillac, R.; Farsang, S.; Zhou, C.; Belenguer, A. M.; Wilson, C. W.; Palmer, H.; Kleppe, A. K. Pressure promoted low-temperature melting of metal–organic frameworks. Nat. Mater. 2019, 18, 370376,  DOI: 10.1038/s41563-019-0317-4
  16. 16
    Choi, J.; Im, J.; Noh, K.; Kim, J.; Vogt, T.; Lee, Y. Universal Gas-Uptake Behavior of a Zeolitic Imidazolate Framework ZIF-8 at High Pressure. J. Phys. Chem. C 2019, 123, 2576925774,  DOI: 10.1021/acs.jpcc.9b08539
  17. 17
    Formalik, F.; Mazur, B.; Fischer, M.; Firlej, L.; Kuchta, B. Phonons and Adsorption-Induced Deformations in ZIFs: Is It Really a Gate Opening?. J. Phys. Chem. C 2021, 125, 79998005,  DOI: 10.1021/acs.jpcc.1c01342
  18. 18
    Iacomi, P.; Maurin, G. ResponZIF Structures: Zeolitic Imidazolate Frameworks as Stimuli-Responsive Materials. ACS Appl. Mater. Interfaces 2021, 13, 5060250642,  DOI: 10.1021/acsami.1c12403
  19. 19
    Eslava, S.; Zhang, L.; Esconjauregui, S.; Yang, J.; Vanstreels, K.; Baklanov, M. R.; Saiz, E. Metal-organic framework ZIF-8 films as low-κ dielectrics in microelectronics. Chem. Mater. 2013, 25, 2733,  DOI: 10.1021/cm302610z
  20. 20
    Pimentel, B. R.; Parulkar, A.; Zhou, E.-k.; Brunelli, N. A.; Lively, R. P. Zeolitic imidazolate frameworks: next-generation materials for energy-efficient gas separations. ChemSusChem 2014, 7, 32023240,  DOI: 10.1002/cssc.201402647
  21. 21
    Hoop, M.; Walde, C. F.; Riccò, R.; Mushtaq, F.; Terzopoulou, A.; Chen, X.-Z.; deMello, A. J.; Doonan, C. J.; Falcaro, P.; Nelson, B. J. Biocompatibility characteristics of the metal organic framework ZIF-8 for therapeutical applications. Appl. Mater. Today 2018, 11, 1321,  DOI: 10.1016/j.apmt.2017.12.014
  22. 22
    Dai, H.; Yuan, X.; Jiang, L.; Wang, H.; Zhang, J.; Zhang, J.; Xiong, T. Recent advances on ZIF-8 composites for adsorption and photocatalytic wastewater pollutant removal: Fabrication, applications and perspective. Coord. Chem. Rev. 2021, 441, 213985,  DOI: 10.1016/j.ccr.2021.213985
  23. 23
    Hu, W.; Pattengale, B.; Huang, J. Zeolitic imidazolate frameworks as intrinsic light harvesting and charge separation materials for photocatalysis. J. Chem. Phys. 2021, 154, 240901,  DOI: 10.1063/5.0048720
  24. 24
    Knebel, A.; Caro, J. Metal–organic frameworks and covalent organic frameworks as disruptive membrane materials for energy-efficient gas separation. Nat. Nanotechnol. 2022, 17, 911923,  DOI: 10.1038/s41565-022-01168-3
  25. 25
    Paul, A.; Banga, I. K.; Muthukumar, S.; Prasad, S. Engineering the ZIF-8 Pore for Electrochemical Sensor Applications - A Mini Review. ACS Omega 2022, 7, 2699327003,  DOI: 10.1021/acsomega.2c00737
  26. 26
    Tsang, C. Y.; Cheung, M. C. Y.; Beyer, S. Assessing the colloidal stability of copper doped ZIF-8 in water and serum. Colloids Surf., A 2023, 656, 130452,  DOI: 10.1016/j.colsurfa.2022.130452
  27. 27
    Wang, H.; Zhao, L.; Xu, W.; Wang, S.; Ding, Q.; Lu, X.; Guo, W. The properties of the bonding between CO and ZIF-8 structures: a density functional theory study. Theor. Chim. Acta 2015, 134, 3139,  DOI: 10.1007/s00214-015-1636-4
  28. 28
    Yilmaz, G.; Peh, S. B.; Zhao, D.; Ho, G. W. Atomic-and Molecular-Level Design of Functional Metal–Organic Frameworks (MOFs) and Derivatives for Energy and Environmental Applications. Adv. Sci. 2019, 6, 1901129,  DOI: 10.1002/advs.201901129
  29. 29
    Liu, Y.; Huo, Y.; Wang, X.; Yu, S.; Ai, Y.; Chen, Z.; Zhang, P.; Chen, L.; Song, G.; Alharbi, N. S. Impact of metal ions and organic ligands on uranium removal properties by zeolitic imidazolate framework materials. J. Clean. Prod. 2021, 278, 123216,  DOI: 10.1016/j.jclepro.2020.123216
  30. 30
    Wang, H.; Pei, X.; Kalmutzki, M. J.; Yang, J.; Yaghi, O. M. Large Cages of Zeolitic Imidazolate Frameworks. Acc. Chem. Res. 2022, 55, 707721,  DOI: 10.1021/acs.accounts.1c00740
  31. 31
    Bhattacharyya, S.; Han, R.; Kim, W.-G.; Chiang, Y.; Jayachandrababu, K. C.; Hungerford, J. T.; Dutzer, M. R.; Ma, C.; Walton, K. S.; Sholl, D. S. Acid Gas Stability of Zeolitic Imidazolate Frameworks: Generalized Kinetic and Thermodynamic Characteristics. Chem. Mater. 2018, 30, 40894101,  DOI: 10.1021/acs.chemmater.8b01394
  32. 32
    Sarkar, S.; Grønbech, T. B. E.; Mamakhel, A.; Bondesgaard, M.; Sugimoto, K.; Nishibori, E.; Iversen, B. B. X-ray Electron Density Study of the Chemical Bonding Origin of Glass Formation in Metal–Organic Frameworks. Angew. Chem., Int. Ed. 2022, 61, e202202742  DOI: 10.1002/anie.202202742
  33. 33
    Butler, K. T.; Worrall, S. D.; Molloy, C. D.; Hendon, C. H.; Attfield, M. P.; Dryfe, R. A.; Walsh, A. Electronic structure design for nanoporous, electrically conductive zeolitic imidazolate frameworks. J. Mater. Chem. C 2017, 5, 77267731,  DOI: 10.1039/C7TC03150E
  34. 34
    Möslein, A. F.; Tan, J.-C. Vibrational Modes and Terahertz Phenomena of the Large-Cage Zeolitic Imidazolate Framework-71. J. Phys. Chem. Lett. 2022, 13, 28382844,  DOI: 10.1021/acs.jpclett.2c00081
  35. 35
    Möslein, A. F.; Donà, L.; Civalleri, B.; Tan, J.-C. Defect Engineering in Metal–Organic Framework Nanocrystals: Implications for Mechanical Properties and Performance. ACS Appl. Mater. Interfaces 2022, 5, 63986409,  DOI: 10.1021/acsanm.2c00493
  36. 36
    Graiver, D.; Farminer, K.; Narayan, R. A review of the fate and effects of silicones in the environment. J. Polym. Environ. 2003, 11, 129136,  DOI: 10.1023/A:1026056129717
  37. 37
    Du, X.; Zhou, M. Strategies to enhance catalytic performance of metal-organic frameworks in sulfate radical-based advanced oxidation processes for organic pollutants removal. Chem. Eng. J. 2021, 403, 126346,  DOI: 10.1016/j.cej.2020.126346
  38. 38
    Lewis, D. W.; Ruiz-Salvador, A. R.; Gómez, A.; Rodriguez-Albelo, L. M.; Coudert, F.-X.; Slater, B.; Cheetham, A. K.; Mellot-Draznieks, C. Zeolitic imidazole frameworks: structural and energetics trends compared with their zeolite analogues. CrystEngComm 2009, 11, 22722276,  DOI: 10.1039/b912997a
  39. 39
    Li, K.; Olson, D. H.; Seidel, J.; Emge, T. J.; Gong, H.; Zeng, H.; Li, J. Zeolitic Imidazolate Frameworks for Kinetic Separation of Propane and Propene. J. Am. Chem. Soc. 2009, 131, 1036810369,  DOI: 10.1021/ja9039983
  40. 40
    Amrouche, H.; Aguado, S.; Pérez-Pellitero, J.; Chizallet, C.; Siperstein, F.; Farrusseng, D.; Bats, N.; Nieto-Draghi, C. Experimental and Computational Study of Functionality Impact on Sodalite–Zeolitic Imidazolate Frameworks for CO2 Separation. J. Phys. Chem. C 2011, 115, 1642516432,  DOI: 10.1021/jp202804g
  41. 41
    Chaplais, G.; Fraux, G.; Paillaud, J.-L.; Marichal, C.; Nouali, H.; Fuchs, A. H.; Coudert, F.-X.; Patarin, J. Impacts of the Imidazolate Linker Substitution (CH3, Cl, or Br) on the Structural and Adsorptive Properties of ZIF-8. J. Phys. Chem. C 2018, 122, 2694526955,  DOI: 10.1021/acs.jpcc.8b08706
  42. 42
    Yagi, R.; Ueda, T. Substitution (CH 3, Cl, or Br) effects of the imidazolate linker on benzene adsorption kinetics for the zeolitic imidazolate framework (ZIF)-8. Phys. Chem. Chem. Phys. 2023, 25, 2058520596,  DOI: 10.1039/D3CP01662E
  43. 43
    Amrouche, H.; Creton, B.; Siperstein, F.; Nieto-Draghi, C. Prediction of thermodynamic properties of adsorbed gases in zeolitic imidazolate frameworks. RSC Adv. 2012, 2, 60286035,  DOI: 10.1039/c2ra00025c
  44. 44
    Dürholt, J. P.; Fraux, G.; Coudert, F.-X.; Schmid, R. Ab Initio Derived Force Fields for Zeolitic Imidazolate Frameworks: MOF-FF for ZIFs. J. Chem. Theory Comput. 2019, 15, 24202432,  DOI: 10.1021/acs.jctc.8b01041
  45. 45
    Morris, W.; Stevens, C. J.; Taylor, R. E.; Dybowski, C.; Yaghi, O. M.; Garcia-Garibay, M. A. NMR and X-ray Study Revealing the Rigidity of Zeolitic Imidazolate Frameworks. J. Phys. Chem. C 2012, 116, 1330713312,  DOI: 10.1021/jp303907p
  46. 46
    Ke, Q.; Duan, Y.; Ji, Y.; Zhao, D.; Zhang, H.; Duan, C.; Li, L.; Wei, Y. Identical Composition and Distinct Performance: How ZIF-8 Polymorphs’ Structures Affect the Adsorption/Separation of Ethane and Ethene. J. Chem. Eng. Data 2021, 66, 34833492,  DOI: 10.1021/acs.jced.1c00322
  47. 47
    Hinuma, Y.; Pizzi, G.; Kumagai, Y.; Oba, F.; Tanaka, I. Band structure diagram paths based on crystallography. Comput. Mater. Sci. 2017, 128, 140184,  DOI: 10.1016/j.commatsci.2016.10.015
  48. 48
    Hohenberg, P.; Kohn, W. Inhomogeneous Electron Gas. Phys. Rev. 1964, 136, B864B871,  DOI: 10.1103/PhysRev.136.B864
  49. 49
    Kohn, W.; Sham, L. J. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A1133A1138,  DOI: 10.1103/PhysRev.140.A1133
  50. 50
    Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 38653868,  DOI: 10.1103/PhysRevLett.77.3865
  51. 51
    Giannozzi, P.; Baseggio, O.; Bonfà, P.; Brunato, D.; Car, R.; Carnimeo, I.; Cavazzoni, C.; de Gironcoli, S.; Delugas, P.; Ferrari Ruffino, F. Quantum ESPRESSO toward the exascale. J. Chem. Phys. 2020, 152, 154105,  DOI: 10.1063/5.0005082
  52. 52
    Blöchl, P. E. Projector augmented-wave method. Phys. Rev. B 1994, 50, 1795317979,  DOI: 10.1103/PhysRevB.50.17953
  53. 53
    Dal Corso, A. Pseudopotentials periodic table: From H to Pu. Comput. Mater. Sci. 2014, 95, 337350,  DOI: 10.1016/j.commatsci.2014.07.043
  54. 54
    Bader, R.; Zou, P. An atomic population as the expectation value of a quantum observable. Chem. Phys. Lett. 1992, 191, 5458,  DOI: 10.1016/0009-2614(92)85367-J
  55. 55
    Otero-de-la-Roza, A.; Johnson, E. R.; Luaña, V. Critic2: A program for real-space analysis of quantum chemical interactions in solids. Comput. Phys. Commun. 2014, 185, 10071018,  DOI: 10.1016/j.cpc.2013.10.026
  56. 56
    Yu, M.; Trinkle, D. R. Accurate and efficient algorithm for Bader charge integration. J. Chem. Phys. 2011, 134, 064111,  DOI: 10.1063/1.3553716
  57. 57
    Guilherme Buzanich, A.; Radtke, M.; Yusenko, K. V.; Stawski, T. M.; Kulow, A.; Cakir, C. T.; Röder, B.; Naese, C.; Britzke, R.; Sintschuk, M. BAMline─A real-life sample materials research beamline. J. Chem. Phys. 2023, 158, 244202,  DOI: 10.1063/5.0157194
  58. 58
    Ravel, B.; Newville, M. ATHENA, ARTEMIS, HEPHAESTUS: data analysis for X-ray absorption spectroscopy using IFEFFIT. J. Synchrotron Radiat. 2005, 12, 537541,  DOI: 10.1107/S0909049505012719
  59. 59
    Cui, K.; Schmidt, J. Enabling efficient and accurate computational studies of MOF reactivity via QM/MM and QM/QM methods. J. Phys. Chem. C 2020, 124, 1055010560,  DOI: 10.1021/acs.jpcc.0c01220
  60. 60
    Hartmann, M.; Clark, T.; van Eldik, R. Hydration and water exchange of zinc (II) ions. Application of density functional theory. J. Am. Chem. Soc. 1997, 119, 78437850,  DOI: 10.1021/ja970483f
  61. 61
    Bock, C. W.; Katz, A. K.; Markham, G. D.; Glusker, J. P. Manganese as a replacement for magnesium and zinc: functional comparison of the divalent ions. J. Am. Chem. Soc. 1999, 121, 73607372,  DOI: 10.1021/ja9906960
  62. 62
    Sandmark, C.; Brändén, C. I.; Kjekshus, A.; Åkeson, Å.; Theorell, H.; Blinc, R.; Paušak, S.; Ehrenberg, L.; Dumanović, J. The Crystal Structure of Hexaimidazole Zinc (II) Dichloride Tetrahydrate, Zn(C3H4N2)6Cl2·4H2O. Acta Chem. Scand. 1967, 21, 993999,  DOI: 10.3891/acta.chem.scand.21-0993
  63. 63
    Torzilli, M. A.; Colquhoun, S.; Doucet, D.; Beer, R. H. The interconversion of dichlorobis (Nn-propylsalicylaldimine) zinc (II) and bis (Nn-propylsalicylaldiminato) zinc (II). Polyhedron 2002, 21, 697704,  DOI: 10.1016/S0277-5387(02)00839-2
  64. 64
    Boada, R.; Chaboy, J.; Hayama, S.; Keenan, L. L.; Freeman, A. A.; Amboage, M.; Díaz-Moreno, S. Unraveling the Molecular Details of the “Gate Opening” Phenomenon in Zif-8 with X-Ray Absorption Spectroscopy. J. Phys. Chem. C 2022, 126, 59355943,  DOI: 10.1021/acs.jpcc.2c00373
  65. 65
    Cocchi, C.; Prezzi, D.; Ruini, A.; Caldas, M. J.; Molinari, E. Optical properties and charge-transfer excitations in edge-functionalized all-graphene nanojunctions. J. Phys. Chem. Lett. 2011, 2, 13151319,  DOI: 10.1021/jz200472a

Cited By

Click to copy section linkSection link copied!
Citation Statements
Explore this article's citation statements on scite.ai

This article is cited by 7 publications.

  1. Mueed Ahmad, Roshan Patel, Dennis T. Lee, Peter Corkery, Andrea Kraetz, Prerna, Samuel A. Tenney, Dmytro Nykypanchuk, Xiao Tong, J. Ilja Siepmann, Michael Tsapatsis, J. Anibal Boscoboinik. ZIF-8 Vibrational Spectra: Peak Assignments and Defect Signals. ACS Applied Materials & Interfaces 2024, 16 (21) , 27887-27897. https://doi.org/10.1021/acsami.4c02396
  2. Holger-Dietrich Saßnick, Fabiana Machado Ferreira De Araujo, Joshua Edzards, Caterina Cocchi. Impact of Ligand Substitution and Metal Node Exchange in the Electronic Properties of Scandium Terephthalate Frameworks. Inorganic Chemistry 2024, 63 (4) , 2098-2108. https://doi.org/10.1021/acs.inorgchem.3c03945
  3. Kevin Dedecker, Martin Drobek, Anne Julbe. Selective adsorption and separation of C 6 hydrocarbons: the role of structural flexibility and functionalization in zeolitic imidazolate frameworks. RSC Applied Interfaces 2025, 2 (2) , 364-372. https://doi.org/10.1039/D4LF00388H
  4. Sihem Bouallag, A. Mougari, M. Zabat, A. Belayadi. Computational study of transition metal coordinated polyaniline: A first principle investigation into tuning the electronic properties of the resulting hybrid material. Physica B: Condensed Matter 2024, 695 , 416597. https://doi.org/10.1016/j.physb.2024.416597
  5. Joshua Edzards, Holger-Dietrich Saßnick, Julia Santana Andreo, Caterina Cocchi. Tuning structural and electronic properties of metal-organic framework 5 by metal substitution and linker functionalization. The Journal of Chemical Physics 2024, 160 (18) https://doi.org/10.1063/5.0203605
  6. Miao Miao Zhao, Yang Yang, Na Du, Yu Ying Zhu, Peng Ren, Fei Yen. Determining the chemical composition of diamagnetic mixed solids via measurements of their magnetic susceptibility. Journal of Materials Chemistry C 2024, 12 (16) , 5877-5885. https://doi.org/10.1039/D4TC00757C
  7. Satyendra Verma, Rinki Brahma, Jubaraj B. Baruah. Consequences of twisting of the flexible arms of imidazole-derived urea in zinc-dicarboxylate coordination polymers. CrystEngComm 2024, 26 (14) , 1976-1985. https://doi.org/10.1039/D4CE00005F

The Journal of Physical Chemistry C

Cite this: J. Phys. Chem. C 2023, 127, 43, 21456–21464
Click to copy citationCitation copied!
https://doi.org/10.1021/acs.jpcc.3c06054
Published October 25, 2023

Copyright © 2023 The Authors. Published by American Chemical Society. This publication is licensed under

CC-BY 4.0 .

Article Views

1304

Altmetric

-

Citations

Learn about these metrics

Article Views are the COUNTER-compliant sum of full text article downloads since November 2008 (both PDF and HTML) across all institutions and individuals. These metrics are regularly updated to reflect usage leading up to the last few days.

Citations are the number of other articles citing this article, calculated by Crossref and updated daily. Find more information about Crossref citation counts.

The Altmetric Attention Score is a quantitative measure of the attention that a research article has received online. Clicking on the donut icon will load a page at altmetric.com with additional details about the score and the social media presence for the given article. Find more information on the Altmetric Attention Score and how the score is calculated.

  • Abstract

    Figure 1

    Figure 1. Ball-and-stick representation of the ZIF-8 compounds considered in this work including the (a) crystal structure with R3 space-group symmetry, (b) isolated metal–organic complex replicated in the crystal, and (c) imidazole ligand. Carbon atoms are depicted in black, hydrogen atoms in white, nitrogen atoms in blue, zinc atoms in gray, and the terminating group in green (CH3 in ZIF-8, H, Cl, and Br in the substituted siblings). The atoms in the ligand are marked in panel c: identical labels indicate equivalent atoms; Ha is the proton added to the imidazole unit and to the complex to passivate one of the two N atoms.

    Figure 2

    Figure 2. Partial charges per atom in the considered ZIF structures. Results reported for the “additional H” concern only the complexes and isolated imidazole units, where protonation of a dangling bond was needed to neutralize the system.

    Figure 3

    Figure 3. Electron deformation density plotted along the molecular plane of a single molecule terminated by CH3, H, Cl, and Br, in the three considered configurations: isolated imidazole (top), complex (middle), and crystal (bottom). Electron depletion is indicated in blue and electron accumulation in red.

    Figure 4

    Figure 4. (a) X-ray absorption from the Zn K-edge and (b) its first derivative measured in zinc(II)2-bromoimidazolate and zinc(II)2-methylimidazolate.

    Figure 5

    Figure 5. PDOS of (a) ZIF-8 and its counterparts with (b) H, (c) Cl, and (d) Br functionalization. The shaded gray area indicates the total density of states (TDOS) summing up all atomic (color code in the legend) contributions. The Fermi energy is set to 0 eV at the top of the valence band and is marked by a vertical dashed bar.

  • References


    This article references 65 other publications.

    1. 1
      Park, K. S.; Ni, Z.; Côté, A. P.; Choi, J. Y.; Huang, R.; Uribe-Romo, F. J.; Chae, H. K.; O’Keeffe, M.; Yaghi, O. M. Exceptional chemical and thermal stability of zeolitic imidazolate frameworks. Proc. Natl. Acad. Sci. U.S.A. 2006, 103, 1018610191,  DOI: 10.1073/pnas.0602439103
    2. 2
      Huang, X.-C.; Lin, Y.-Y.; Zhang, J.-P.; Chen, X.-M. Ligand-directed strategy for zeolite-type metal–organic frameworks: zinc (II) imidazolates with unusual zeolitic topologies. Angew. Chem., Int. Ed. 2006, 45, 15571559,  DOI: 10.1002/anie.200503778
    3. 3
      Bauman, J. E., Jr; Wang, J. C. Imidazole complexes of nickel (II), copper (II), zinc (II), and silver (I). Inorg. Chem. 1964, 3, 368373,  DOI: 10.1021/ic50013a014
    4. 4
      Yaghi, O. M.; Li, G.; Li, H. Selective binding and removal of guests in a microporous metal–organic framework. Nature 1995, 378, 703706,  DOI: 10.1038/378703a0
    5. 5
      Li, H.; Eddaoudi, M.; O’Keeffe, M.; Yaghi, O. M. Design and synthesis of an exceptionally stable and highly porous metal-organic framework. Nature 1999, 402, 276279,  DOI: 10.1038/46248
    6. 6
      Beyer, S.; Prinz, C.; Schürmann, R.; Feldmann, I.; Zimathies, A.; Blocki, A. M.; Bald, I.; Schneider, R. J.; Emmerling, F. Ultra-sonication of ZIF-67 crystals results in ZIF-67 nano-flakes. ChemistrySelect 2016, 1, 59055908,  DOI: 10.1002/slct.201601513
    7. 7
      Beyer, S.; Schürmann, R.; Feldmann, I.; Blocki, A.; Bald, I.; Schneider, R. J.; Emmerling, F. Maintaining stable Zeolitic Imidazolate Framework (ZIF) templates during polyelectrolyte multilayer coating. Colloids Interface Sci. Ser. 2018, 22, 1417,  DOI: 10.1016/j.colcom.2017.11.004
    8. 8
      Buzanich, A. G.; Kulow, A.; Kabelitz, A.; Grunewald, C.; Seidel, R.; Chapartegui-Arias, A.; Radtke, M.; Reinholz, U.; Emmerling, F.; Beyer, S. Observation of early ZIF-8 crystallization stages with X-ray absorption spectroscopy. Soft Matter 2021, 17, 331334,  DOI: 10.1039/D0SM01356K
    9. 9
      Phan, A.; Doonan, C. J.; Uribe-Romo, F. J.; Knobler, C. B.; O’Keeffe, M.; Yaghi, O. M. Synthesis, Structure, and Carbon Dioxide Capture Properties of Zeolitic Imidazolate Frameworks. Acc. Chem. Res. 2010, 43, 5867,  DOI: 10.1021/ar900116g
    10. 10
      Tan, J. C.; Bennett, T. D.; Cheetham, A. K. Chemical structure, network topology, and porosity effects on the mechanical properties of Zeolitic Imidazolate Frameworks. Proc. Natl. Acad. Sci. U.S.A. 2010, 107, 99389943,  DOI: 10.1073/pnas.1003205107
    11. 11
      Chen, B.; Yang, Z.; Zhu, Y.; Xia, Y. Zeolitic imidazolate framework materials: recent progress in synthesis and applications. J. Mater. Chem. A 2014, 2, 1681116831,  DOI: 10.1039/C4TA02984D
    12. 12
      Spencer, E. C.; Angel, R. J.; Ross, N. L.; Hanson, B. E.; Howard, J. A. K. Pressure-Induced Cooperative Bond Rearrangement in a Zinc Imidazolate Framework: A High-Pressure Single-Crystal X-Ray Diffraction Study. J. Am. Chem. Soc. 2009, 131, 40224026,  DOI: 10.1021/ja808531m
    13. 13
      Hu, Y.; Kazemian, H.; Rohani, S.; Huang, Y.; Song, Y. In situ high pressure study of ZIF-8 by FTIR spectroscopy. Chem. Commun. 2011, 47, 1269412696,  DOI: 10.1039/c1cc15525c
    14. 14
      Hu, Y.; Liu, Z.; Xu, J.; Huang, Y.; Song, Y. Evidence of pressure enhanced CO2 storage in ZIF-8 probed by FTIR spectroscopy. J. Am. Chem. Soc. 2013, 135, 92879290,  DOI: 10.1021/ja403635b
    15. 15
      Widmer, R. N.; Lampronti, G. I.; Anzellini, S.; Gaillac, R.; Farsang, S.; Zhou, C.; Belenguer, A. M.; Wilson, C. W.; Palmer, H.; Kleppe, A. K. Pressure promoted low-temperature melting of metal–organic frameworks. Nat. Mater. 2019, 18, 370376,  DOI: 10.1038/s41563-019-0317-4
    16. 16
      Choi, J.; Im, J.; Noh, K.; Kim, J.; Vogt, T.; Lee, Y. Universal Gas-Uptake Behavior of a Zeolitic Imidazolate Framework ZIF-8 at High Pressure. J. Phys. Chem. C 2019, 123, 2576925774,  DOI: 10.1021/acs.jpcc.9b08539
    17. 17
      Formalik, F.; Mazur, B.; Fischer, M.; Firlej, L.; Kuchta, B. Phonons and Adsorption-Induced Deformations in ZIFs: Is It Really a Gate Opening?. J. Phys. Chem. C 2021, 125, 79998005,  DOI: 10.1021/acs.jpcc.1c01342
    18. 18
      Iacomi, P.; Maurin, G. ResponZIF Structures: Zeolitic Imidazolate Frameworks as Stimuli-Responsive Materials. ACS Appl. Mater. Interfaces 2021, 13, 5060250642,  DOI: 10.1021/acsami.1c12403
    19. 19
      Eslava, S.; Zhang, L.; Esconjauregui, S.; Yang, J.; Vanstreels, K.; Baklanov, M. R.; Saiz, E. Metal-organic framework ZIF-8 films as low-κ dielectrics in microelectronics. Chem. Mater. 2013, 25, 2733,  DOI: 10.1021/cm302610z
    20. 20
      Pimentel, B. R.; Parulkar, A.; Zhou, E.-k.; Brunelli, N. A.; Lively, R. P. Zeolitic imidazolate frameworks: next-generation materials for energy-efficient gas separations. ChemSusChem 2014, 7, 32023240,  DOI: 10.1002/cssc.201402647
    21. 21
      Hoop, M.; Walde, C. F.; Riccò, R.; Mushtaq, F.; Terzopoulou, A.; Chen, X.-Z.; deMello, A. J.; Doonan, C. J.; Falcaro, P.; Nelson, B. J. Biocompatibility characteristics of the metal organic framework ZIF-8 for therapeutical applications. Appl. Mater. Today 2018, 11, 1321,  DOI: 10.1016/j.apmt.2017.12.014
    22. 22
      Dai, H.; Yuan, X.; Jiang, L.; Wang, H.; Zhang, J.; Zhang, J.; Xiong, T. Recent advances on ZIF-8 composites for adsorption and photocatalytic wastewater pollutant removal: Fabrication, applications and perspective. Coord. Chem. Rev. 2021, 441, 213985,  DOI: 10.1016/j.ccr.2021.213985
    23. 23
      Hu, W.; Pattengale, B.; Huang, J. Zeolitic imidazolate frameworks as intrinsic light harvesting and charge separation materials for photocatalysis. J. Chem. Phys. 2021, 154, 240901,  DOI: 10.1063/5.0048720
    24. 24
      Knebel, A.; Caro, J. Metal–organic frameworks and covalent organic frameworks as disruptive membrane materials for energy-efficient gas separation. Nat. Nanotechnol. 2022, 17, 911923,  DOI: 10.1038/s41565-022-01168-3
    25. 25
      Paul, A.; Banga, I. K.; Muthukumar, S.; Prasad, S. Engineering the ZIF-8 Pore for Electrochemical Sensor Applications - A Mini Review. ACS Omega 2022, 7, 2699327003,  DOI: 10.1021/acsomega.2c00737
    26. 26
      Tsang, C. Y.; Cheung, M. C. Y.; Beyer, S. Assessing the colloidal stability of copper doped ZIF-8 in water and serum. Colloids Surf., A 2023, 656, 130452,  DOI: 10.1016/j.colsurfa.2022.130452
    27. 27
      Wang, H.; Zhao, L.; Xu, W.; Wang, S.; Ding, Q.; Lu, X.; Guo, W. The properties of the bonding between CO and ZIF-8 structures: a density functional theory study. Theor. Chim. Acta 2015, 134, 3139,  DOI: 10.1007/s00214-015-1636-4
    28. 28
      Yilmaz, G.; Peh, S. B.; Zhao, D.; Ho, G. W. Atomic-and Molecular-Level Design of Functional Metal–Organic Frameworks (MOFs) and Derivatives for Energy and Environmental Applications. Adv. Sci. 2019, 6, 1901129,  DOI: 10.1002/advs.201901129
    29. 29
      Liu, Y.; Huo, Y.; Wang, X.; Yu, S.; Ai, Y.; Chen, Z.; Zhang, P.; Chen, L.; Song, G.; Alharbi, N. S. Impact of metal ions and organic ligands on uranium removal properties by zeolitic imidazolate framework materials. J. Clean. Prod. 2021, 278, 123216,  DOI: 10.1016/j.jclepro.2020.123216
    30. 30
      Wang, H.; Pei, X.; Kalmutzki, M. J.; Yang, J.; Yaghi, O. M. Large Cages of Zeolitic Imidazolate Frameworks. Acc. Chem. Res. 2022, 55, 707721,  DOI: 10.1021/acs.accounts.1c00740
    31. 31
      Bhattacharyya, S.; Han, R.; Kim, W.-G.; Chiang, Y.; Jayachandrababu, K. C.; Hungerford, J. T.; Dutzer, M. R.; Ma, C.; Walton, K. S.; Sholl, D. S. Acid Gas Stability of Zeolitic Imidazolate Frameworks: Generalized Kinetic and Thermodynamic Characteristics. Chem. Mater. 2018, 30, 40894101,  DOI: 10.1021/acs.chemmater.8b01394
    32. 32
      Sarkar, S.; Grønbech, T. B. E.; Mamakhel, A.; Bondesgaard, M.; Sugimoto, K.; Nishibori, E.; Iversen, B. B. X-ray Electron Density Study of the Chemical Bonding Origin of Glass Formation in Metal–Organic Frameworks. Angew. Chem., Int. Ed. 2022, 61, e202202742  DOI: 10.1002/anie.202202742
    33. 33
      Butler, K. T.; Worrall, S. D.; Molloy, C. D.; Hendon, C. H.; Attfield, M. P.; Dryfe, R. A.; Walsh, A. Electronic structure design for nanoporous, electrically conductive zeolitic imidazolate frameworks. J. Mater. Chem. C 2017, 5, 77267731,  DOI: 10.1039/C7TC03150E
    34. 34
      Möslein, A. F.; Tan, J.-C. Vibrational Modes and Terahertz Phenomena of the Large-Cage Zeolitic Imidazolate Framework-71. J. Phys. Chem. Lett. 2022, 13, 28382844,  DOI: 10.1021/acs.jpclett.2c00081
    35. 35
      Möslein, A. F.; Donà, L.; Civalleri, B.; Tan, J.-C. Defect Engineering in Metal–Organic Framework Nanocrystals: Implications for Mechanical Properties and Performance. ACS Appl. Mater. Interfaces 2022, 5, 63986409,  DOI: 10.1021/acsanm.2c00493
    36. 36
      Graiver, D.; Farminer, K.; Narayan, R. A review of the fate and effects of silicones in the environment. J. Polym. Environ. 2003, 11, 129136,  DOI: 10.1023/A:1026056129717
    37. 37
      Du, X.; Zhou, M. Strategies to enhance catalytic performance of metal-organic frameworks in sulfate radical-based advanced oxidation processes for organic pollutants removal. Chem. Eng. J. 2021, 403, 126346,  DOI: 10.1016/j.cej.2020.126346
    38. 38
      Lewis, D. W.; Ruiz-Salvador, A. R.; Gómez, A.; Rodriguez-Albelo, L. M.; Coudert, F.-X.; Slater, B.; Cheetham, A. K.; Mellot-Draznieks, C. Zeolitic imidazole frameworks: structural and energetics trends compared with their zeolite analogues. CrystEngComm 2009, 11, 22722276,  DOI: 10.1039/b912997a
    39. 39
      Li, K.; Olson, D. H.; Seidel, J.; Emge, T. J.; Gong, H.; Zeng, H.; Li, J. Zeolitic Imidazolate Frameworks for Kinetic Separation of Propane and Propene. J. Am. Chem. Soc. 2009, 131, 1036810369,  DOI: 10.1021/ja9039983
    40. 40
      Amrouche, H.; Aguado, S.; Pérez-Pellitero, J.; Chizallet, C.; Siperstein, F.; Farrusseng, D.; Bats, N.; Nieto-Draghi, C. Experimental and Computational Study of Functionality Impact on Sodalite–Zeolitic Imidazolate Frameworks for CO2 Separation. J. Phys. Chem. C 2011, 115, 1642516432,  DOI: 10.1021/jp202804g
    41. 41
      Chaplais, G.; Fraux, G.; Paillaud, J.-L.; Marichal, C.; Nouali, H.; Fuchs, A. H.; Coudert, F.-X.; Patarin, J. Impacts of the Imidazolate Linker Substitution (CH3, Cl, or Br) on the Structural and Adsorptive Properties of ZIF-8. J. Phys. Chem. C 2018, 122, 2694526955,  DOI: 10.1021/acs.jpcc.8b08706
    42. 42
      Yagi, R.; Ueda, T. Substitution (CH 3, Cl, or Br) effects of the imidazolate linker on benzene adsorption kinetics for the zeolitic imidazolate framework (ZIF)-8. Phys. Chem. Chem. Phys. 2023, 25, 2058520596,  DOI: 10.1039/D3CP01662E
    43. 43
      Amrouche, H.; Creton, B.; Siperstein, F.; Nieto-Draghi, C. Prediction of thermodynamic properties of adsorbed gases in zeolitic imidazolate frameworks. RSC Adv. 2012, 2, 60286035,  DOI: 10.1039/c2ra00025c
    44. 44
      Dürholt, J. P.; Fraux, G.; Coudert, F.-X.; Schmid, R. Ab Initio Derived Force Fields for Zeolitic Imidazolate Frameworks: MOF-FF for ZIFs. J. Chem. Theory Comput. 2019, 15, 24202432,  DOI: 10.1021/acs.jctc.8b01041
    45. 45
      Morris, W.; Stevens, C. J.; Taylor, R. E.; Dybowski, C.; Yaghi, O. M.; Garcia-Garibay, M. A. NMR and X-ray Study Revealing the Rigidity of Zeolitic Imidazolate Frameworks. J. Phys. Chem. C 2012, 116, 1330713312,  DOI: 10.1021/jp303907p
    46. 46
      Ke, Q.; Duan, Y.; Ji, Y.; Zhao, D.; Zhang, H.; Duan, C.; Li, L.; Wei, Y. Identical Composition and Distinct Performance: How ZIF-8 Polymorphs’ Structures Affect the Adsorption/Separation of Ethane and Ethene. J. Chem. Eng. Data 2021, 66, 34833492,  DOI: 10.1021/acs.jced.1c00322
    47. 47
      Hinuma, Y.; Pizzi, G.; Kumagai, Y.; Oba, F.; Tanaka, I. Band structure diagram paths based on crystallography. Comput. Mater. Sci. 2017, 128, 140184,  DOI: 10.1016/j.commatsci.2016.10.015
    48. 48
      Hohenberg, P.; Kohn, W. Inhomogeneous Electron Gas. Phys. Rev. 1964, 136, B864B871,  DOI: 10.1103/PhysRev.136.B864
    49. 49
      Kohn, W.; Sham, L. J. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A1133A1138,  DOI: 10.1103/PhysRev.140.A1133
    50. 50
      Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 38653868,  DOI: 10.1103/PhysRevLett.77.3865
    51. 51
      Giannozzi, P.; Baseggio, O.; Bonfà, P.; Brunato, D.; Car, R.; Carnimeo, I.; Cavazzoni, C.; de Gironcoli, S.; Delugas, P.; Ferrari Ruffino, F. Quantum ESPRESSO toward the exascale. J. Chem. Phys. 2020, 152, 154105,  DOI: 10.1063/5.0005082
    52. 52
      Blöchl, P. E. Projector augmented-wave method. Phys. Rev. B 1994, 50, 1795317979,  DOI: 10.1103/PhysRevB.50.17953
    53. 53
      Dal Corso, A. Pseudopotentials periodic table: From H to Pu. Comput. Mater. Sci. 2014, 95, 337350,  DOI: 10.1016/j.commatsci.2014.07.043
    54. 54
      Bader, R.; Zou, P. An atomic population as the expectation value of a quantum observable. Chem. Phys. Lett. 1992, 191, 5458,  DOI: 10.1016/0009-2614(92)85367-J
    55. 55
      Otero-de-la-Roza, A.; Johnson, E. R.; Luaña, V. Critic2: A program for real-space analysis of quantum chemical interactions in solids. Comput. Phys. Commun. 2014, 185, 10071018,  DOI: 10.1016/j.cpc.2013.10.026
    56. 56
      Yu, M.; Trinkle, D. R. Accurate and efficient algorithm for Bader charge integration. J. Chem. Phys. 2011, 134, 064111,  DOI: 10.1063/1.3553716
    57. 57
      Guilherme Buzanich, A.; Radtke, M.; Yusenko, K. V.; Stawski, T. M.; Kulow, A.; Cakir, C. T.; Röder, B.; Naese, C.; Britzke, R.; Sintschuk, M. BAMline─A real-life sample materials research beamline. J. Chem. Phys. 2023, 158, 244202,  DOI: 10.1063/5.0157194
    58. 58
      Ravel, B.; Newville, M. ATHENA, ARTEMIS, HEPHAESTUS: data analysis for X-ray absorption spectroscopy using IFEFFIT. J. Synchrotron Radiat. 2005, 12, 537541,  DOI: 10.1107/S0909049505012719
    59. 59
      Cui, K.; Schmidt, J. Enabling efficient and accurate computational studies of MOF reactivity via QM/MM and QM/QM methods. J. Phys. Chem. C 2020, 124, 1055010560,  DOI: 10.1021/acs.jpcc.0c01220
    60. 60
      Hartmann, M.; Clark, T.; van Eldik, R. Hydration and water exchange of zinc (II) ions. Application of density functional theory. J. Am. Chem. Soc. 1997, 119, 78437850,  DOI: 10.1021/ja970483f
    61. 61
      Bock, C. W.; Katz, A. K.; Markham, G. D.; Glusker, J. P. Manganese as a replacement for magnesium and zinc: functional comparison of the divalent ions. J. Am. Chem. Soc. 1999, 121, 73607372,  DOI: 10.1021/ja9906960
    62. 62
      Sandmark, C.; Brändén, C. I.; Kjekshus, A.; Åkeson, Å.; Theorell, H.; Blinc, R.; Paušak, S.; Ehrenberg, L.; Dumanović, J. The Crystal Structure of Hexaimidazole Zinc (II) Dichloride Tetrahydrate, Zn(C3H4N2)6Cl2·4H2O. Acta Chem. Scand. 1967, 21, 993999,  DOI: 10.3891/acta.chem.scand.21-0993
    63. 63
      Torzilli, M. A.; Colquhoun, S.; Doucet, D.; Beer, R. H. The interconversion of dichlorobis (Nn-propylsalicylaldimine) zinc (II) and bis (Nn-propylsalicylaldiminato) zinc (II). Polyhedron 2002, 21, 697704,  DOI: 10.1016/S0277-5387(02)00839-2
    64. 64
      Boada, R.; Chaboy, J.; Hayama, S.; Keenan, L. L.; Freeman, A. A.; Amboage, M.; Díaz-Moreno, S. Unraveling the Molecular Details of the “Gate Opening” Phenomenon in Zif-8 with X-Ray Absorption Spectroscopy. J. Phys. Chem. C 2022, 126, 59355943,  DOI: 10.1021/acs.jpcc.2c00373
    65. 65
      Cocchi, C.; Prezzi, D.; Ruini, A.; Caldas, M. J.; Molinari, E. Optical properties and charge-transfer excitations in edge-functionalized all-graphene nanojunctions. J. Phys. Chem. Lett. 2011, 2, 13151319,  DOI: 10.1021/jz200472a
  • Supporting Information

    Supporting Information


    The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.jpcc.3c06054.

    • Lattice parameters of all the investigated structures; average bond lengths of all the considered structures before and after optimization; partial charges and band gaps, including the energies of the frontier levels; band structure plots (PDF)


    Terms & Conditions

    Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.