A novel series of isoreticular metal organic frameworks: realizing metastable structures by liquid phase epitaxy.

A novel class of metal organic frameworks (MOFs) has been synthesized from Cu-acetate and dicarboxylic acids using liquid phase epitaxy. The SURMOF-2 isoreticular series exhibits P4 symmetry, for the longest linker a channel-size of 3 × 3 nm(2) is obtained, one of the largest values reported for any MOF so far. High quality, ab-initio electronic structure calculations confirm the stability of a regular packing of (Cu(++))(2)- carboxylate paddle-wheel planes with P4 symmetry and reveal, that the SURMOF-2 structures are in fact metastable, with a fairly large activation barrier for the transition to the bulk MOF-2 structures exhibiting a lower, twofold (P2 or C2) symmetry. The theoretical calculations also allow identifying the mechanism for the low-temperature epitaxial growth process and to explain, why a synthesis of this highly interesting, new class of high-symmetry, metastable MOFs is not possible using the conventional solvothermal process.

Highly crystalline molecular scaffolds with nm-sized pores offer an attractive basis for the fabrication of novel materials. The scaffolds alone already exhibit attractive properties, e.g., the safe storage of small molecules like H2, CO2 and CH4 1. In addition to the simple release of stored particles, changes in the optical and electronic properties of these nanomaterials upon loading their porous systems with guest molecules offer interesting applications, e.g., for sensors2 or electrochemistry3. For the construction of new nanomaterials, the voids within the framework of nanostructures may be loaded with nm-sized objects, such as inorganic clusters, larger molecules, and even small proteins4, a process that holds great potential, as, for example, in drug release5 or the design of novel battery materials6. Furthermore, meta-crystals may be prepared by loading metal nanoparticles into the pores of a three-dimensional scaffold to provide materials with a number of attractive applications ranging from plasmonics789 to fundamental investigations of the electronic structure and transport properties of the meta-crystals10.

In the last two decades, numerous studies have shown that MOFs, also termed porous coordination polymers (PCPs), fulfill many of the aforementioned criteria. Although originally developed for the storage of hydrogen, e.g., in tanks of automobiles, MOFs have recently found more technologically advanced applications as sensors2, matrices for drug release5, as membranes for enantiomer separation11 and for proton conductance in fuel cells12. Although pore-sizes available within MOFs13 are already sufficiently large to host metal-nanoclusters such as Au55 1415, the actual realization of meta-crystals requires, in addition to structural requirements, a strategy for the controlled loading of the nanoparticles (NPs) into the molecular framework. Simply adding NPs to the solutions before starting the solvothermal synthesis of MOFs has been reported16, as well as the controlled encapsulation of nanoparticles17, but these procedures do not allow for the fabrication of nanoparticle superlattices with translational order. The layer-by-layer growing process of SURMOFs allows the uptake of nanosized objects during synthesis, including the fabrication of compositional gradients of different NPs18. Pure metal, metal oxide or covalent cluster NPs have already been proposed for the fabrication of novel electronic devices10, but one may speculate that also the embedding of larger one-dimensional objects such as nanowires, carbon or inorganic nanotubes will be of interest.

The introduction of procedures for epitaxial growth of MOFs on functionalized substrates19 was a major step towards controlled functionalization of frameworks. This liquid phase epitaxy (LPE) process allowed the preparation of compositional gradients2021 and provided a new strategy to load nano-objects into predefined pores.

Unfortunately, the LPE process has so far been only demonstrated for a fairly small number of MOFs11192223, none of which exhibits the required pore sizes. In addition, for the fabrication of meta-crystals, the architecture of the network should be sufficiently adjustable to realize pores of different sizes. There should also be a straightforward way to functionalize the framework itself in order to tailor the interaction of the walls with the targeted guest objects. Ideally, such a strategy would be based on an isoreticular series of MOFs in which the pore size can be determined by a choice from a homologous series of ligands with different lengths13.

