Loading of ionic compounds into metal – organic frameworks : a joint theoretical and experimental study for the case of La 3 + †

Crystalline, highly orientated surface-anchored MOF thin films, grown on Au substrates, were prepared using liquid-phase epitaxy (LPE). The successful loading of La ions into the Cu3(BTC)2 (HKUST-1) SURMOFs (surface-mounted metal–organic frameworks) was monitored using X-ray diffraction (XRD). Theoretical calculations using classical force-field based Monte Carlo simulations yield a structure with two La ions within the large Cu3(BTC)2 pores, in full agreement with experimental results on the composition of these films and the relative intensities of the XRD peaks. Implications of these findings for using MOF thin films for electronic applications are briefly discussed.


Introduction
Among highly porous crystalline materials, presently metalorganic frameworks (MOFs) [1][2][3] attract increasing amount of attention resulting from their flexible pore sizes and large loading capacity for guest species, 4 in particular metals and metal ions. 5uring the past decade, a huge number of metal@MOF systems (for metals like Au, 6 Zn, 7 Pd 8 ) have been reported, with the metal guests rendering fascinating properties [e.g., such as hydrogen storage, 9,10 catalysis (e.g.CO oxidation, 11,12 alcohol oxidation 13 )] to the host MOF material.Particularly attractive are the cases where the loading of the metal guest species not only changes the crystal structure of the host material but also modifies their properties (e.g.electrical conductivity, 14 capacity of ion absorption, 15 optical properties). 16,17For a number of applications, e.g.electric/ electronic devices, 18 information storage, 19 electrochemistry, 20 heterogeneous catalysts, 21 MOF thin films, 22,23 the ability to load this MOF-coatings with metals is important.
In the present work we focus on a special type of MOF coatings, surface-mounted metal-organic frameworks (SURMOFs).They represent a particularly important class of MOF thin films since they are crystalline, oriented and exhibit a well-defined thickness.
They can be manufactured with thicknesses in the micrometerregime, allowing the application of standard methods to determine their mechanical, 24 optical 17 and electrical and electrochemical 14,20 properties.The loading of Eu 3+ ions into HKUST-I SURMOFs has been shown to result in interesting optical properties, it could be shown that optical excitations resulting from absorption in the MOF ligands could be transferred to the embedded Eu ions. 25n the present paper we investigate the loading of La 3+ cations into HKUST-1 or Cu 3 (BTC) 2 (BTC, 1,3,5-benzenetricarboxylic acid), a popular MOF-material with fairly large (B1 nm) pore sizes. 26Since we target applications in electrochemistry, information storage and sensors, here we focus on MOF thin films, SURMOFs.Investigations are carried out using X-ray diffraction (XRD), while quantitative information on the Cu/La-ratio is obtained from ICP-OES (Inductively Coupled Plasma Optical Emission Spectrometry) and XPS (X-ray Photoelectron Spectrum) data.In addition, force-field calculations are carried out to obtain the precise positions of the metal cations and their corresponding counter ions in the lattice.Interestingly, these calculations reveal a rather pronounced binding energy of the metal ions within the MOF cavities brought about by the fairly substantial cations-BTC ligand interaction.Knowledge about the precise structure of metal ion-loaded MOF thin films is crucial for understanding electrical transport, 14 electrochemical behaviour 20 and metal ion diffusion 19 in MOFs and SURMOFs.by thermal vapor deposition of 100 nm gold (99.995%,Chempur) onto polished silicon wafers (Wacker) precovered with a 5 nm titanium adhesion layer.Evaporation was carried out at a pressure of 2 Â 10 À7 mbar and a deposition rate of 0.5 nm s À1 , yielding a root-mean-square roughness of the Au films of about 1 nm.The gold substrates were kept in an argon atmosphere until use.

Preparation of COOH-terminated self-assembled monolayers (SAMs)
The 16-mercaptohexadecanoic acid (MHDA, 99%, Aldrich) solution was prepared by dissolving MHDA in a 5% (by volume) solution of acetic acid in ethanol (Merck) to reach the desired concentration of 20 mM. 27A clean gold substrate was placed in this solution for 72 h and then rinsed with ethanol and gently dried under nitrogen flux.

