Uncovering growth species of multivariate MOFs in liquid phase by mass spectrometry

Jinli Han, Suming Chen, Xiaohun Zhou, Hexiang Denga,,∗

a Key Laboratory of Biomedical Polymers-Ministry of Education, College of Chemistry and Molecular Sciences, Wuhan University, Wuhan 430072, China

b Division of Advanced Nanomaterials, Suzhou Institute of Nano-tech and Nano-bionics, Chinese Academy of Sciences, Suzhou 215125, China

c The Institute for Advanced Studies, Wuhan University, Wuhan 430072, China

ABSTRACT The unveiling of MOF growth mechanism is hampered by the lack of fundamental knowledge about the very early stage of nucleation, especially the form and ratio of molecular species in the solution for crystal growth.Herein, we report the detection of growth species for a series of MOFs with mono-linker, Cu-MOF-2-BDC and Cu-MOF-2-NDC, and two linkers, MTV-MOF-2-(C4H4)x, by high resolution ESI-MS, where a large variety of Cu-containing species are identified unambiguously.The solvent molecules such as H2O,methanol and DMF participate in the formation of these species, other than ethanol.Furthermore, in the growth solution of MTV-MOF-2-(C4H4)x, growth species containing two different organic linkers are observed.The feeding ratio is not the only factor controlling the distribution of growth species for MTVMOFs, but also the solvent involves in coordination, an aspect usually overlooked previously.

Keywords:Metal-organic framework Multivariate MOFs Mass spectrometry Growth species Self-assembly

Metal-organic framework (MOF) is a fast developing class of porous crystals constructed by linking organic linkers and metal containing joints by strong coordination bonds.These frameworks offer molecularly defined pore spaces within which matter can be manipulated and controlled with the implications on their gas storage and separation, chemical sensing and energy applications[1–6].Multivariate MOF (MTV-MOF), is constructed by mixing organic linkers of identical size and coordination but carrying different functional groups.These functional groups exist in the same MOF single crystal without altering the framework topology, therefore the properties of MTV-MOFs are much more tuneable than the mono-linker counterparts [7], giving opportunities to create new pore environments for potential applications in catalytic and biological properties [8–10].However, the spatial arrangements of multiple types of functional groups or metal clusters in bulk crystals are complicated, which is reflected in the discrepancy between feeding ratio and experimental ratio for the synthesis [11].The spatial arrangements of organic linkers and metals are gradually unveiled by solid-state nuclear magnetic resonance (NMR) spectroscopy [12], X-ray photoelectron spectroscopy and atom probe tomography [13,14], however, it remains unclear that how these different organic linkers assemble in the solution and form molecular species at early nucleation stage in MTV-MOF growth solution.

Recently, bothex situandin situtechniques have been developed to investigate MOF crystallization mechanisms, with the attempt to offer details to understand MOF formation, morphology evolution, and structure transition [15].Electron microscopy[16–18], X-ray diffraction [19,20], and atomic force microscopy were also applied for the study of crystal nucleation and growth(Scheme 1) [21,22].However, these techniques did not distinguish the molecular species in the growth solution, especially for bulk crystal of MTV-MOFs (Scheme 1).Qualitative analysis of these chemicals in growth solution has the potential to offer valuable information about basic structural units, which is favorable for understanding the initial coordination process of metal ions and organic linkers.Thanks to the developments of electrospray ionization mass spectrometry (ESI-MS), the information of molecular species in solution was obtained clearly for MOFs with monolinker [23–27].Here, this technique is demonstrated to effectively distinguish the different type of growth species with their characteristic molecular weights in the growth solution of MTV-MOFs(Scheme 1).

Scheme 1.Experimental approaches to study MTV-MOFs.Shown is the advantage of ESI-MS technique to uncover the types of growth species existing in MTV-MOF growth solution.

