Electronic Supporting Information (ESI) to the article NMR relaxation and modelling study of the dynamics of SF6 and Xe in porous organic cages

Note on cage nomenclature and chirality. Helicity, or axial chirality, is an intrinsic property of the organic cage CC3. In this study, a homochiral diamine, (R, R)-1,2-cyclohexanediamine was used in the cage synthesis. The cage used throughout the gas study was therefore the homochiral form, CC3-R. Identical results would, presumably, be obtained with the opposite enantiomer, CC3S. Gas adsorption measurements suggested that the racemic equivalent, CC3-(R, S), also has similar properties with respect to gas adsorption. This latter cage would be less expensive to produce on a large scale because it can be prepared from the racemic cyclohexanediamine. However, the chiral forms, CC3-R and CC3-S, can form large single crystals with enhanced long-range order with respect to the less expensive racemic cage, CC3-(R, S). This is because the racemic form is much less soluble than the homochiral form and precipitates rapidly, preventing the growth of large, highquality single crystals. For this reason, the homochiral form of CC3 was used for these fundamental studies, in particular because it allows direct analysis of a highly crystal structure with fewer defects.


Introduction
Porous materials are solids that exhibit networks of pores and cavities and that have a much larger internal surface area as compared to the external surface area. Porous solids are ubiquitous and have widespread applications, such as molecular separation, purification, catalysis, chromatography, fluid transport, energy storage etc. 1,2 They also provide promising opportunities to capture greenhouse gases by adsorption, 3 such as CO 2 , CH 4 , as well as noble gases 4,5 (Xe, Ar, Kr). Various extended porous networks have been developed for these applications, such as activated carbons, 6 zeolites, 7 metal organic frameworks (MOFs), 8 covalent organic frameworks (COFs) 9 and porous polymers. 10 Nuclear magnetic resonance (NMR) relaxation and diffusion experiments provide extraordinarily versatile information about dynamics of fluid molecules in porous media. 11,12 Combined with a molecular level modelling, a detailed picture about the gas and framework dynamics can be revealed. 13 The modelling provides deep explanations for the observations and, overall, the analysis may guide in further system design.
This study focuses on the dynamics of sulphur hexafluoride (SF 6 ) and Xe gas in porous organic cages. SF 6 is a gas that has long atmospheric life span. It absorbs infrared light that makes it a potent greenhouse gas. 14 SF 6 acts as an excellent insulator for electrical equipment, thermos-acoustic insulator for windows, contrast agent in medical applications, and etchant plasma in semi-conductor industry. 15,16 However, due to its high global warming potential-23 900 times higher than CO 2 over a hundred year span 17 -it should be handled careful and recovered and stored in industrial usage. 16 Many techniques, including pressure swing adsorption/desorption processes, have been proposed to separate SF 6 from other gas mixtures and prevent its release to the atmosphere but the process requires high size selectivity and adsorption capacity of adsorbents. 18 Several extended porous solids 19,20 have been utilised for the separation and adsorption of SF 6 . Moreover, cucurbituril in aqueous solution 21 has also been utilised for the reversible trapping of SF 6 .
Xenon is used in light appliances due to its glow effects and in medical imaging field. The inertness of Xe and its low abundance in air (0.087 ppm by volume) 22 makes it commercially expensive. Due to its inert nature, a precise size selectivity is required to capture xenon gas in porous materials. 23 Porous organic cages are discrete cage molecules having permanent internal voids accessed by well-defined windows that provide porosity for guest molecules. 24 These cages provide a way for the analytical separation of rare gases and chiral molecules by selective bindings. 25 It has been shown that one of the organic cage molecules, CC3, 26 exhibits extraordinary characteristics for the selective separation of rare gases and SF 6 over gas mixtures. 26,27 The tetrahedral CC3 cage is composed of imine bonds that connect the rigid aromatic rings to the flexible cyclohexane linkers (Fig. 1A). These cages pack into a crystalline structure to form an interconnected 3D pore structure accessible via windows of the cages (Fig. 1B). The size of cage cavities is 4.4 Å while the narrowest point between the cage and window cavities is 3.6 Å in diameter. Adsorbed xenon (with diameter of 4.3 Å) was shown to have fast exchange, on the order of 10 8 s À1 , between these cavities, 28 and even SF 6 with kinetic diameter of 5.5 Å can enter the cage through the window cavity because of thermal and vibrational motion of the porous organic crystals. 27 The fluorine atoms in SF 6 have an NMR active 19 F isotope with spin-1/2, high gyromagnetic ratio and 100% natural abundance. The presence of six magnetically equivalent 19 F atoms in each SF 6 molecule makes the spin density of the gas high. The stable isotopes of S have a spin zero and natural abundance of 99.2%. Chemically SF 6 is inert, non-toxic, non-explosive, incombustible and thermally stable, which makes it a suitable NMR probe in various applications in chemistry, materials science and medicine. 29 It has been utilized for the investigation of porous media, 30 MRI of lungs, 31 liquid crystals, 32 and ethylene-propylene rubbers. 29 In addition, it has been used also as a spy molecule for highlighting the cavities in proteins. 33 Xenon has an NMR active 129 Xe isotope with spin-1/2 nucleus and relatively high natural abundance (26%). It is an inert, monoatomic noble gas with an easily polarizable electron cloud, which renders its chemical shifts extremely sensitive to the local environments. The sensitivity of the adsorbed 129 Xe can be enhanced by several orders of magnitude by applying a spinexchange optical pumping (SEOP) hyperpolarization technique. Therefore, it is widely used in NMR and magnetic resonance imaging (MRI) applications in chemistry, biochemistry, materials science and medicine. 34 It has been used in MRI of lungs, 35 microfluidic flow imaging, 36 investigation of liquid crystals, 37 polymers, 38 cements and shales 39 and ionic liquids, 40 for the determination of pore sizes of porous networks 41,42 and also as an biosensors. 43,44 In our previous work, 28 we performed a thorough 129 Xe NMR analysis of adsorption and dynamical phenomena of xenon in CC3. We calculated chemical shifts of 129 Xe in cage and window cavities, and, based on a simple two-site exchange model, we were able to interpret observed exchange averaged chemical shifts and T 2 relaxation times as well as extract populations of cage and window cavities and exchange rates.
In the present work, we extend the previous work by investigating T 1 , T 2 and diffusion of both Xe and SF 6 in homochiral CC3-R material (see Note on cage nomenclature and chirality in ESI †). Quantitative insights on the gas-mobility in CC3-R are deduced by analysing the data with a relaxation and diffusion model. We perform a global analysis of T 1 , T 2 and diffusion data using an extensive set of relaxation mechanisms for 129 Xe and SF 6 . The model benefits from a reduced number of adjustable parameters by incorporating quantum chemical information and configurational sampling. Three main aspects are covered here: first, extraction detailed information on the dynamics of the xenon and SF 6 gases in CC3-R; second, exploration of a model development where a reduced set of adjustable parameters are required by making use of potential energy surface (PES) to precompute relaxation model parameters; finally, we show how global analysis in the parameter estimation is performed with Bayesian Monte Carlo approach. 45 2 Materials and methods