Fabrication of Cu 3 (BTC) 2 SURMOFs
All Cu 3 (BTC) 2 thin films used in the present work were grown on modified Au substrates using the liquid-phase epitaxy (LPE) method, which yields highly oriented and well-defined MOF coatings and can be applied to a variety of substrates.In contrast to the conventional bulk solvothermal MOF synthesis, where reactants are mixed to initiate the reaction, the layer-bylayer growth of Cu 3 (BTC) 2 SURMOFs proceeds in two steps comprising a sequential spray solution of the individual reactants.For the present experiments the SURMOFs were grown on a COOH-terminated organic surface on Au.
The SURMOFs were fabricated using the following diluted ethanolic solutions: copper acetate hydrate (1 mM, Aldrich) and BTC (1,3,5-benzenetricarboxylic acid) (0.2 mM, Aldrich).For the SURMOF deposition, we used a spray system (Fig. S1, ESI †), as described in detail in an earlier publication. 28The spray times were 15 s for the copper acetate solution and 25 s for the BTC solution.Each spray step was followed by a rinsing step (3 s) with pure ethanol to remove residual reactants.A total of 20 growth cycles were used for all SURMOFs investigated in this work.Before further experiment, all SURMOF samples were activated by ultrasound in dichloromethane solution for 5 min to remove residual solvent from the SURMOF pores and finally checked by XRD.
Preparation of the La(OTf ) 3 @Cu 3 (BTC) 2 SURMOF First, a Cu 3 (BTC) 2 thin film on the gold substrate was put into a 250 ml flask and then evacuated to 0.2 kPa at room temperature for 30 min.Subsequently, the sample was immersed in a freshly prepared solution of La(OTf ) 3 [lanthanum(III) trifluoromethanesulfonate] in ethanol (1 mM, Aldrich) kept at 65 1C.After an immersion time of 12 h the sample was removed from the solution, rinsed with pure ethanol, and finally dried in a flux of nitrogen gas.
The characteristic R-SO 3 À vibration at 1042 cm À1 appears after immersing the Cu 3 (BTC) 2 SURMOFs in a La(OTf ) 3 solution (Fig. S2, ESI †), this observation provides strong support for the hypothesis presented above, namely, that the (OTf ) anion is adsorbed in the SURMOFs.Upon immersing the La(OTf ) 3 @Cu 3 (BTC) 2 SURMOF in dichloromethane under ultrasonic conditions for 15 min, the XRD pattern has no significant change (Fig. S4, ESI †) suggesting an irreversible loading of the SURMOF with La(OTf ) 3 .

Modelling
Molecular modeling simulations were performed using a Monte Carlo algorithm implemented in the simulation package SIMONA. 30n order to avoid the trapping of the system in metastable conformations, the basin hopping approach [31][32][33] with consecutive simulated annealing cycles was used to obtain the lowest energy for the system.Each simulated annealing cycle consists of a series of Metropolis Monte Carlo steps starting at high temperature, where the molecule can escape local minima and traverse energy barriers.Throughout the simulation the temperature is lowered according to: where T n is the temperature at step n, T E is the end temperature and N is the total number of steps in simulated annealing run.The interaction between the atoms was modeled by the classical electrostatic potential: and the Lennard-Jones interaction which includes terms accounting for van-der-Walls attraction and Pauli repulsion: We used standard parameters of the Lennard-Jones interaction for all atoms except La.In order to determine the LJ parameters for La 3+ we performed Metropolis Monte Carlo simulations of La 3+ with 3(OTf ) at constant low temperature (T = 1 K) and varied the La-parameters until the distance between oxygen and La 3+ reached the equilibrium distance obtained from ab initio calculations.The resulting parameters for La 3+ are e = 0.15 kcal mol À1 and s = 3.425 Å. Partial charges were calculated using the ESP fit implemented in the DFT code TURBOMOLE 34 using SV(P) basis and b3-lyp functional (see ESI † for details, Fig. S5).
To screen possible configurations a system consisting of two neighboring cavities (the medium and the large cavity of the MOF) was constructed.Each cavity was loaded with up to four La(OTf ) 3 molecules resulting in 24 different model systems.
To prevent the ions from leaving the complex we applied a strong restraining force to particles outside the cavities: U confine (q) = k(q À x 0 ) 2 , q 4 x 0 or q o Àx 0 , q = x, y, z The simulations consisted of 10 simulated annealing runs with N = 100 000 Metropolis Monte-Carlo steps each, starting at T S = 750 K and ending at T E = 300 K using rigid body translation for La 3+ and (OTf ) and in addition rigid body rotation and rotation around dihedral angles in the (OTf ) ions.The step size for the random rigid translations was drawn from a random uniform distribution between 0 and 1 Å, the angle of the rigid rotations from a random uniform distribution between 01 and 1801.