Here, we present a novel isoreticular MOF series with the required flexibility as regards pore-sizes and functionalization, which is well suited for the LPE process. This class, denoted as SURMOF-2, is derived from MOF-2, one of the simplest framework architectures2425. MOF-2 is based on paddle-wheel (pw) units formed by attaching 4 dicarboxylate groups to Cu++ or Zn++-dimers, yielding planar sheets with 4-fold symmetry (see Fig. 1). These planes are held together by strong carboxylate/metal-bonds26. Conventional solvothermal synthesis yields stacks of pw-planes shifted relative to each other, with a corresponding reduction in symmetry to yield a P2 or C2 structure242728.

Figure 1

Schematic representation of the synthesis and formation of the SURMOF-2 analogues.

Results

Experiments

The relative shifts between the pw-planes can be avoided when using the recently developed liquid phase epitaxy (LPE) method22. As demonstrated by the XRD-data shown in Fig. 2, for a number of different dicarboxylic acid ligands, the LPE-process yields surface attached metal organic frameworks (SURMOFs) with P4 symmetry29, except for the 2,6-NDC ligand, which exhibits a P2 symmetry as a result of the non-linearity of the carboxyl functional groups on the ligand. Solvothermal synthesis conditions using either DMF or ethanol as solvent yield bulk MOF powder, which does not exhibit the high-symmetry P4 structure. Instead the structure obtained for the QPDC ligand corresponds to the P2 structure also seen for MOF-2 (see Fig. S3, supporting information). The ligand PPDC, pentaphenyldicarboxylic acid, with a length of 2.5 nm, to our knowledge, is among the longest ligands which have so far been used successfully for a MOF synthesis, only in the very recent work by Deng et al.30 this length has been exceeded. A detailed analysis of the diffraction data allowed us to propose the structures shown in Fig. 1. This novel isoreticular series of MOFs, hereafter termed SURMOF-2, also consists of pw-planes, which – in contrast to MOF-2 – are stacked directly on-top of each other. This absence of a lateral shift when stacking the pw-planes yields 1d-pores of quadratic cross section with a diagonal up to 4 nm (see Fig. 2). Importantly, for the SURMOF-2 series, interpenetration is absent. For many known isoreticular MOF series, the formation of larger and larger pores is limited by this phenomenon: if the pores become too large, a second or even a third 3d-lattice will be formed inside the first lattice yielding interpenetrated networks, which exhibit the expected larger lattice constant but with rather small pore sizes31, the hosting of NPs thus becomes impossible. For SURMOF-2, like for other MOF-structures with the appropriate topology (non self-complementrary nets, see Ref. 30 and references cited there) the presence of 1d-pores of quadratic cross section is not compatible with a second, interwoven network and, as a result, interpenetration is suppressed.

Figure 2

Out of plane XRD data of Cu-BDC, Cu-2,6-NDC, Cu-BPDC, Cu-TPDC, Cu-QPDC, Cu-P(EP)2DC and Cu-PPDC (upper left), schematic representation (upper right) and proposed structures of SURMOF-2 analogues (lower panel).

All the SURMOF-2 are grown on –COOH terminated SAM surface using the LPE method.

Although the solubility of the long dicarboxylic acids TPDC, QPDC and PPDC, (see Fig. 1) used as ligands for the SURMOF-2 fabrication is fairly small, we encountered no problems in the LPE process. This can be understood by considering that already small concentrations of dicarboxylic acids are sufficient for the formation of a single monolayer on the substrate, the elementary step of the LPE synthesis process32. Even with the longest dicarboxylic acid available to us, PPDC, growth occurred in a straightforward fashion and optimization of the growth process was not necessary.

Theory