Samples
Two CC3-R samples with different SF 6 loadings were prepared. The sample was inserted into a 5 mm glass tube and SF 6 gas was condensed in it by using liquid nitrogen. Thereafter, the tube was sealed with flame. The SF 6 : CC3-R molar ratio for the low loading (LL) and middle loading (ML) samples were 0.18 : 1, 0.53 : 1, respectively. Three CC3-R samples with different xenon loadings were prepared as well in the similar way as the SF 6 samples by using enriched 129 Xe gas (91%). The Xe : CC3-R molar ratios were 0.10 : 1, 0.52 : 1 and 2.4 : 1 for the low loading (LL), middle loading (ML) and high loading (HL) samples, respectively. The last sample represents an almost fully saturated material, in which nearly all three binding sites (one in each cage cavity plus four shared (equals to two) in the surrounding window cavities) are occupied by xenon. The adsorption rate of Xe sample was remarkably higher as compared to larger SF 6 molecules. 19   of accumulated scans were 60, 55 and 150 s, and 16, 1 and 8 for the LL, ML and HL samples, respectively. Diffusion delay D was varied from 0.05 to 25 s and length of the gradient pulsed from 0.3 to 2 ms.

Xenon
An introduction of 129 Xe relaxation mechanisms is provided, following the scheme in Fig. 2A. The relaxation contributions may be discussed in terms of amplitudes governed by the spininteraction strengths of the mechanisms (amplitude) and modulation, i.e. time constants of molecular processes causing the relaxation, summarized in Table 1. The details of the computation of relaxation rate amplitudes are described in Section 3.4. The ambition is to find as complete relaxation representation as possible with limited number of adjustable parameters that explains the experimental features in Fig. 3C and D.
The diffusion of gas molecules is a molecular process causing intermolecular dipole-dipole relaxation (DDinter) due to the interaction with the protons of the organic structure. This is illustrated with motion between unit cells in Fig. 2A, lower panel. The DDinter 129 Xe-H relaxation model was originally derived for liquid state relaxation. 46,47 In the implementation, the protons of the CC3-R structure are assumed to be fixed in space with the proton spin density of N H (computed in ESI †). The model assumes hard sphere diffusion with the spin located at the centre of a sphere, which is reasonable for 129 Xe. The amplitude of the DDinter interaction (see Table 1) depends on the length of closest distance (d HXe ) between 1 H and 129 Xe spins (in the liquid state derivation), 47 which can be explicitly computed by eqn (4). The DDinter mechanism is dependent on diffusion correlation time t D = d 2 HXe /D, where D is the experimentally known 129 Xe translational diffusion constant (see Fig. 3D).
The chemical shift difference of 129 Xe in the cage and window cavities (Dd = d C À d W ) 28 provides the shift-ex mechanism 48 that has a dominant role in T 2 relaxation. The amplitude (Table 1) is determined by Dd, the cavity molar fractions of 129 Xe (X S , S = C, W) and the strength of external magnetic field (B 0 ). The relaxation is modulated by the exchange rate (k ex ).
There are two local dipole-dipole mechanisms. The first one, DDHXe-ex, is due to exchange process and dependent on the b ex HXe interaction (see eqn (7)). The second one (DDHXe-L) is modulated by characteristic time t L due to local cavity motion (see Fig. 2A) and is dependent on the in-cavity dipole-dipole interaction (b S HXe , S = C, W; the subscripts refer to the 129 Xe-1 H pairs in window and cage cavities, see eqn (7)).
The CSA-L mechanism is due to the same local cavity motion (t L ) and have amplitude dependent on the chemical shielding anisotropy (Ds S ), populations (X S ) and B 0 .
The relaxation mechanisms are assumed additive (statistically uncorrelated) such that cavity specific relaxation rates (T i À1 ) m m = {DDHXe-L, CSA-L}, (i = {1, 2}) are population weighted: The cage-window exchange process is fast compared to the NMR timescale (milliseconds), based on the spectral observations in Section 4.1. These approximations are discussed further in Section 3.4.1. The mechanisms related to exchange (shift-ex and DDHXe-ex) are weighted with the product X C X W (cf. eqn (2) and (6)).