In this study, we report the use of high resolution ESI-MS for identification of metal-containing molecular species of different types for MOFs with mono-linker, Cu-MOF-2-BDC, Cu-MOF-2-NDC and Cu-MOF-2-BDC-Cl2, and two linkers, MTV-MOF-2-(C4H4)x(xis varied between 0 and 1), where multiple different types of metalcontaining molecular species were assigned accurately in their growth solutions, including CuL (L: organic linker), CuL2, Cu2L,Cu2L2and so on.Growth species of the same type but with different molecular weight were also observed, due to the coordination of various solvent molecules.Moreover, the molecular species were detected in MTV-MOF growth solution, and two different organic linkers could coordinate with the same metal ion to give different Cux(LinkerA)y(LinkerB)z-type molecular species.In addition,the discrepancy between feeding ratio and experimental ratio of different organic linkers is found to be related to the contents of different growth species of the same type, which is influenced by the solvents and coordination ability of organic linkers during the synthesis.

Specifically, we use a series of two dimensional (2D) porous MOFs with isoretucular structure (the same topology), Cu-MOF-2, to illustrate the ESI-MS method.The structure and property of MOFs were characterized prior to the ESI-MS analysis.Cu-MOF-2-BDC [28], constructed from paddle-while secondary building units (SBUs) and benzene-1,4-dicarboxylate (BDC) ligands to form a simple framework (Fig.S1 in Supporting information), is known to have isoreticular structures of various organic linker types,e.g., naphthalene-1,4-dicarboxylic acid (NDC) for Cu-MOF-2-NDC and 2,5-dichloroterephthalic acid (BDC-Cl2) for Cu-MOF-2-BDC-Cl2[29,30].Cu-MOF-2-BDC crystals were synthesized using the similar procedure as reported with modifications for improving crystal crystallization rate.Mixed solvents (methanol, ethanol andN,N-dimethylformamide (DMF) with volumetric ratio of 1:1:4)were used rather than only pure DMF.The morphology of the Cu-MOF-2-BDC crystals was revealed by scanning electron microscopy(SEM), which showed a regular square shaped plate (Fig.S2 in Supporting information).The powder X-ray diffraction (PXRD) patterns of the as-synthesized crystals had a slightly shift comparing with the pattern simulated from single crystal structure data, due to the effects of the coordination molecules, where the solvents in our system were different from only pure DMF as solvent in previous study (Fig.S5 in Supporting information) [28].Similar results were also observed in Cu-MOF-2-NDC system, where the regular square plated-shaped crystals were exhibited in SEM image (Fig.S3 in Supporting information), and PXRD patterns were in good agreement with simulated pattern demonstrating that this MOF possessed the well-defined 2D network structure, same like Cu-MOF-2-BDC (Fig.S6 in Supporting information).It was worth noting that the PXRD patterns of activated Cu-MOF-2-BDC and Cu-MOF-2-NDC also matched well, due to their identical topology (Figs.S1 and S7 in Supporting information).But their PXRD patterns differed from the corresponding as-synthesized samples, consistent with previous studies [28].The structural transformation can be attributed the removal of solvent molecule from MOF crystal structure, which is common in MOF studies.The calculated Brunauer-Emmett-Teller(BET) surface areas of Cu-MOF-2-BDC and Cu-MOF-2-NDC were 150 and 12 m2/g, respectively.The higher permanent porosity of the guest-free MOFs was observed in Cu-MOF-2-BDC due to the smaller steric hindrance of BDC in comparison to NDC (Figs.S1 and S8 in Supporting information).