The P4 structures proposed for the SURMOF-2 series, except for the 2,6-NDC ligand, differ markedly from the bulk MOF-2 structures obtained from conventional solvothermal synthesis242728. To understand this unexpected difference, and in particular the absence of relative shifts between the pw-planes in the proposed SURMOF-2 structure, we performed extensive quantum chemical calculations, employing an approximate variant of density-functional theory (DFT), London dispersion-corrected self-consistent charge density-functional based tight-binding (DFTB)333435. We studied a periodic model of MOF-2 as well as different SURMOF derivatives.

Periodic structure DFTB calculations (Fig. S1) confirmed that the P2/m structure, proposed by Yaghi et al. for the bulk MOF-2 material24, is indeed the energetically most favorable configuration for MOF-2, while the C2/m structure, that has also been reported for bulk MOF-22728 is slightly higher in energy (41 meV per formula unit, see Table 1). The optimum interlayer distance between the pw-planes in the P4 structure is 5.8 Å, in very good agreement with the experimental value of 5.6 Å (as obtained from the in-plane XRD-data, see Fig. S2). It is important to note that at that distance the linkers cannot remain in rectangular position as in MOF-2, but are twisted in order to avoid steric hindrance and to optimize linker-linker interactions.

Table 1

Stacking energy and geometries of P4 SURMOF-2 derivatives

	Symmetry	a = c	b	Stacking Energy*
Cu2(BDC)2	C2	11.19	5.0	−0.76
Cu2(BDC)2	P2	11.19	5.4	−0.8
Cu2(BDC)2	P4	11.19	5.8	−0.59
Cu2(BDC)2	P2	13.35	5.6	−0.4
Cu2(BPDC)2	P4	15.49	5.9	−0.68
Cu2(TPDC)2	P4	19.84	5.9	−0.91
Cu2(QPDC)2	P4	24.24	5.9	−1.21
Cu2(P(EP)2DC)2	P4	25.12	5.2	−1.73
Cu2(PPDC)2	P4	28.59	5.9	−1.45

*Stacking energies (DFTB level, in eV) of SURMOFs. The energies were calculated within periodic boundary conditions and are given per formula unit.

The P4/mmm structure, as proposed for SURMOF-2, is less stable by 200 meV per formula unit as compared to the P2/m configuration. Although the crystal energy for the P4 stacking is substantially smaller than that for the P2 packing, simulated annealing (DFTB model using Born-Oppenheimer Molecular Dynamics) carried out at 300 K for 10 ps demonstrated that the P4/mmm structure corresponds to a local minimum. This observation is confirmed by the calculation of the stacking energy of MOF-2, where a shifting of the layers is simulated along the P2-P4-C2 path (Figure 3). On the basis of these calculations, we thus propose that SURMOF-2 adopts the metastable P4 structure.

Figure 3

Energy for the relative shift of each layer in two different directions: horizontal (P4 to P2) and diagonal (P4 to C2) at a fixed interlayer distance of 5.6 Å.

Structures have been partially optimized and energies are given per formula unit with respect to fully dissociated planes.

In Table 1, the interlayer interaction energies are summarized. They range from 590 meV per formula unit with respect to fully dissociated MOF-2 layers for Cu2(bdc)2 and successively increase for longer linkers of the same type, reaching 1.45 eV per formula unit for Cu2(ppdc)2, indicating that the linkers play an important role in the overall stabilization of SURMOF-2 with large pore sizes. Even stronger interlayer interactions are found for different linker topologies (PPDC). Additional calculations carried out reveal that the presence of ligands like DMF (dimethylformamide) or water destabilizes the P4 structure by inducing shifts between the pw-planes. A detailed computational analysis of the role of the individual building blocks of SURMOF-2 for the crystal growth and stabilization as well as the role of solvent molecules will be published elsewhere.

Discussion