SF 6
The relaxation model for SF 6 is illustrated in the scheme in Fig. 2B and the symbols for amplitudes and modulation are shown in Table 2. The intermolecular mechanism DDinter contains the same parameters as in the case of xenon. Since the diffusion measurements were performed at higher temperatures than the relaxation experiments, one Arrhenius model was used for modelling diffusion coefficients in the entire temperature region of the diffusion and relaxation experiments. Three dipole-dipole mechanisms are included, DDFH-ex, DDFH-R and DDFF-R, where (b ex HF ), F-H (b C HF ) and F-F (b C FF ) couplings are considered, respectively. The first accounts for the two-site exchange modulation. The following two mechanisms are local and assumed to be modulated by the rotational motion with correlation time (t R ) in the cage cavity (C), weighted by the molar fraction X C . Simulation of SF 6 metadynamics in cage structure suggest there is a significant energy barrier and a rotation in passing the window cavity. 27 The barrier translates to a low window population and motivates not to include a specific window site relaxation process. The passing through the window cavity is accounted for locally with DDFH-ex and via the proton concentration dependent diffusion process DDinter. The local DD amplitudes have an additional parameter scaleR as compared to 129 Xe. Since there is only QC equilibrium and not temperature dependent trajectory information simulated, this parameter serves to rescale F-H and F-F distances in order to resemble the configuration averaged values (see ''r'' in eqn (7)).
The relatively complex local DDFH interaction is thus modulated by two correlation times due to rotation and site exchange. In the analysis (Section 4.3) the latter is approximated with t D . An additional fast (sub-picoseconds) liberation process can be anticipated, however, is not explicitly included other than in the final DD amplitude. For this process, an effect on the T 1 model cannot be ruled out, but a negligible influence on T 2 is estimated due to the short timescale.
The chemical shielding anisotropy (CSA-R) with respect to zz axis (along a bond) is large for SF 6 as shown by the quantum chemical (QC) calculations (see Section 3.4.2) and the relaxation is modulated by the correlation time (t R ) in the cage cavity. For the window cavity no large amplitude modulation is expected, and no CSA mechanism is included. Potential DD-CSA cross correlations are not included in the model as discussed in the ESI † (see Section S2.1).
For similar reason, nuclear spin-rotation (SR) mechanism is estimated for SF 6 rotation in the cage cavity following the Hubbard SR model. 49 In this case, the modulation is dependent on the inverse rotational correlation time and thus has a reverse temperature dependence. The SR tensors are estimated by using QC calculations (Section 3.4.2).
Strictly speaking, contrary to Xe, SF 6 molecule is not a perfect sphere, as the model assumes, and this has some effects on its dynamics in CC3. For example, it has been shown by simulations that, when SF 6 exits the cage, it has to align its outermost triangular face with the triangular face of the window. Thereafter a rotation of about 601 is required so that the second face can align the cage window and pass through. 27 However, the details of the dynamics in the window cavity is not really important in the relaxation modelling, as the population of SF 6 in the window cavity is very low. Furthermore, the transport of SF 6 enters the relaxation model via the diffusion coefficients. Thus, the details of SF 6 deformations are not explicitly required. Therefore, the spherical approximation should be well justified.