In order to study the relationship between crystal structure and molecular species in MOF crystal formation process,it is critical to provide standard samples for comparison.First of all, the growth solution of pure Cu-MOF-2-BDC was detected by high-resolution ESI-MS, where a large variety of Cucontaining species were identified unambiguously.The mother liquor for MOF synthesis was filtered with nylon filters before injecting into the high-resolution ESI-MS.Multiple new peaks were identified in either the positive- or negative-ion spectra of mother liquor, other than in reference solutions (solvent, metal ions and organic linkers solutions), demonstrating the self-assembling of multiple Cu-containing molecular species in MOF growth solution (Figs.S11 and S13 in Supporting information).Accurate assignments of these growth species reveal their detailed chemical formulas with the reference to the theoretical isotopic distributions (Fig.1, Fig.S12, Tables S2 and S3 in Supporting information).The relatively more abundant CuBDCcontaining molecular species in Cu-MOF-2-BDC mother liquor were [Cu(C8H5O4)(HCON(CH3)2)2]+, [Cu(C8H5O4)(HCON(CH3)2)]+,and [Cu(C8H5O4)(H2O)(HCON(CH3)2)]+in the positive-ion spectra, as well as [Cu(C8H5O4)(NO3)2]−and [Cu(C8H5O4)2(NO3)]−in the negative-ion spectra, respectively.The most abundant of these molecular species were [Cu(C8H5O4)(HCON(CH3)2)2]+and[Cu(C8H5O4)(NO3)2]−, in the positive-ion and negative-ion spectra,respectively.Surprisingly, the organic linker BDC containing two -COOH functional groups was only deprotonated on the side where the Cu2+was coordinated, while the exposed -COOH on the other side was not deprotonated, giving a CuL type species, for instance[Cu(C8H5O4)(NO3)2]−in negative ion spectra (Fig.1A, I).If the -COOH on both sides of BDC were deprotonated, these two -COOH would have been coordinated both with Cu2+at once to form Cu2L type species, such as [Cu2(C8H4O4)(NO3)3]−(Fig.1A, III).Similarly,the Cu coordination sites in CuL type species could also capture one organic linker to form CuL2type species (Fig.1A, II).Furthermore, Cu2L2type species,e.g., [Cu2(C8H4O4)(C8H5O4)(NO3)2]−could be formed from CuL2species through removing one proton H and combing with Cu2+(Fig.1A, VI).However, the low relative abundance of SBUs (Cu2L4, X in Fig.1A) implies they are unlikely to be the direct growth units for the growth of MOFs.

Moreover, small solvent molecules such as DMF, methanol,H2O and anion NO3−also coordinated with metal-containing species to give multiple growth species of same type, but with different molecular weight, [Cu2(C8H4O4)(NO3)3]−and[Cu2(C8H4O4)(CH3OH)(NO3)3]−(Fig.1A, III and IV) for instance.However, another common solvent molecule, ethanol, didn’t participate in the formation of these species, indicating its weaker coordination to the metal ions.This provides valuable guidance to the choice of solvents for the synthesis of MOF crystals and their activation.

Fig.1.High resolution ESI-MS spectra and assignments of different molecular species in Cu-MOF-2-BDC and Cu-MOF-2-NDC growth solutions.(A) Negative-ion spectrum of Cu-MOF-2-BDC obtained with full range of species and schematic diagrams of a series of Cu-containing molecular species (Cu2+:BDC=1:1, magnified image for blue dashed box from m/z 460–1130 for optimum display of speciation).The inset provides confirmation of the chemical composition for the [Cu(C8H5O4)(NO3)2]−species centered at m/z ∼352 through comparison of the measured and theoretical isotopic distributions.(B) Negative-ion spectrum of Cu-MOF-2-NDC obtained with full range of species and schematic diagrams of a series of Cu-containing molecular species (Cu2+:NDC=1:1, magnified image for blue dashed box from m/z 460–1130 for optimum display of speciation).The inset provides confirmation of the chemical composition for the [Cu(C12H7O4)(NO3)2]−species centered at m/z ∼402 through comparison of the measured and theoretical isotopic distributions.

In addition to Cu-MOF-2-BDC, the metal-containing molecular species of other isoreticular MOF, Cu-MOF-2-NDC, were also evaluated, where the types of molecular species were similar to those for Cu-MOF-2-BDC.It is worth mentioning that the exact masses of BDC and NDC are respectively 166.0266 and 216.0423, and the exact mass of NDC is 50.0157 greater than that of BDC.Thus, it is facile to assign these same types of molecular species of Cu-MOF-2-NDC from ESI-MS spectra based on the results of Cu-MOF-2-BDC.Specifically, the exact masses of CuL-type and CuL2-type species in Cu-MOF-2-NDC growth solution ought to be 50.0157 and 100.0314 greater than those of Cu-MOF-2-BDC, respectively(Tables S2 and S3).Indeed, all types of molecular species in Cu-MOF-2-BDC growth solution were also detected in Cu-MOF-2-NDC mother liquor unambiguously (Table S3).Similarly, the most abundant of these species were [Cu(C12H7O4)(HCON(CH3)2)2]+and[Cu(C12H7O4)(NO3)2]−, in the positive-ion and negative-ion spectra, respectively, and the measured isotopic distributions of each species matched well with the theoretical spectra, confirming the accuracy of molecular species assignments (Fig.1B).

