The Lewis acidic and basic character of the internal HKUST-1 surface determined by inverse gas chromatography †

438 | CrystEngComm, 2015, 17, 438–447 This journal is © The R Institute of Physical Chemistry, Technical University Bergakademie Freiberg, Leipziger Straße 29, 09599 Freiberg, Germany. E-mail: florian.mertens@chemie.tu-freiberg.de † Electronic supplementary information (ESI) available: Additional experimental details, retention times, linear regression functions, determined enthalpies, entropies, and other quantities are listed in the ESI. The characterization of HKUST-1 used for the adsorption experiments by PXRD, BET is also described in the ESI. See DOI: 10.1039/c4ce01327a Cite this: CrystEngComm, 2015, 17, 438


Introduction
Until now, the synthesis of a very large number of metalorganic frameworks with different and partially exotic topologies, that display, if the inapplicability of the standard BET analysis for microporous systems is ignored, BET values of up to 10 000 m 2 g −1 , is still the main objective of many research groups. [1][2][3][4][5][6][7] Beside synthesis, many potential applications of this class of materials connected to the structural variability in combination with crystallinity and permanent porosity as well as the possibility to tune pore size and geometry for specific applications are conceivable and were intensively discussed in recent years. [8][9][10][11][12][13][14][15][16][17] A specific field of MOF application that allows to utilize these features are the various types of chromatography, such as liquid chromatography 16,17 and gas chromatography. [18][19][20] Recently, the separation of different analyte mixtures, such as those of xylene isomers, 19,21 n-alkanes, 20,22,23 polychlorinated biphenyls, 24 polycyclic aromatic hydrocarbons, 24 and branched alkanes 23,25 by MOF based capillary GC was demonstrated. In contrast to the interest in the separations, there have been almost no studies using the so called inverse gas chromatography with MOFs, a term that is used for GC investigations which do not primarily focus on separation problems but on the evaluation of physico-chemical properties of the material which the stationary phase is made of. 20,21,[26][27][28][29][30] HKUST-1 is one of the best investigated metalorganic frameworks in respect to crystal structure, properties, and possible applications. 21,[31][32][33][34] This MOF, originally reported by Chui et al., 35 consists of copperĲII)-dimer units, where each copper ion is coordinated to four oxygen atoms of the organic linker benzene-1,3,5-tricarboxylate, forming a three dimensional crystalline and porous framework with three cavities possessing the diameters 4.7 Å, 10.3 Å, and 12.1 Å. 36 One coordination site of the metal ions is unsaturated and therefore free to accommodate Lewis basic molecules, like water. 35 The coordinated water molecules can easily be removed at mild conditions (by heating in vacuum or inert gas) from the coordination shell of the Cu 2+ ions without any loss of porosity and crystallinity creating an unsaturated coordination site, a procedure that is called activation. 37 Beside these Lewis acidic centers of the copper(II)-SBUs, HKUST-1 possesses also Lewis basic areas like the π-electron rich benzene rings of the linker molecules benzene-1,3,5-tricarboxylate.
Molecules, like aromatic compounds, can be coordinated by weak interactions such as π-π stacking and C-H-π interactions as largely explored in the field of supramolecular chemistry. 38,39 The interaction between analytes and the crystalline lattice can be classified in several ways depending on the properties of the analyte. Analytes possessing a permanent dipole moment are characterized in addition to the nonspecific part of the van der Waals interaction (interactions with a 1/r 6 dependence), which is the London-or dispersion interaction, 40 by the so called specific part which should be in the case of HKUST-1 mainly the Debye-interaction. 41 This type of interaction is defined as the angle averaged interaction between a permanent and an induced dipole. The other contribution to the specific part of the van der Waals interaction is the so called Keesom interaction, 42 which results from averaged interaction between two permanent dipoles.
Since the linker molecules and metal centers in many MOFs can be exchanged and new functionalities can be introduced to some degree in a predictable manner (crystal engineering), it appears to be important that the effect of such modifications can be quantitatively evaluated. In this contribution we therefore report a procedure how specific and nonspecific interaction between the Lewis acidic (open metal sites of the copper-SBUs) and the Lewis basic (π-electron rich aromatic rings) parts of HKUST-1 and different analytes can be analyzed by means of inverse gas chromatography.