Spin-system dynamics
The relaxation models for 129 Xe and SF 6 consider a system with single spin in case of CSA, shift-ex and SR and a two spin system for dipole-dipole relaxation, in the standard perturbative form. 48 For completeness the 1/T 1 and 1/T 2 components are provided by eqn (S1)-(S10) in ESI. † Particular important for 129 Xe 1/T 2 is the shift-ex mechanism (eqn (S9) in ESI †): with amplitude parameters (see Table 1) and the spectral density J ex (0) at zero frequency. The latter depends on correlation time t ex = 1/k ex and in general resonance frequency o A for spin A. The spectral densities for local DD, CSA and shift-ex are computed as a Lorentzian from the correlation times t Y , Y = L, R, ex: The DDinter, 47,48 (see page 61 of ref. 48 and ESI †) has the non-Lorentzian spectral density associated with translational diffusion. The spectral density for SR is in the extreme narrowing regime, 49 and thus independent of resonance field. 129 Xe model includes a set of four fitted parameters (x) associated with two Arrhenius models for k ex and t L , respectively: x Xe = {t À1 ex0 , E ex , t À1 L0 , E L }. SF 6 study has seven fitted parameters x SF 6

Computational methods
The computational methods involving QC and PES are outlined together with parameter estimation procedure. The QC and PES computed system properties are directly related to the 129 Xe and SF 6 relaxation models discussed in previous section.  In this work we extend the use of the previously developed PES 28 to reduce the number of adjustable parameters in the relaxation model. The PES allows for sampling of spatial xenon configurations at the chosen temperature in the cage and windows cavities. Furthermore, a property surface of non-relativistic (NR) nuclear shielding tensor was developed, 28 from which the evaluation of shielding tensor at the spatial location is done at low computational cost. As explained in detail in ref. 28 a static reference nuclear shielding tensor was computed for xenon at the centre of both cage and window cavity (see Fig. 1B) including spin-orbit (SO) and scalar relativistic corrections at SO-ZORA/DFT(BHandHLYP) level.
The list of 129 Xe relaxation mechanisms is provided in Table 1. In this section the amplitude parametrisation is considered. The computed relaxation amplitudes follow from the hypothesis that there is statistically independent cavity dipole-dipole interaction and exchange and diffusion processes. Reasons for this separation may originate from a separation of timescale as well as approximately orthogonal processes. The latter processes are thus sensitive to averaged interaction tensors. Furthermore, our previous work 28 provides the mole fractions of Xe in the cage and window cavities (X C , X W ) that are here used as fixed parameters. Thus, the current relaxation model provides a validation of the two-site model used in our previous work. 28 In order to compute the estimates of shielding tensors, the Monte Carlo NVT sampling 28 was performed at 298 K. The calculations provided spatial Xe configurations and shielding trajectories in the cage and window cavities (not including the loading effect of the neighbouring cavities). The trajectories allow for relaxation modelling with molecular details but only with few adjustable parameters.
For DDinter, the average (denoted with hÁi i ) of the proton (i) to Xe distance at the equilibrium is computed as: where d S,i (S = C, W) are the H-Xe distances, and X S are the molar fractions known from experiments. 28 The inverse power of 6 weights the distances relevant to DD contribution. Equilibrium, as opposed to sampled Xe-H distances, are motivated by the fact that the local dipole-dipole interaction is accounted for with a separate mechanism. Eqn (4) gives d HXe : 5.52, 5.68, 5.72 Å, for LL, ML and HL samples, respectively. The shift-ex mechanism has a chemical shift difference sampled with the configuration trajectory as: with REL terms denote relativistic effect at SO-ZORA level (SO-ZORA-NR) at the static reference geometries. 28 Here, it is added on top of the sampled nonrelativistic (NR) chemical shift difference. The relaxation rate can be calculated after mean square averaging of Dd (see eqn (2)). The averaging leads to Dd = 181.9 ppm (the latter static relativistic contribution being 32.2 ppm). According to the simulation, the 129 Xe chemical shifts both in cage and window cavity increase slightly with temperature. 28 However, the difference changes only by 1 ppm over the relevant 80 K temperature interval and could affect relaxation amplitude (and rate) less than 1.5%. Thus, Dd evaluated at 298 K by eqn (5) was used at all the temperatures. To arrive at local DDHXe-ex amplitude for the cage-window exchange process, the starting point is the standard deviation of [(b C HA ) 2 À (b W HA ) 2 ] and, for comparison, computing the ensemble average as a population weighted effective distance: where SD denotes sampled ( j) standard deviation over equilibrium Xe-H i distances (i). The result of eqn (6) is r ave : 5.40, 5.36, 5.51 Å for the LL, ML and HL samples, respectively. The interaction (with r = r ave ) is given by: where m 0 is the permeability of vacuum, g A is the gyromagnetic ratio of nuclei A and h is Planck constant in angular units. The local cavity interaction (b S HXe , S = C, W) is provided in similar way from the effective distance: computed with SD for each sampled Xe i -H j distance. It is noted that population weight is taken in eqn (1) for individual cavity contributions. The results are r C : 3.05 Å, r W : 2.46 Å and they take the place of r in eqn (7). The CSA-L is obtained from shielding tensors trajectory (see page 106 of ref. 48) for the cage and window cavity in crystal frame and compute: resulting in Ds S : 12.7 and 21.3 ppm for S = W and C, respectively. 3.4.2 SF 6 quantum chemistry. The QC calculations to optimize SF 6 structure inside the cage/window cavities were carried out using Turbomole 50 code at non-relativistic density functional theory (DFT) level with hybrid B3LYP 51-53 functional. The DFT-D3 dispersion correction 54 was used in order to realistically describe the attraction between F atoms of SF 6 and H atoms in surrounding cage/window cavity structure. The def2-TZVP 55 basis set was used for the atoms of SF 6 molecule, while for the atoms in the surrounding cage/window cavities we used the smaller def2-SVP. 55 Only the atomic positions of SF 6 were optimized, while the coordinates of cage/window atoms were fixed to the DFT/PBE-TS optimized crystal structure (see details and structures in ESI †) that was used already in Xe NMR modeling of ref. 28.
The nuclear shielding tensors calculations with the Turbomole code were carried out using the B3LYP functional. For S atom in SF 6 molecule we used def2-TZVP basis set but for F atoms a completeness-optimized co-r primitive basis set (18s10p6d2f) was used, which provides result for all NMR parameters practically at complete basis set limit for all NMR parameters. 56 The def2-SVP basis set was used for the atoms in the surrounding cage/window structure. The obtained shielding tensors were symmetrized and diagonalized with python code to yield the isotropic, anisotropic and asymmetric parts in principal axis system (PAS) of the tensor for the relaxation models. The computed chemical shielding anisotropy (CSA) in the cage cavity was Ds C = 390 ppm and the chemical shift difference between the cage and window cavities was Dd = 10.7 ppm (see Table S5 in ESI †).
The nuclear spin-rotation (SR) constant (mean of tensor trace) was computed with DFT vacuum optimized geometry (for details see Section S4 in ESI †). The purpose of SR tensor is to explore the order of magnitude contribution of SR-relaxation of SF 6 in the cage.
3.4.3 SF 6 relaxation amplitudes. For SF 6 , an optimized structure in CC3-R was computed, allowing for some information on amplitudes of relaxation contributions (see Table 2).
The relaxation due to F-H dipole-dipole interaction in the cage cavity has a modulation due to overall rotation of SF 6 with correlation time t R . The amplitude (b C HF ) 2 may be computed following eqn (7) and (8), however, with the SD estimation from the fluoride equilibrium positions as opposed to NVT configurations. The effective distance was found to be r C : 2.72 Å. For the F-F interaction (b 2 FF ) the sum over the five intramolecular 1/d 6 FF distances were computed leading to an effective distance r FF = 1.8 Å, to be compared with the shortest equilibrium F-F distance of 2.2 Å. Thus, the sum of F-F relaxation channels provides an efficient relaxation mechanism. Due to the approximation of equilibrium distances these distances were rescaled with a fitted parameter scaleR. Amplitude for DDFH-ex is obtained from eqn (6), with the mean computed over F-distances. Related to nuclear spin-rotation (SR) and CSA the QC values (Section 3.4.2) where used.
3.4.4 Parameter estimation by MCMC. First, we summarize parameter estimation with Markov-chain Monte Carlo 57 (MCMC) and then discuss the relaxation model results for 129 Xe and SF 6 in Section 4. The deviation from experimental data is defined as the w 2 : where N exp is the number of experimental data, (s 2 i ) is the uncertainty (variance) of an experimental observation, i, estimated for the experimental data (see Section S4 in ESI †). The model is parametrized with the sets x Xe and x SF 6 defined in Section 3.3. For xenon the w 2 is computed from the {1/T 1 , 1/T 2 } data with diffusion constant interpolated from measurements. In the SF 6 study, diffusion experiments were recorded at higher temperature and w 2 was computed from the whole set {1/T 1 , 1/T 2 , Diff}, and diffusion was estimated with the Arrhenius model in the set x SF 6 . Thus, one activation energy is estimated based on data of the whole temperature interval.
The parameters in x are estimated together with confidence intervals by sampling the so-called posterior probability of the hypothesis x. The form of this distribution is given in Section S3 in ESI. † It consists of the likelihood function, i.e. the negative exponential of eqn (10) constrained by the prior assumptions of x. In practise, the distribution of x is obtained by generating MCMC-trajectories for the four and seven-dimensional parameter space for 129 Xe and SF 6 , respectively.
First the parameters are assumed to be of uniform distribution prior to the study, with their boundaries (cf. Tables S1 and S2 in ESI †), chosen based on physical constraints, for instance shortest plausible atom distance or alternatively, if the computational cost is too large, the boundary domain is reduced. The MCMC sampling was done following the Metropolis-Hastings algorithm, 45,57 where initial configurations are generated uniformly and kept within the boundaries (cf. Tables S1 and S2 in ESI †) by implementing reflective boundary condition. After calibrating MCMC step length such that approximately 50% of steps are accepted and simulating an initial ''burn in'' period the productive estimation follows. For Xe and SF 6 , 25 and 5 independent trajectories were simulated, with parameters initiated from uniformly distributed random numbers within the boundaries. In order to reduce chance of getting stuck in local minima during the ''burn in'' sequence an inverse ''temperature'' (w 2 -bÁw 2 , b = {10 À7 , 10 À3 , 1}) was used. Hence, first an equilibration at ''high temperature'' (small b) is performed, followed by the sampling of the distribution related to eqn (10) (i.e. b = 1).
Consensus is confirmed between the trajectories (with different starting parameters) in that the distribution of w 2 is similar for each trajectory as well as the trajectories predicts similar parameter distribution. Thus, from multiple trajectories it is seen that the final w 2 -distribution does not depend on the initial condition. 45 Production run of the trajectories are performed with 10 6 MCMC steps, as deemed necessary considering the parameter selfcorrelation in the MCMC-trajectory. The statistically best estimates of parameters are given by the mean over MCMC configurations. Their corresponding 95% confidence interval are estimated by numerically integrating each parameter cumulative probability distribution. Scripts are implemented in MATLAB.