Furthermore, this finding that the molecular species in the growth solution of these two isoreticular MOFs were similar to each other provides the chemical bases for the construction of MTV-MOFs with functionality variance [7].MTV-MOF-2-(C4H4)xwas chosen here for the study of growth species and the exploration of their the self-assembly from metal ions and linkers.Prior to study it, the characterizations of this MTV-MOF were carried out.In MTV-MOF-2-(C4H4)x, although the different functional groups (-H in BDC, -C4H4in NDC) appended to the same backbones, the framework topology did not change based on their PXRD patterns (Fig.S7).The crystals were synthesized using the similar experimental procedures of Cu-MOF-2-BDC and the feeding ratio of BDC to NDC was 1 (Section 1 in Supporting information).The MTV-MOF-2-(C4H4)xwas a regular square shaped plate,similar to that of Cu-MOF-2-NDC (Figs.S3 and S4 in Supporting information).The calculated BET surface area of MTV-MOF-2-(C4H4)xwas 94 m2/g, in between those of Cu-MOF-2-BDC and Cu-MOF-2-NDC (Fig.S8 in Supporting information).The exact ratio of organic linkers, BDC to NDC, in MTV-MOF-2-(C4H4)xcrystals was further obtained from the1H NMR spectroscopy of an HCl-digested solution of bulk samples, which was in a good agreement with elemental analysis results (Figs.S9, S10 and Table S1 in Supporting information).Nevertheless, a critical remain unanswered is why these different organic linkers can be self-assembly into the same MOF backbone without changing the framework topology.

In the attempt to answer this challenging question, the mother liquor of MTV-MOF-2-(C4H4)xwas also detected by high resolution ESI-MS, where the molecular species in the growth solution were assigned precisely to give accurate chemical formulae(Fig.2).In addition to the same Cu-containing species observed in the study of pure Cu-MOF-2, new molecular species containing two different organic linkers were observed.Typically, most of molecular species of Cu-MOF-2-BDC and Cu-MOF-2-NDC were observed in growth solution of MTV-MOF-2-(C4H4)xexcept some CuBDC-containing growth species, indicating that the coordination mode of the metal ion and the organic linkers remained unchanged under the same synthetic condition (Figs.2A and C).But the strength of coordination formed using different organic linkers was different.The more abundant CuL-containing molecular species in growth solutions were [Cu(C12H7O4)(HCON(CH3)2)2]+and [Cu(C8H5O4)(HCON(CH3)2)2]+in the positive-ion spectra, as well as [Cu(C12H7O4)(NO3)2]−and [Cu(C8H5O4)(NO3)2]−in the negative-ion spectra, while the most abundant species were[Cu(C12H7O4)(HCON(CH3)2)2]+and [Cu(C12H7O4)(NO3)2]−, in the positive-ion and negative-ion spectra of MTV-MOF-2-(C4H4)x, respectively.The results showed that the relative content of CuNDCcontaining molecular species in growth solution was higher than that of CuBDC-containing molecular species when the feeding ratio of BDC to NDC was 1, without considering the difference in ionization efficiency of similar growth species.This indicated that the strength of coordination formed by NDC is greater than that of BDC.This is in good accordance to the experimental ratio of BDC and NDC in bulk crystal obtained from the1H NMR spectroscopy,where the experimental ratio of BDC to NDC was less than 1 with equal amount of linkers in starting materials (Fig.3A and Fig.S10).Moreover, two different organic linkers, BDC and NDC can coordinate with the same Cu2+to give the CuxLBDCLNDC-type species,which were detected in positive and negative-ion spectra of MTVMOF-2-(C4H4)x(II, VI, VII, VIII, IX, X and XI in Fig.2C).Therefore,CuL2-type species possessed three forms with different molecular weight, including Cu(LBDC)2, Cu(LNDC)2and Cu(LBDC)(LNDC), and the relative content of these molecular species varied with the feeding ratio of BDC to NDC (Figs.2B and 3B).Based on the results of CuL2-type species, similarly, the Cu2L3-type species should include four different forms, if BDC and NDC coordinated with the same Cu2+.Indeed, Cu2(LBDC)2(LNDC), Cu2(LBDC)(LNDC)2and Cu2(LNDC)3were detected except Cu2(LBDC)3in growth solution of MTV-MOF-2-(C4H4)x.The absence of Cu2(LBDC)3is likely attributed to the weak coordination for BDC, further confirming the stronger coordination of NDC with Cu2+in comparison to that of BDC (Fig.2C,VIII).