Experiments
The out-of-plane XRD data for the pristine HKUST-1 SURMOFs are shown in Fig. 1c (black).The (002) and (004) peaks are welldefined and sharp, their relative intensity agrees with that observed for bulk MOF powders obtained from simulations using the well-established bulk structure of HKSUST-1 (Fig. 1a).The absence of diffraction peaks for other crystallographic directions in the out-of-plane data reveals that that the SURMOF growth proceeds only along the (001) direction on the MHDA SAM.The in-plane XRD data showing only the (200) and (220) direction also support this orientation structure (Fig. 1c, red).XRD data recorded after immersion into the La(OTf ) 3 solution are shown in Fig. 1c (out-of-plane, blue; in-plane, dark yellow).No new diffraction peaks appeared revealing that the structure of the MOF lattice has remained unchanged.However, the form factor is strongly different, after loading the ratio of (002)/(004) (out-of-plane data) has decreased from 1.1 for the pristine film to 0.7 for the La loaded film.Also the in-plane XRD data reveal pronounced changes: the ratio of (200)/(400) has dropped from 1.6 to 0.1.Longer immersion times did not lead to a further change of the relative XRD peak intensities.
The same loading process has also been carried out for powder Cu 3 (BTC) 2 MOF.The changes in relative peak intensities observed for the powders were similar to that seen for the SURMOFs (Fig. S6, ESI †).The ratio of (002)/(004) has dropped from 2.15 to 1.21 after loading La(OTf ) 3 .
This significant change in relative intensities of XRD-peaks directly demonstrates that La 3+ ions are loaded in virtually every pore of the SURMOF (and MOF powders), a decoration of the outer surface only would not be compatible with the change of diffraction peak relative intensities.
The results of the XPS analysis are shown in Fig. 2. Clearly, after loading the data recorded for the La(OTf ) 3 @Cu 3 (BTC) 2 SURMOFs show distinctive La peaks. 35A quantitative analysis reveals that the total copper and lanthanum content of the SURMOFs after loading are 4.3% and 0.8%, respectively (Table S2, ESI †).
The metal ion concentration curves from inductively coupled plasma optical emission spectrometry (ICP-OES) measurement are shown in Fig. 3.By calculation of the integrated area of the measured curves, the copper and lanthanum concentration after loading the sample are 6.6 mg L À1 and 0.9 mg L À1 , respectively (Table S3, ESI †), the calculated ratio of Cu/La amounts to 7.3 : 1, which is close to the value calculated from XPS data.