Overview of NMR observations
19 F and 129 Xe NMR spectra of CC3-R with different SF 6 and xenon loadings measured at 295 and 260 K at 14.1 T are shown in Fig. 3A and B. The spectra show only one, relatively narrow peak due to fast thermal motion; the random changes in the 19 F and 129 Xe frequencies due to exchange between cage and window cavity, modulation of dipole-dipole, chemical shift anisotropy interactions, etc. occur much faster than the millisecond NMRtimescale. The chemical shift of 129 Xe increases significantly with loading, mostly due to increasing relative occupancy of window cavities, as the chemical shift of 129 Xe in the window cavity is much larger (about 210 ppm) than that in the cage cavity (about 20 ppm). 28 In the LL and ML samples, the 129 Xe chemical shift increases with temperature due to increased relative occupancy of window cavities, while the shift remains constant in the HL sample. 28 Contrary to the xenon samples, the 19 F chemical shift of SF 6 does not change almost at all with temperature or loading (see also Fig. S1A in ESI †). According to the quantum chemical calculations described in Section 3.4.2, the 19 F chemical shift difference between the cage and window cavities (about 11 ppm, see also Table S5 in ESI †) is much smaller than in the case of 129 Xe (about 190 ppm). Furthermore, SF 6 (5.5 Å) is larger than Xe (4.3 Å), and therefore it strongly favors the larger cage cavity over the window cavity, and changes in the temperature or loading does not change much the relative occupancies of SF 6 in the cavities. Consequently, 19 F chemical shift remains almost constant.
The transverse relaxation rate of both 19 F and 129 Xe decreases with increasing temperature due to the increased mobility (see Fig. 3C). T 2 of 129 Xe is also strongly dependent on the magnetic field strength whereas T 2 of 19 F is not. The predominant mechanism of T 2 relaxation of 129 Xe is the changing chemical shift due to the exchange between the cage and window cavities, 28 while in the case of 19 F it is intermolecular dipolar mechanism. The former depends on the square of the field strength (cf. eqn (2)), while the latter has vanishing field dependence at the slow 19 F mobility and high fields of this study (see eqn (S2) in ESI †). The T 2 of 19 F is approximately equal for the LL and ML loading (see Fig. S1C in ESI †) suggesting the same homogeneous broadening in their spectra. However, the spectral line width in Fig. 3A is almost double for the ML sample. Considering the few ppm larger chemical shift for gas phase SF 6 (QC estimate Table S6 in ESI †), this is consistent with a inhomogeneously broadened ML spectrum due to unresolved gas phase. The longitudinal SF 6 relaxation rate goes down by a factor of 5 with increased loading (see Fig. S1B in ESI †), most probably due to the same gas pool. Thus, the T 1 values, chemical shifts and linewidths are implying that not all the SF 6 molecules are absorbed in CC3-R (see also Section S1.1 in ESI †). As it is an additional complication to consider the effect of non-adsorbed gas, the quantitative modeling was conducted only for the LL SF 6 sample. Fig. 3D shows the self-diffusion coefficients of SF 6 and xenon in CC3-R. The diffusion coefficient of SF 6 is very low, about 4.8 Â 10 À13 m 2 s À1 at 335 K, and it is almost the same for the LL and ML samples (see also Fig. S2 in ESI †). This reflects very slow diffusional motion of the large SF 6 molecule in the small cavities of CC3-R. The diffusion coefficient of xenon in the HL sample is more than an order of magnitude higher than the diffusion coefficient of SF 6 because of the smaller size of xenon. On the other hand, the diffusion of xenon in the LL and ML samples is more than an order of magnitude larger than in the HL loading sample, because diffusion of xenon is significantly restricted in the almost saturated HL sample as xenon atoms cannot pass each other in the channels. Fig. 4 provides the experimental 129 Xe 1/T 2 and 1/T 1 relaxation rates in left and right columns, respectively, together with the model estimates. Panels A-B, C-D and E-F display the results of HL, ML and LL samples, respectively. The dominant mechanism for T 2 relaxation is the shift-ex mechanism for all the three samples, and the contributions of dipolar mechanisms (DDinter, DDHXe-L and DDHXe-ex) as well as CSA-L are negligible. For the DDinter mechanism, the diffusion constants were interpolated from the experimental data (see Fig. S6 in ESI †). Thus, the DDinter mechanism does not have any adjustable parameters. This mechanism has negligible T 2 influence due the fast diffusion. The most significant mechanisms to explain 1/T 1 are the DDinter, DDHXe-L and CSA-L mechanisms. Table 3 shows the resulting parameters of the Arrhenius model in bold (columns 4-7) with 95% confidence interval in brackets. The values of the cavity correlation time t L and exchange rate k ex at 298 K are given in columns 2 and 3 of Table 3. The exchange rates are in close proximity to our previous estimation, in which the chemical shift exchange mechanism was taken into account. 28 The exchange rate of the HL sample is 3-5 smaller than in the LL and ML samples, because almost all cage and window cavities are occupied by xenon atoms, restricting significantly the exchange between the cavities. Likewise, the cavity correlation time representing local dynamics show factor 70 slower dynamics at the HL sample than at the other samples. The activation energy is lower for the exchange process as compared to the local cage dynamics.