Theoretical backgroundinverse gas chromatography
Inverse gas chromatography 43,44 has been widely applied in the determination of physico-chemical parameters in very different systems. From the measurement of the temperature and carrier gas flow dependent retention times, quantities characterizing the stationary phaseanalyte interaction such as the thermodynamic adsorption enthalpy, ΔH ads , and adsorption entropy, ΔS ads , as well as kinetic parameters, such as mass transport and diffusion coefficients can be obtained. Once a routinely applicable procedure to prepare chromatographic columns and capillaries is established, IGC is an efficient method to investigate adsorption phenomena in porous media and thus in metal-organic frameworks as well.

Determination of thermodynamic parameter values of the adsorption process
By plotting the logarithm of the retention volume gained from isothermal gas chromatographic measurements versus the inverse temperature, the enthalpy and entropy of adsorption are easily obtained from if the first term can be neglected. A condition that is usually fulfilled in gas chromatographic experiments (n s being the total amount of the adsorbed analyte). 45,46 A mathematical derivation of the relationship (1) in case of gas-solid chromatography can be found in the paper by Katsanos et al. 45 and in the ESI. † The retention volume V N was obtained from the adjusted retention times t r ′ = (t r − t m ), the gas flow F a , and the application of the two correction terms reflecting the adjustment of the gas to the column temperature and the pressure drop experienced during the passage through the column: 47,48 T a stands for ambient temperature and j is the James-Martin gas compressibility correction factor.
Starting from the desire to obtain heats of adsorption that can be compared to calorimetric measurements, which are usually performed under isochoric conditions, the heat of adsorption derived from eqn (1), which represents an enthalpy change, needs to be converted to a heat of adsorption that represents an internal energy change by |ΔH iso | = |ΔU diff | + RT av (for details see ESI †).

IGC for the determination of surface interaction Gibbs free energies
Following the terminology used by Fowkes 49 the interaction of an adsorbate can be split in the so called nonspecific interaction Gibbs free energy ΔG nonspec caused by the Londondispersion interaction and a term denoted as specific interaction Gibbs free energy ΔG spec caused by the Keesom-and Debye-interactions as well as all Gibbs free energy contributions caused by acid-basic interactions, hydrogen bonding, π-π-interactions, and steric effects: Measuring the retention times as a function of the temperature, the interaction Gibbs free energy can be determined by calculating the net retention volume according to where C is a constant that depends on the reference state. 50,51 V N is calculated by using eqn (2).
Calculation of the nonspecific ΔG nonspec and specific component ΔG spec of the interaction Gibbs free energy of adsorption ΔG interact In case that a polar analyte interacts with the surface of solids, dispersive and specific interactions occur and the problem arises how to split the nonspecific part of the interaction Gibbs free energy from the specific part. Both types of interactions, the specific and the nonspecific one, cannot be investigated separately. In the literature different approaches to solve this problem can be found 50,[52][53][54][55][56] that are all based on the same evaluation model. At first, the (linear) relationship between the Gibbs free energy of interaction, ΔG interact , of n-alkanes possessing a specific property ( f.e. the polarizability) that can be expected to scale with the nonspecific interaction strength is determined as a reference function. n-Alkanes are chosen because they display only nonspecific interactions. The choice of the specific property is the quantity that distinguishes the different reported approaches. In the next step, the reference function value at the value of the specific property of the polar analyte to be investigated is taken and the difference to ΔG interact of the analyte to be investigated is calculated. Under the assumption that the dependence of the nonspecific interaction does strictly correlate with the specific property of the analytes, no matter what sort of molecule it is, the calculated difference is fully ascribed to the specific interactions. Fig. 1 illustrates the described procedure.
Using the polarizability α as the specific property bears the advantage that this quantity practically does not vary with temperature, 52 which is why it was chosen in this work for the analysis. The relationship between polarizability, the London dispersion interaction, and the corresponding −ΔG L can be described by the expression (see Donnet et al. 52 ): From eqn (3) and (4) then follows: , is characteristic for a given solid surface and is related to the London dispersive interaction component of this surface. 52 The required electronic deformation polarizabilities can be taken from suitable references, 57 calculated either by using the Debye equation or by the increment method according to Miller et al. 58 Characterization of the surface donor-acceptor properties by IGC Plotting the interaction Gibbs free energy ΔG spec as a function of temperature, the specific enthalpy ΔH spec and the specific entropy ΔS spec of adsorption can be obtained from the intercept and the slope of the Gibbs Helmholtz equation: The specific enthalpy of adsorption of an analyte on a specific surface of a solid can be correlated to the acid-base properties of this solid, here the stationary phase of the chromatographic capillaries used. 50 The quantitative description of the donor and acceptor properties of the surface requires the determination of the acid and base constant K A and K B . The relation to the specific interaction is given according to Gutmann 59 by with DN and AN denoting the respective Gutmann donor and acceptor numbers. The values of these quantities are known for many solvents (analytes). 59 If the specific enthalpy of adsorption is known, the values for K A and K B , can be calculated. In order to increase the precision, it is recommendable to use a multitude of analytes and to determine K A and from the slope and the intercept of the linear regression plot:

HKUST-1 preparation
For the infrared spectroscopy and powder X-ray crystallography investigations of the analyte -HKUST-1 interaction, the use of a completely desolvated HKUST-1 material possessing high surface area is required. In a typical synthesis a solution of 0.210 g (1 mmol) benzene-1,3,5-tricarboxylic acid (ABCR, 98%) dissolved in 25 mL ethanol (VWR, >99.8%) was dropwise added at room temperature to a solution of 0.299 g (1.5 mmol) copper(II)-acetate monohydrate (Acros Organics, >98%) in 25 mL deionized water. Afterwards the suspension was stirred for 1 hour and then the blue precipitate was filtered off and washed with ethanol and diethyl ether. For activation, the HKUST-1 material was heated to 373 K under vacuum conditions. From the procedure described, a dark blue powder was obtained. Powder X-ray diffraction, IR-spectroscopy, and the modified BET surface area measurements for microporous systems were conducted. The results can be found in the ESI. †

Capillary coating
All experiments were conducted with a gas chromatographic silica capillary of 30 m length and an inner diameter of 0.25 mm. In order to apply the coating procedure used, the silica surface of the column has to be terminated with hydroxyl groups, which was achieved by first pumping a 1.0 M solution of sodium hydroxide through the capillary for a period of 60 min followed by an intensive rinsing with water until the pH-value of the solution at the outlet indicated neutrality. The same procedure was repeated with 0.1 M solution of hydrochloric acid. Finally the column was dried in a stream of argon (excess pressure at the inlet 0.2 bar) at a temperature of 423 K for 24 hour. 21 After this surface activation procedure was conducted, HKUST-1 was deposited by a cyclic layer deposition procedure based on the controlled SBU-approach. For this procedure, 21 a solution, called the precursor solution, of 299 mg (1.5 mol) copperĲII)-acetate monohydrate in 150 mL ethanol and a solution, called the linker solution, of 210 mg (1 mmol) benzene-1,3,5-tricarboxylic acid dissolved in the same amount of ethanol were prepared. These solutions were pumped with an excess pressure at the inlet of 5 bar alternately through the capillary using a self-built coating apparatus (ESI †). Between every cycling step the liquid film was conditioned for 2 min followed by a rinsing or purging procedure with ethanol executed by the coating apparatus in order to avoid the formation of small HKUST-1 particles inside the steel pipes. The whole coating procedure progressed by the following sequence, which was 15 times repeated: 1 mL precursor solution (5 bar) → 2 min conditioning → cleaning of the coating apparatus → 1 mL precursor solution (5 bar) → 2 min conditioning → cleaning of the coating apparatus.
After this procedure the film was homogenized and strongly fixed onto the capillary walls by the following conditioning sequence: 2 hours argon (0.2 bar excess pressure) at 333 K → 2 hours argon (0.2 bar excess pressure) at 353 K → 2 hours argon (0.2 bar excess pressure) at 373 K.

Chromatographic setup
All GC experiments were performed using an HP 5890 series II gas chromatograph equipped with a flame ionization detector (FID). Helium (Praxair, 99.999%) was the carrier gas in all cases. The data evaluation was conducted with the EZChrome software from Agilent Technologies.
The chromatograms for the determination of the adsorption enthalpies and entropies of cyclohexane and benzene were recorded using pure substances only and taken at a constant inlet gas flow of 5.31 mL min −1 at different temperature that varied from 443.15 K to 473.15 K. The chromatograms for the determination of the adsorption enthalpies and entropies of the ethers were also recorded using pure substances but at a constant inlet gas flow of 4.15 mL min −1 and temperatures that varied from 493.15 K to 523.15 K. During the gas chromatographic measurements the ambient temperature and pressure was 298 K and 0.969 bar, respectively. The dead time was determined by the retention time of methane.