The flat, well-defined organic surfaces used as the templating (or nucleating) substrate in the LPE growth process provide an excellent explanation for why SURMOF-2 grows with the highly symmetrical P4 structure instead of adopting the energetically more favorable, low symmetry bulk P2 (or C2) symmetry3236. The carboxylic acid groups exposed on the surface force the first layer of deposited Cu-dimers into a coplanar arrangement, which is not compatible with the bulk P2 (or C2) structure, but, rather, nucleates the P4 packing. In turn, the carboxylate-groups in the first layer of deposited dicarboxylic acids provide the same constraints for the next layer (see Fig. 1). As a result of the layer-by-layer method employed for further SURMOF-2 growth, the same boundary conditions apply for all subsequent layers. The organic surface thus acts as a seed structure to yield the metastable P4 packing, not an unusual motif in epitaxial growth37.

Furthermore, the calculations allow us to predict that it will be possible to grow SURMOF structures with even larger linker molecules as used here. The π stacking will further stabilize the P4 symmetry as the interaction energy per formula unit is more and more dominated by the linkers (Table 1). At present, we cannot predict the maximum possible size of the pores, one may speculate that by adding additional functionalities to the linkers the linker-linker interaction may by sufficiently large to exceed 10 nm × 10 nm and thus reach a size regime which is of interest for the realization of plasmonic structures by embedding metal nanoparticles6.

Methods

Computational details

All structures were created using a preliminary version of our topological framework creator software, which allows the creation of topological network models in terms of secondary building units and their replacement by individual molecules to create the coordinates of virtually any framework material. The generated starting coordinates, including their corresponding lattice parameters, have been hierarchically fully optimized using the Universal Force Field (UFF)38, followed by the dispersion-corrected self-consistent-charge density-functional based tight-binding (DFTB) method3334353940 that we have recently validated against experimental structures of HKUST-1, MOF-5, MOF-177, DUT-6 and MIL-53(Al), as well as for their performance to compute the adsorption of water and carbon monoxide34. For all calculations, we employed the deMonNano software41. We have chosen periodic boundary conditions for all calculations, and the repeated slab method has been employed to compute the properties of the single layers in order to evaluate the stacking energy. A conjugate-gradient scheme was employed for geometry optimization of atomic coordinates and cell dimensions. The atomic force tolerance was set to 3 × 10−4 eV/Å.

The confined conditions at the SURMOF-substrate interface was modeled by fixing the corresponding coordinate in the computer simulations. All calculated structures have been substantiated by simulated annealing calculations, where a Born-Oppenheimer Molecular Dynamics trajectory at 300 K has been computed for 10 ps without geometry constrains. All structures remained in P4 topology.

Experimental methods

The SURMOF-2 samples were grown on Au substrates (100-nm Au/5-nm Ti deposited on Si wafers) using a high-throughput approach, spray method42. The gold substrates were functionalized by self-assembled monolayers, SAMs, of 16-mercaptohexadecanoic acid (MHDA). These substrates were mounted on a sample holder and subsequently sprayed with a 1 mM solution of Cu2(CH3COO)4·H2O in ethanol and an ethanolic solution of SURMOF-2 organic linkers (0.1 mM of BDC, 2,6-NDC, BPDC and saturated solutions of TPDC, QPDC, PPDC, P(EP)2DC) at room temperature. The number of spray cycles employed was different for different linkers (PPDC: 90, TPDC: 80, QPDC: 50 and P(EP)2DC: 50, BDC, 2,6-NDC and BPDC: all 20). After synthesis all samples were characterized with X-ray diffraction (XRD). The solvothermal synthesis condition for (Cu-QPDC) MOF-2 were as those reported for (Cu, BDC) MOF-224. Ethanol and DMF were used as sovents.

Author Contributions

J.L., O.S., H.K.A. prepared the SURMOFs investigated in this study. S.B., S.G., X.F. and K.M. synthesized the organothiols used to fabricate the SURMOFs. B.L., I.-B. M. and T.H. carried out the calculations. Data analysis was carried out by P.W., H.G. and C.W. The work was directed by A.G., K.M., T.H. and C.W. All authors contributed equally in writing the manuscript.