Xenon
The 1/T 1 increases with temperature for the HL sample and decreases for the ML and LL samples (see Fig. 4B, D and F). The higher temperature interval for HL cannot be the single explanation. It is found that both diffusion and local cavity dynamics are slower in the HL sample. It is noted that T 1 is particularly sensitive to dynamics close to resonance frequency (see eqn (3), with t L o 0 B 1). Both the net transport (t À1 D factor 3-5) and local cavity dynamics (cavity correlation time t L factor 70) is decreased by the high loading. The slow local dynamics at HL and the cause of 1/T 1 change of temperature dependence may be rationalized with high pressure increasing the influence of slow crystallite undulations on the shape of cavity. Hence, only at the timescale of these modes can the Xe gas fully explore the cavity, leading to the longer correlation time and increased activation energy. Table 3 shows transition from sub-ns to 25 ns, which is still too fast to provide a large change in 1/T 2 .
In Fig. 5A-C the exchange and diffusion mobility rates k ex and t À1 D of xenon are displayed versus temperature. Both have similar magnitude; however, the temperature dependence is less pronounced for k ex . This is also reflected in the activation energies, which are larger for diffusion than for exchange; the experimentally determined diffusion activation energies (see Fig. S6 in ESI †) are 33 AE 1, 11 AE 0.5 and 19 AE 5 kJ mol À1 for the HL, ML and LL samples, respectively, which are larger than the E ex values shown in Table 3.