Inductively coupled plasma optical emission spectroscopy (ICP-OES)
For the determination of a film thickness profile on the basis of the deposited amount of HKUST-1, a capillary with an inner diameter d c of 0.25 mm and a length L of 30 m was coated according to the procedure described above. After purging of the coated layer with deionized water, a piece of approximately 10 mm was cut off from each side of the column for collecting scanning electron microscope (SEM) images and energy dispersive X-ray spectroscopic (EDX) data. Afterwards the capillary was cut into pieces of 2 m each. The deposited HKUST-1 was removed from every section by very pure concentrated hydrochloric acid (suprapore, Merck) and the concentration of the copperĲII)-ions in the obtained solutions was determined using a Varian Vista-Pro-CCD Simultaneous ICP-OES apparatus. The amount of HKUST-1 in each of the original capillary sections was then calculated from the obtained concentrations.

Characterization of the HKUST-1 coatings inside a capillary
In contrast to the published procedures to deposit a thin film of a metal-organic framework in capillaries made from fused silica by other authors, we chose the approach to coat GC capillaries with HKUST-1 based on the controlled SBU approach (CSA). Previous processes were based either on the solvothermal synthesis of little MOF-crystallites and the following deposition of the particles from a suspension onto the capillary wall via the interaction with the surface hydroxyl groups 19,22 or on a cyclic layer deposition technique in which the MOF structure is formed directly during the coating procedure inside the capillary via a precipitations reaction anchoring the MOF layer via a bifunctional self-assembled monolayer (SAM) onto the fused silica. 20,21,60 Because the pipes of the coating apparatus become easily congested with the precipitate applying this method, we switched to the simpler anchoring technique using the hydroxyl groups of the fused silica surface in combination with the CSA method as described in the section Materials and methods.
After the coating of the capillary, SEM-images from the inner capillary surface near the inlet (Fig. 2a) and near the outlet (Fig. 2b) were taken from which respective film thicknesses of 200-300 nm and 150-200 nm were determined, indicating that more HKUST-1 is deposited at the beginning of the capillary than at the end. In order to track the decrease of the film thickness, the amount of HKUST-1 was determined as described earlier by acid treatment and ICP-OES measurements. The MOF coating thickness profile obtained is displayed in Fig. 2c and shows the presence of an approximately constant film thickness in the capillary sections from 2 m to 26 m measured from the capillary entry. The increased amount at the end of the capillary in such systems is in the literature often attributed to the solution acceleration at the end of the capillary. 19 Because of the comparably small amount of HKUST-1 on the fused silica surface, the small dimension, and the curved surface of the capillary, a direct structural verification of HKUST-1 on the substrate wall by PXRD was not possible. We therefore, as an indirect verification, separately mixed test samples of the reaction solutions to be used for the coating experiment and investigated the precipitate by PXRD. In addition the reaction solutions leaving the capillary outlet were collected in one flask and the precipitate formed there was investigated by PXRD as well. In both cases the diffractograms of HKUST-1 were obtained (see ESI †).

Gas chromatographic separation
With the coated capillary (L = 30 m, d c = 0.25 mm), a mixture of benzene (bp = 80.1°C) and cyclohexane (bp = 80.7°C) was separated at 473 K with a helium flow of 4.37 mL min −1 . In order to avoid overcharging effects, only the methane enriched gas phase over a mixture of these liquid analytes was injected into the gas chromatograph. As can be seen in Fig. 3, the identification of the two detected peaks is easily carried out by the comparison of their retention times with those of the pure substances. The chromatogram shows that benzene displays a higher adjusted retention time (t r ′ = 0.88 min) than its completely hydrogenated analog (t r ′ = 0.67 min). If the separation mechanism is mainly dominated by differing London dispersive interactions, the order of the analyte retention times follows the order of their boiling points. From the fact that this order is reversed, one can assume that a type of specific interaction, such as π-π stacking interactions between the Lewis basic π-electron rings of the HKUST-1 linker, benzene-1,3,5tricarboxylate, and the π-electron density of the analyte benzene may be present what would not be the case with the nonaromatic cyclohexane.
As an example demonstrating the influence of the Lewis acidity of the open metal sites of HKUST-1 on chromatographic separations, we chose ethers with different alkyl chains as Lewis basic analyte portfolio. The ether oxygen atoms with their high electron densities should be capable of interacting with the free unsaturated coordination sites of the copper ions of the SBU of the investigated stationary phase. The corresponding chromatogram of a mixture of diisopropyl ether, diethyl ether, tetrahydrofuran, and di-n-propyl ether is also shown in Fig. 3. If the ether molecules were to interact mainly via dispersive interactions, they should more or less separate according to their boiling points, which would be in the order diethyl ether → tetrahydrofuran → diisopropyl ether → di-n-propyl ether in contrast to the measured sequence of diisoproyl ether → diethyl ether → tetrahydrofuran → di-n-propyl ether. Of particular interest is the difference between tetrahydrofuran and diisopropyl ether, which have similar boiling points, i.e. 65.8°C and 69°C, respectively. These significant differences point to the relevance of specific interactions. Considering the geometry of the analytes and the accessibility of the free electron pair of the oxygen atoms of the different analytes for the interaction with the open metal site of the MOF-SBU, it is conceivable that the electron density of the oxygen atom of the diisopropyl group is shielded more intensively by the branched alkyl groups than in the cases of diethyl ether and tetrahydrofuran. Similarly, one can assume that the electron density of the tetrahydrofuran oxygen is also less shielded than the one of diethyl ether because of the ring structure of the THF molecule.