Simulations
In order to complement the experimental data to obtain detailed information regarding the positions of the molecules within the MOF pores we have carried out force-field based Monte Carlo simulations for different numbers of La(OTf ) 3 embedded in Cu 3 (BTC) 2 .First, simulations of up to four La(OTf ) 3 moieties within a single Cu 3 (BTC) 2 unit cell were performed and the binding energies were calculated.To confine the La(OTf ) 3 moieties to the cavity an additional harmonic penalty potential (confining potential) was introduced for particles outside the cavity.Results for binding energies and the confining harmonic potential of the lowest energy state are reported in Table 1.The binding energy decreases with an increase in the number of moieties for up to four molecules, as repulsive interactions between the molecules compete with attractive interactions with the MOF.We note that the confining potential at the energy minimum is increased significantly when more than three La(OTf ) 3 are loaded in each per cavity, which indicates that not more than 3 molecules can fit one pore.There is a substantial binding energy associated with the loading, which mainly results from the dispersive interaction of the La 3+ ion, which is moved close to the central phenyl ring of one of the BTC ligands at the energy minimum (Fig. S7, ESI †).Because this interaction may not be fully represented in the force field, we have performed additional DFT calculations to obtain the binding between La 3+ and the phenyl rings and obtained a binding energy of À7.84 eV which is approximately 300 times k B T at room temperature.
In order to compare with the experimental data both cavities were loaded with different numbers of La(OTf ) 3 and simulation with 10 SA-cycles with 100 000 MC steps each, starting at T = 750 K and cooling down to T = 300 K, was carried out.The optimized structure with possible locations and numbers (1, 2, 3, 4) of La(OTf ) 3 in the large pore of Cu 3 (BTC) 2 MOF is presented in Fig. 4a and Fig. S6 (ESI †).The simulated out-ofplane XRD data for a (001) growth direction for the optimized structures obtained from the calculations are shown in Fig. 4b.
The increase of the (004) peak relative to the (002) peak in the calculated XRD data when adding the first La ion close to the Cu 2+ ions can be rationalized by considering the fact that this structure can (to first approximation) be described as an F lattice, and the intensity of the (002) peaks has to increase because placing the La 3+ ion close to the Cu ions effectively increases the electron density around these Cu position.When additional La 3+ ions are placed more close to the middle of the major pore (as in the structures of the 2 and 3 ions by the structure optimizations) the F-lattice is effectively changed into an I-lattice (Fig. S6, ESI †).The extinction rules for I-lattices lead to a decrease of the (002) peak intensity, as is evident from the figure for the case of 2 and 3 La 3+ ions.Further increase of the La content in the unit cell leads to an increase of the intensity.This can be explained by the higher atomic scattering factor of the La and its position.The former Cu F-lattice transforms into a La F-lattice (or even P-lattice) with some admixtures of Cu in the I-lattice positions.The atomic scattering factor for Cu   This journal is © the Owner Societies 2014 cannot fully compensate the La contributions to the scattered intensity.The effect explaining the intensity decrease for the 1La and 2La loading state applies also for the higher loading states for La, hence, an initial Cu dominated structure is not transformed into a La dominated structure.This result fits very well with the experiment data.The results for the two and three ions are very similar in the out-of-plane data.By comparing with the results of structure optimizations, loading with a single La ion can be ruled out because for the poor agreement of the simulated with the experimental XRD data.On the other hand a loading with 3 ions per HKUST-1 unit cell can be ruled out because of the poor agreement with the ICP-OES data.We thus conclude that in the loaded structure 2La 3+ ions are hosted within the pore of Cu 3 (BTC) 2 .This finding is fully consistent with all available experimental and theoretical information.

Conclusions
In conclusion, La(OTf ) 3 was successfully loaded into Cu 3 (BTC) 2 SURMOF via a solution impregnation method as evidenced by XPS, ICP-OES and IR data.By additionally considering the results of theoretical structure optimizations using a refined force-field analysis, we find a loading with two La-ions per unit cell to give the best agreements.This structure is stabilized by pronounced dispersive interactions between the La ion and the phenyl ring of a BTC ligand.The resultant loading of La (and other noble metal) species into the porous solid is of pronounced interest with regard to using MOF thin films for applications.Loading with metal ions and subsequently using redox-chemistry to create new species inside the MOF have recently received substantial attention in the context of information storage. 19n future work we plan to study the electrochemical properties of the loaded SURMOFs, which also offer rich potential for catalysis applications.