SF 6
The experimental and modelled 1/T 2 and 1/T 1 values of SF 6 are shown in Fig. 6A and B, respectively. The corresponding individual mechanism are displayed in Fig. 6C and D. The 1/T 2 data are explained by the DDinter mechanism. For 1/T 1 DDFF-R together with DDFH-R and CSA-R play significant roles.
The estimate of SR and DDinter is shown to have negligible contribution. The diffusion data is displayed in Fig. 6E. Due to the very slow diffusion of SF 6 in CC3-R, we were able to measure diffusion coefficients only at high temperatures and the relaxation modeling extends the diffusion coefficient range to the lower temperature range. Diffusion is here modeled with one activation energy over the whole temperature range. Considering the ML T 1 experimental data (see Fig. S1B in ESI †) it is seen that a factor of 5 change is seen in increasing from LL to ML. This together change in lineshape is explained in terms of some free SF 6 gas present in ML sample (see Section 4.1).
In Table 4 the estimated parameters are provided with 95% confidence interval. The closest contact H-F spin distance is estimated to be 4.1 Å. This is relatively close to the static QC configuration estimate, 5.7 Å, following eqn (4). For the local DD mechanisms (DDFH and DDFF) the equilibrium distances were scaled down by 43% in the model; this is expected to be within a plausible range of distance deviations due to thermal motion. The SF 6 molecule has a rotational correlation time of Table 3 Resulting parameters of the Arrhenius model for xenon in CC3-R (bold) with 95% confidence intervals based on w 2 -MCMC simulation in brackets. The cavity correlation time t L (common for the cage and window cavities) and the exchange rate k ex at 298 K are given in columns 2 and 3 Gas load t L (ns) k ex (10 8 s À1 ) t À1 ex0 (10 8 s À1 ) E ex (kJ mol À1 ) t À1 L0 (10 10 s À1 ) E L (kJ mol À1 ) HL 25 [23,29]    6.7 ps and a diffusion constant of 1.08 Â 10 À13 m 2 s À1 at 298 K. These are computed from the estimated Arrhenius parameters (see Table 4). From Fig. 5D it is noted that SF 6 has more than two orders of magnitude slower diffusion as compared to the xenon samples. The model predicts 0.999 gas molar fraction of the cage cavity. The computed 19 F chemical shift difference between the cage and window cavities of 10.7 ppm together with the uneven populations gives a small shift-ex contribution (Fig. 6C). The Xe study suggest that exchange rate is in the regime of inverse diffusion correlation time (Fig. 5). Hence, k ex is approximated as t À1 D in the SF 6 shift-ex and DDFH-ex to keep adjustable parameters at a minimum.