Fig.2.High resolution ESI-MS spectra and assignments of molecular species in MTV-MOF-2-(C4H4)x growth solutions.(A) Negative-ion spectrum of MTV-MOF-2-(C4H4)x obtained with full range of species (Cu2+:organic linkers=1:1 and BDC:NDC=1:1, magnified image for blue dashed box from m/z 440–1130 for optimum display of speciation).The inset provides confirmation of the chemical composition for the [Cu(C8H5O4)(C12H7O4)(NO3)]−species centered at m/z ∼505 through comparison of the measured and theoretical isotopic distributions.(B) Negative-ion spectra of CuL2-type and Cu2L2-type species.(C) The detailed assignments of all growth species in negative-ion spectra of MTV-MOF-2-(C4H4)x.Green and orange backgrounds in (C) correspond to the species in (B).

In order to further study the impact from different feeding ratio of organic linkers for these different CuL2-type species, high resolution ESI-MS spectra were also obtained for MTV-MOF-2-(C4H4)xwith different feeding ratios of BDC to NDC (Figs.S15 and S17 in Supporting information).There was no obvious change in the types of growth species, however, the relative content of these growth species varied drastically (Figs.3B-D and Figs.S15-S17 in Supporting information).In order to compare the changing trends of growth species under different conditions, CuL2-type species was chosen here to illustrate the influence of different feeding ratios, where the relative abundance of the most abundant of molecular species in each spectrum was defined as 100% (Fig.3B).Specifically, when the feeding ratio of BDC to NDC was 3:1,the fraction of Cu(LBDC)2species present the most within CuL2-type species, while that of Cu(LNDC)2species was the least, indicating that the probability of Cu2+coordinating with BDC was greater than that with NDC, when BDC is more abundant in the solution (Fig.3B).These results were in stark contrast to those of the feeding ratio of BDC to NDC was 1:3, where the fractions of Cu(LBDC)2and Cu(LNDC)2species were respectively the least and most within CuL2-type species.It was not surprising that NDCcontaining molecular species were major species in the case of excess NDC.However, when the feeding ratio of BDC to NDC was 1, the fraction of Cu(LNDC)2species was also much greater than that of Cu(LBDC)2species, demonstrated that Cu2+tended to coordinate with NDC rather than BDC in growth solution.Moreover,other molecular species such as Cu2L2-type species also fellow above rules, further confirming the coordination of NDC with Cu2+is stronger than that of BDC.This finding was also revealed from the positive-ion spectra, where the fraction of Cu2(LBDC)2species decreased with the ratio of NDC increasing (Fig.S16).This demonstrated that BDC in Cu(LBDC)2and Cu(LBDC)(LNDC) species prefers to be replaced by NDC to respectively form Cu(LBDC)(LNDC) and Cu(LNDC)2species in the case of abundant NDC, confirming the difference in coordination ability of similar organic linkers.More interestingly, the small solvent also plays a critical role during the assembling process, while the addition of methanol favors to the formation of CuNDC containing species at low NDC to BDC ratio,but this role is not obvious at high NDC to BDC ratio (Fig.3A).This shows that solvent molecule participate in the MOF formation process, and will change the content of linkers in the resulting MTVMOFs.In short, the solvent is considered as a critical factor to control the ratio of linkers in MTV-MOFs, which was usually ignored in the early conducted research studies.