Thermodynamic investigation of the adsorption process on HKUST-1
The thermodynamic parameters ΔG ads , ΔU diff , ΔH iso , and ΔS ads obtained from the van't Hoff plots in Fig. 4 are summarized in Table 1. Benzene displays a higher absolute value of the enthalpy of adsorption |ΔH ads | than cyclohexane indicating a stronger interaction presumably because of the π-π stacking interactions between the aromatic rings of the HKUST-1 linker and the benzene molecule. Given that π-π stacking interactions are comparable weak, they are app. 8-12 kJ mol −1 , the difference of 4 kJ mole appears to be in accordance with this mode of interaction. 61 The absolute value of ΔS ads for the adsorption of benzene on HKUST-1 is also slightly higher than for cyclohexane, which is somewhat counterintuitive given the more flexible nature of the cyclohexane molecule compared to the benzene one and the small difference in the adsorption enthalpies. If one adopts the assumption of the presence of π-π stacking interactions, the possibly very good alignment of the benzene ring of the MOF linker and the analyte may explain the observed ΔS ads values. In the case of the ether separations, it is noteworthy that, although diisopropyl ether eluates earlier than diethyl ether, its enthalpy of adsorption is higher. Since the retention time is determined by the adsorption Gibbs free energy ΔG ads , this effect must be caused by entropic contributions. It is comprehensible that the alkyl chains of diisopropyl ether, which are longer than the ones of diethyl ether and also branched, possess more degrees of freedom that are affected/restricted by the adsorption process. In the case of the two subsequently eluated compounds tetrahydrofuran and di-n-propyl ether,  both the enthalpy and entropy differences support the observed order.
Determination of the specific und nonspecific interaction Gibbs free energies of ether-adsorption on HKUST-1 Conducting analyses as described in the theoretical background section for the adsorption of the above listed ethers on the HKUST-1 surface, a division of the interaction Gibbs free energy in a part which describes the specific interaction, like donor-acceptor contacts, and a part caused by nonspecific, dispersive London-interactions can be obtained. For these analyses the method by Donnet et al., 52 i.e. choosing the scaled electronic polarizability Ĳhν) 0,5 α O,L as the specific property, was applied. A graphic summary for the four ethers investigated is given in Fig. 1 bottom.
For the estimate of the contribution of the nonspecific interaction of the ether analytes, the general relationship between the scaled electronic polarizability of molecules and nonspecific interactions with HKUST-1 was established using the unbranched alkanes n-pentane, n-hexane, and n-heptane. It is noteworthy that ΔG interact changes to negative values for longer chains. Since both the adsorption enthalpy ΔH ads and the adsorption entropy ΔS ads can be expected to be negative, the gain in enthalpy (absolute value) must outweigh the contribution to ΔG interact by the change of ΔS ads caused by every additional -CH 2unit at the given temperature.
The values for |−ΔG interact |, |−ΔG specific |, and |−ΔG nonspecific | of the four investigated ethers at 493.15 K are compiled in Table 2. The specific Gibbs free energies of adsorption of diethyl ether and tetrahydrofuran are higher than the one of diisopropyl ether.
According to Table 2, the eluation order of the first three compounds, which is given by |−ΔG interact |, is mainly determined by the specific interaction whereas the nonspecific part by itself would cause a reversed order. If one considers the geometric factors of the analysts, the result can be interpreted in the way that THF has a more exposed, i.e. less shielded by side chains, oxygen atom that can therefore more strongly interact with the open metal sites. In addition, this adsorption process is also entropically supported since the ring is already strongly restricted in its motion compared to the freely moving side chains of diethyl ether and diisopropyl ether. Hence, the adsorption process is accompanied with a rather low entropy change. In respect to the sequence of the first three eluated compounds, di-n-propyl breaks the trend of increasing specific interaction and decreasing nonspecific interaction. The specific enthalpy of adsorption for all ethers is app. between 50-100% of the nonspecific one (see Table 3). Due to the lack of data on the interaction of HKUST-1 with other oxygen containing organic molecules, one can only compare the interaction strength with small gas molecules such as water (51 kJ mol −1 ), CO 2 (20 kJ mol −1 ), and CO (24 kJ mol −1 ), 61 which is in the same range.
Determination of the donor-acceptor character of the HKUST-1 surface A quantitative measure for the donor and acceptor properties, i.e. the acid base characteristics, of the surface of the stationary phase, here HKUST-1, can be gained by determining the acid constant K A and base constant K B . In order to calculate these constants, the specific enthalpy of adsorption and the donor and acceptor numbers of Gutmann must be known. 50,59 Only for benzene, diethyl ether, and tetrahydrofuran we found corresponding values in the literature and thus the analysis had to be limited to these compounds. Using the Gibbs-Helmholtz equation directly or a van't Hoff analysis, the specific enthalpy and entropy of adsorption can be determined by a linear fitting procedure with respect to temperature of −ΔG spec . The values of −ΔG spec , −ΔH spec , and −ΔS spec for benzene calculated from the interaction Gibbs free energy of adsorption in the temperature range from 493.15 K to 523.15 K are given in the ESI. † It is important to note that the split of −ΔG spec into the two different parts, −ΔH spec and −ΔS spec , via temperature dependent measurements cures the problem of potential entropy contamination between the reference compounds and the investigated analytes. The term entropy contamination is supposed to indicate that although it can be assumed that the nonspecific interaction enthalpy scales generically with the specific property, here the polarizability, this cannot necessarily be assumed for the entropy. Therefore the particular adsorption entropies of the reference molecules influence the split of −ΔG interact into −ΔG spec and −ΔG nonspec . By subsequently splitting −ΔG spec into −ΔH spec and −ΔS spec the adsorption entropy of the reference molecules will only affect the determined value of −ΔS spec , and will leave −ΔH spec unaffected. As a consequence, the determination of K A and K B will be free of entropy contamination effects. With the values −ΔH spec and −ΔS spec of benzene, diethyl ether, and tetrahydrofuran and the corresponding acceptor and donor numbers (see ESI †), the constants K A and K B are available, as discussed in the theoretical section, from the corresponding linear fit, which is graphically depicted in Fig. 5 The study indicates that MIl-100 (Fe, Cr, Al) binds nitrogen based contaminants better than HKUST-1 while in the case sulphur ones the opposite appears to be true. It is therefore difficult to rank the MOFs in respect to their Lewis acidity. In addition, the results also proved to be very solvent dependent. In this respect, the presented IGC method displays itself, given the use of noble gases as analyte carriers, as a more general means for the determination of the Lewis acidity and basicity of MOFs and will allow to quantify subtle differences of MOF materials properties.

Conclusion
Beside the continued demonstration that HKUST-1 is a suitable material for gas chromatographic applications, the influence of donor and acceptor areas of the MOF-lattice on the adsorption of different analytes was studied using inverse gas chromatography. The determination of a London force governed reference system and the determination of the so called interaction Gibbs free energies of a set of analytes allowed the separation of their interactions into nonspecific and specific contributions. An analysis based on the use of the Gutmann donor and acceptor numbers of a specific analyte portfolio allowed a quantitative judgement of the Lewis acidity and basicity of the MOF material HKUST-1. The analysis applied to HKUST-1 can be seen as a procedure that can be used for the design of MOF materials with intentionally modified properties.

Acknowledgements
Financial support by the Deutsche Forschungsgemeinschaft is gratefully acknowledged (Priority Program 1362 "Porous Metal-Organic Frameworks"). A. M. thanks Stephan Glante for his experimental help with the chromatographic measurements of the ethers.