Conclusions
In this work, adsorption and dynamics of SF 6 and xenon gases in CC3-R organic porous material was studied by 19 F and 129 Xe relaxation and diffusion NMR. The experiments show that the mobility of xenon decreases significantly close to the saturation of the cage and window cavity sites. Diffusion of SF 6 in CC3-R is two orders of magnitude slower than that of xenon due to its larger size. Much more detailed understanding about the dynamics and relaxation mechanism was achieved by building a common model for T 1 , T 2 and diffusion data. The model confirmed that cage to window exchange is a completely dominating mechanisms for 129 Xe T 2 relaxation (Fig. 4) and the exchange rate is about 10 8 s À1 (Fig. 5). The T 1 relaxation is dominated by the diffusion modulated dipole-dipole (DDinter) relaxation as well as the CSA and DDHXe-L relaxation due to the local cavity mobility (t L ). This is in contrast to the low loading SF 6 T 2 data, where the dominating mechanism is diffusion modulated dipole-dipole relaxation (Fig. 6); for T 1 the local tumbling of SF 6 in the cage cavity is the key dynamics entering the dipole-dipole (DDFF) and CSA mechanisms. The detailed information about gas transport in organic cages is expected to stimulate the design of these systems further for potential applications such as selective separation of noble gases and SF 6 . Table 4 Estimated parameters (bold) with 95% confidence interval (brackets) for SF 6 in CC3-R. The t R and D are illustrated at 298 K, extracted from the estimated activation energies, E R = 7.3 AE 0.07 and E D = 22.6 AE 0.2 kJ mol À1 , respectively