Fig.3.The change trends of NDC actual ratio in its bulk crystals and relative contents of different CuL2-type and Cu2L2-type species in MTV-MOF-2-(C4H4)x growth solution.(A) The change trends of NDC actual ratio with different feeding ratio in MTV-MOF-2-(C4H4)x bulk crystals.(B-D) Relative contents of CuL2-type and Cu2L2-type species in growth solution at different feeding ratios of BDC to NDC.

Furthermore, the initial growth solution of MTV-MOF-2-(C4H4)xwas also assessed by ESI-MS at room temperature, where the types of molecular species in growth solution were similar to those of sample synthesized 1 h at 70 °C, but the relative abundance of molecular species in these two experiments was different, confirming that the temperature and reaction time has an effect on selfassembly of species.The solution was prepared by mixing metal ions and organic linkers directly at room temperature, the similar negative-ion spectrum demonstrated that the formation these molecular species would be formed once the metal ions contacted with the organic linkers, which is a fast process (Fig.S17).

In addition, the ESI-MS was further successfully applicable to probe another isoreticular MOF,e.g., Cu-MOF-2-BDC-Cl2, where the types of molecular species in growth solution were similar to those of other two isoreticular MOFs.The MOF growth solution was preparedviaa similar condition of Cu-MOF-2-BDC (Section 1 in Supporting information).Specifically, the exact mass of BDC-Cl2is 233.9487 and 67.9221 greater than that of BDC, thus these same types of growth species in spectra of isoreticular Cu-MOF-2-BDCCl2could be assign facilely.Indeed, multiple types of molecular species in Cu-MOF-2-BDC growth solution were also detected in Cu-MOF-2-BDC-Cl2growth solution with organic linkers BDC were replaced by BDC-Cl2and the most abundant of these species were[Cu(C8H3Cl2O4)(HCON(CH3)2)2]+and [Cu(C8H3Cl2O4)(NO3)2]−, in the positive-ion and negative-ion spectra, respectively (Figs.S18 and S19 in Supporting information).The similar results of molecular species detected by high resolution ESI-MS demonstrated the universality of this useful technique for tracking molecular species in MOF growth solution.

In summary, a large variety of metal-containing growth species of a series of MOFs with mono-linker and two linkers were identified unambiguously by high resolution ESI-MS.Small solvent molecules such as H2O, methanol and DMF participated in the coordination of growing species.The relative content of different organic linkers in bulk MTV-MOF crystals, originated from content of organic linkers in various molecular species in the growth solution, was influenced by the feeding ratio and the strength of coordination bond form by different linkers.One important observation is that the solvent molecules participate in the MOF formation process, and will change the content of linkers in the resulting MTV-MOFs.These findings, based on exploration in coordination molecules of metal ions and organic linkers by high resolution ESI-MS, provide new insight to understand the formation of isoreticular MOFs and MTV-MOFs.

Declaration of competing interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Acknowledgments

We thank W.Wang, S.Li and J.Liu for their valuable advice on this research.We also thank the Core Research Facilities of College of Chemistry and Molecular Sciences at Wuhan University for assistance on material characterizations, The State Key Laboratory of Coordination Chemistry (SKLCC) at Nanjing University for high resolution ESI-MS characterizations.We acknowledge support from the National Natural Science Foundation of China (Nos.91545205,91622103, 21971199 and 22025106), National Key Research and Development Project (Nos.2020YFB1505700, 2016YFA0200700 and 2018YFA-0704000), Innovation Team of Wuhan University (No.2042017kf0232) and the Fundamental Research Funds for the Central Universities (No.2042020kf0214).

Supplementary materials

Supplementary material associated with this article can be found, in the online version, at doi:10.1016/j.cclet.2021.11.035.