Adsorption and Diffusion in Siliceous Zeolites

1School of Chemical and Biomolecular Engineering, Georgia Institute of Technology, ... to represent the adsorbent structure and relies on electronic s...
11 downloads 0 Views 2MB Size
Subscriber access provided by Kaohsiung Medical University

C: Surfaces, Interfaces, Porous Materials, and Catalysis 4

First-Principles-Derived Force Fields for CH Adsorption and Diffusion in Siliceous Zeolites Hanjun Fang, Rohan Awati, Salah Boulfelfel, Peter I. Ravikovitch, and David S. Sholl J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.8b03267 • Publication Date (Web): 24 May 2018 Downloaded from http://pubs.acs.org on May 24, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

First-Principles-Derived Force Fields for CH4 Adsorption and Diffusion in Siliceous Zeolites Hanjun Fang,1§ Rohan Awati,1§ Salah E. Boulfelfel,1§ Peter I. Ravikovitch,2 and David S. Sholl1* 1

School of Chemical and Biomolecular Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0100, United States 2

Corporate Strategic Research, ExxonMobil Research and Engineering, 1545 Route 22 East, Annandale, New Jersey 08801, United States

Abstract Most previous studies on development of force fields for molecules in porous materials focus on prediction of adsorption properties. However, accurately reproducing adsorption data is not sufficient to guarantee the accuracy of other properties such as diffusivities. We demonstrate an approach to develop force fields based on periodic first-principles calculations that can accurately predict both adsorption and diffusion properties in crystalline nanoporous materials using CH4 in siliceous zeolites as an example. First, multiple dispersion-corrected density functional theory (DFT) methods were tested for describing CH4 in siliceous chabazite, and the two configurations (CH4 interacting with framework wall and 8-ring window) that are relevant to CH4 adsorption and diffusion were investigated. By comparing with the results from high level random phase approximation (RPA) results, DFT/CC (coupled cluster) was found to be the optimum method for force field development. The DFT/CC-derived force fields that accurately predict both adsorption and diffusion of CH4 in several commonly studied siliceous zeolites are then developed.

1 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 28

1. Introduction As an increasingly valuable complement to experiment, computational modeling has been used to interpret experimental observations and predict properties ahead of experiment for guest molecules in porous materials such as zeolites, metal-organic frameworks (MOFs), and carbon nanotubes. Firstprinciples quantum mechanical (QM) calculations are capable of accurately determining optimized geometries, binding energies, reaction energy profiles, and electron densities, etc. However, they cannot be readily used to predict macroscopic properties like adsorption isotherms, loading-dependent isosteric heats of adsorption, and molecular diffusivities. Instead, molecular simulations such as Grand Canonical Monte Carlo (GCMC) and Molecular Dynamics (MD) have been widely applied for predicting these properties.1,2 Successful implementation of GCMC and MD simulations requires accurate force fields (FFs) to describe the interactions of guest molecules with porous frameworks. Developing accurate and transferable FFs for adsorbed molecules is challenging. While generic FFs and experimentally-derived FFs work well for many simple systems, they often fail to describe interactions in more complex porous materials.3 Since first-principles QM approaches are capable of accurately predicting intermolecular interactions, deriving FFs from QM data without experimental input is a promising solution. Significant recent progress has occurred in deriving FFs from QM calculations.3-9 We previously used this methodology to develop FFs for adsorption of CO2 in siliceous and cationic zeolites.10-12 This approach uses the fully periodic framework to represent the adsorbent structure and relies on electronic structure calculations for hundreds of adsorption configurations randomly scattered throughout the framework. These FFs accurately predict experimentally observed adsorption properties and show good transferability across different zeolite topologies. In this work, we extend our approach to CH4 in siliceous zeolites. Previous FFs for CH4 and higher alkanes in zeolites have been developed by fitting adsorption isotherms and/or heats of adsorption to experimental data.13,14 However, accurately reproducing adsorption data is not sufficient to guarantee the accuracy of other properties. For example, Jee et al. found that the FFs derived from experimental isotherms for CH4 and CO2 in the silica zeolite DDR cannot accurately describe molecular diffusion in this material.15 This occurs because the adsorption isotherm is dominated by configurations with molecules close to the energy minima, while diffusion is controlled by transition states that would lie at considerably higher energies. Our goal in this work is to develop the first-principles-derived FFs that can accurately predict both adsorption and diffusion properties for CH4 in siliceous zeolites. In section 2, we assess the accuracy of various dispersion-corrected DFT methods by computing adsorption energies of CH4 on siliceous chabazite (Si-CHA), where high-level QM results are already available and can be used for comparison. 2 ACS Paragon Plus Environment

Page 3 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

On the basis of these benchmark calculations, we choose an optimum DFT method to develop FFs that account for CH4–zeolite interactions. To test the accuracy and transferability of the resulting FFs, in section 3 we perform GCMC and MD simulations to predict adsorption and diffusion properties of CH4 in multiple siliceous zeolites and compare the results with the experimental measurements. A comparison between our developed FFs and those used in literature for CH4 diffusion in the siliceous zeolites is also performed. Finally, our conclusions and plans for future work are discussed in section 4.

2. Methods 2.1 First-principles Methods for Describing CH4 in Zeolites The accuracy of a first-principles-derived FF is controlled by the accuracy of the QM method upon which it is based.3 Similar to our previous work of FF development for CO2 in zeolites,10 we first performed benchmark calculations using different first-principles QM methods for CH4 in a periodic zeolite structure, aiming to find a QM method that can accurately describe CH4–zeolite interactions with affordable computational cost. The performance of different QM methods is evaluated via the comparison with reliable experimental data or high-level QM results. We chose adsorption of CH4 in Si-CHA (R 3 m space group) as the benchmark system,16 because the computational cost for this material using its primitive unit cell is lower than that of many other zeolites. For this reason, CHA zeolites have been often used for the benchmarking of QM methods.17-23 Among these studies, Göltl et al. investigated the adsorption of CH4 and other small alkanes in purely siliceous and Na-exchanged CHA at different levels of theory, including the standard density functional PBE, dispersion-corrected PBE-D2, vdW-DF, and the random phase approximation (RPA) in combination with the adiabatic-connection fluctuation-dissipation theorem (RPA-ACFDT) using wavefunctions.17 They found that the PBE method underestimates the strength of the interactions between these small alkanes and zeolites, while PBE-D2 and vdW-DF strongly overestimate the adsorption energies. The most accurate description is achieved at the RPA level with HF exchange energies (RPAHF).17 Later, they found similar trends for the small alkanes in protonated CHA. 18-20 The overestimation of adsorption energies using PBE-D2 or -D3 and a variety of vdW-DF approaches has also been observed for CH4 in silicalite, HZSM-5, MOR, ZSM-12, and SAPO-34 zeolites.21,23-25 Although the RPA-HF method performs well for CH4 and other small alkanes in zeolites, the computational cost of this method is still very high for the purpose of FF development, where hundreds of periodic calculations need to be performed. Bludský, Nachtigall, and co-workers proposed the DFT/CC (density functional theory/coupled cluster) method for an accurate description of weakly bound molecular systems.26 This method is based on a correction to DFT, ∆EDFT/CC, defined as the difference between DFT and accurate CCSD(T) interaction energies. The DFT/CC method has been successfully applied to a 3 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 28

variety of systems where dispersion interactions are important including adsorption of various molecules on graphene and graphite surfaces,27,28 water in the MOF CuBTC,29 as well as CO2 in proton and alkaline cation exchanged zeolites.30-33 Chen et al. derived an ab initio potential energy surface (PES) for CH4 interacting with MOF CuBTC at the DFT/CC level, then implemented the PES in GCMC simulations, and obtained good agreement between calculations and experiments for adsorption isotherms and adsorption mechanism.34 Recently, Rubeš et al. used DFT/CC-based method for the test on isosteric heats of adsorption of CH4 in several siliceous zeolites, and they found that DFT/CC and DFT/CC-based FF yield accurate isosteric heats of adsorption for most of the zeolites they studied when compared with experimental results.35 In their work, however, prediction of adsorption isotherms and diffusivities was not included. We have previously applied the DFT/CC method in FF development for CO2 in siliceous and cationic zeolites, and the resulting FFs work well for these systems.11,12 In this work, we examine the performance of the DFT/CC method for CH4 in Si-CHA. Based on the work by Bludský et al.,26 we developed CC corrections for CH4 in siliceous zeolite and included the details in Supporting Information. We also evaluated some other recently proposed dispersion-corrected DFT methods such as the Tkatchenko-Scheffler

(DFT-TS)

method,36,37

the

many-body

dispersion

energy

method

(MBD@rsSCS),38,39 and the density-dependent energy correction method (DFT-dDsC)40,41. DFT-TS and DFT-dDsC are similar to DFT-D2, but their dispersion coefficients and damping function are chargedensity dependent, and therefore they account for the influence of the local chemical environment. The MBD@rsSCS method represents a conceptually clean approach for including infinite-order many-body effects in the dispersion energy calculations. We refer the readers to the references for the detailed information about these methods.

2.2 Benchmark Calculations for CH4 in Siliceous Chabazite (Si-CHA) It should be noted that the benchmark studies by Göltl et al. focus on the most energetically stable configurations for CH4 (and other small alkanes) in CHA zeolites, which are relevant for predicting adsorption properties of guest molecules in zeolites.17-20 For reliably predicting diffusion properties, the transition states (TS) through which the guest molecule diffuses from a cage to the neighboring cage or cross the channel also need to be taken into account in the benchmark calculations. For CH4 in Si-CHA, the most energetically stable state is CH4 interacting with the framework wall inside the super-cage,17 while the TS is where CH4 is located near the center of an 8-ring window.15,42 Figures 1a and 1b show these two configurations, respectively. We obtained the optimized geometries of these two states at different levels of theory, including PBE, PBE-D2, PBE-D3, PBE-TS/HI, PBE-MBD@rsSCS, PBE-dDsC, vdW-DF, and vdW-DF2, which have been implemented in the Vienna Ab initio simulation package 4 ACS Paragon Plus Environment

Page 5 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(VASP).43,44 Currently RPA-HF method cannot be used for calculating atomic forces, and therefore cannot readily be used for geometrical optimization. The RPA-HF energies were on the basis of the geometries optimized at other levels of theory. The same is true for the DFT/CC method. The results are summarized in Table 1, including the geometrical parameters and adsorption energies (∆Eads) at different levels of theory. The difference in energy between the two configurations (∆∆Eads) can be viewed as the enthalpic contribution to the diffusion barrier. The results from Göltl et al. for CH4 interacting with framework wall inside the super-cage are also included in Table 1 for validation purpose.17 Our results are generally consistent with theirs. Some deviation exists for both the nearest atomic distances and adsorption energies, and this is probably because of minor differences in the local minimum states in the two sets of calculations. The adsorption energies (∆Eads) from RPA-HF are similar among different optimized geometries, in the range of –12.1 to –14.7 kJ/mol for CH4 in the super-cage, and –0.2 to –2.2 kJ/mol in the 8-ring window. For the diffusion barriers (∆∆Eads), the difference is even smaller, within 1 kJ/mol. Previous studies indicate that caution must be used in comparing the calculated adsorption energies of guest molecule on zeolites with experimentally measured heats of adsorption.10,19,45 The QM calculations only consider the global or local minimum configurations of guest molecules, while experimental heats of adsorption are an average of all possible configurations of guest molecules in the pore, even at low loading, because of thermal effects. A more straightforward way to evaluate the accuracy of the DFT methods would be to compare the DFT results with those from higher-level QM methods. In this case, the results from the RPA-HF method, which has been demonstrated to reliably describe adsorption of CH4 and other short-chain alkanes in CHA zeolites, can be used as the reference values.17-19 As expected, PBE significantly underestimates adsorption energies for CH4 in both the supercage and 8-ring window. PBE-D2, PBE-D3, PBE-MBD@rsSCS, PBE-dDsC, vdW-DF, and vdW-DF2 overestimate adsorption energies for both configurations. PBE-TS/HI performs well for CH4 in the 8-ring window, but overestimates the adsorption energy for CH4 in the super-cage. It is interesting that the diffusion barriers predicted from PBE-D3, PBE-MBD@rsSCS, PBE-dDsC, and vdW-DF2 are close to those predicted from RPA-HF, with differences within 2~3 kJ/mol. DFT/CC gives adsorption energies in range of –16.4 to –17.2 kJ/mol for CH4 in the super-cage, and –4.3 to –5.7 kJ/mol for CH4 in the 8-ring window, which indicate that the DFT/CC energies are not sensitive to the methods used in geometry optimizations. Our results for CH4 in the super-cage are close to that from Rubeš et al, which is –16.8 kJ/mol at the DFT/CC level of theory for CH4 in Si-CHA.35 When compared with the RPA-HF results, DFT/CC overestimates adsorption energies by 3-4 kJ/mol for both CH4 configurations. The diffusion barriers predicted from DFT/CC, however, are almost identical to those from RPA-HF. These benchmark 5 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 28

calculations indicate that DFT/CC gives good agreement compared to RPA-HF. Considering the trade of between accuracy and computational cost, DFT/CC was used in following force field development.

Figure 1: Two adsorption configurations of CH4 in siliceous CHA: (a) in the super-cage, and (b) in the 8ring window. They are viewed along two directions. C, H, Si, and O atoms are depicted as gray, white, dark cyan, and red spheres, respectively. For clearer visualization, the atoms of CH4 and framework atoms that are close to CH4 are represented by ball-and-stick model, while the rest atoms by line model.

6 ACS Paragon Plus Environment

Page 7 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Table 1: Adsorption energies (in kJ/mol), ∆Eads, of one CH4 molecule in siliceous CHA calculated at different levels of theory, and the nearest distance (in Å), Rmin (C–Si), between the carbon atom of CH4 and silicon atoms of zeolite. The adsorption configurations are shown in Figure 1. The energy difference from the super-cage to the 8-ring, ∆∆Eads, is also listed. When a geometry was obtained from optimization at a different level of theory, the functional used for optimization is given in parentheses. Rmin (C–Si) ∆Eads ∆∆Eads Super-cage 8-Ring Super-cage 8-Ring PBE 4.44 (4.25)a 3.84 -2.9 (-2.4)a 21.2 24.1 a PBE-D2 3.99 (4.15) 3.79 -21.3 (-18.3)a -12.0 9.3 PBE-D3 4.38 3.81 -22.7 -8.7 14.0 PBE-TS/HI 4.21 3.81 -22.8 0.4 23.2 PBE-MBD@rsSCS 4.20 3.81 -22.1 -9.8 12.3 PBE-dDsC 4.21 3.79 -25.8 -12.4 13.4 vdW-DF 4.26 (4.27)a 3.86 -32.6 (-33.2)a -7.4 25.2 vdW-DF2 4.22 3.83 -23.9 -9.3 14.6 DFT/CC (PBE geo.) -16.4 -4.8 11.6 DFT/CC (PBE-D2 geo.) -17.0 -4.4 12.6 DFT/CC (PBE-D3 geo.) -16.6 -4.3 12.3 DFT/CC (vdW-DF geo.) -17.0 -4.8 12.3 DFT/CC (vdW-DF2 geo.) -17.2 -5.7 11.5 RPA-HF (PBE geo.) -12.1 (-14.0)a -0.2 11.9 RPA-HF (PBE-D2 geo.) -14.4 (-8.2)a -1.7 12.7 RPA-HF (PBE-D3 geo.) -14.3 -1.5 12.8 a RPA-HF (vdW-DF geo.) -14.7 (-14.5) -2.2 12.5 RPA-HF (vdW-DF2 geo.) -12.6 -0.6 12.0 a

the values in parentheses are from reference17

2.3. Force Field Model and Computational Method It is important to use a CH4–CH4 potential that correctly captures the phase behavior of pure CH4, so we used the TraPPE potential,46 which was developed for this purpose. In order to fully specify the energy of adsorbed molecules, we only need to define the CH4–zeolite interactions. In previous studies of CH4 in zeolites, FF parameters are often derived by fitting isotherms to experimental data, and CH4 molecule is treated as the united-atom (UA) model, where the carbon atom and its bonded hydrogen atoms are grouped to form one unit. In this case no Coulomb interactions are included between CH4 and the zeolite. It is also relatively common to include only vdW interactions with framework oxygen atoms, and to assume that interactions with silicon atoms are included through the effective interaction potential with O atoms.14,47 The UA model for CH4 would be a good approximation to simulate molecular systems in which the intermolecular motion is much more important than the intramolecular motion. Treating the vdW interaction from silicon atoms through the effective potential with oxygen atoms is reasonable since silicon atoms are screened by the bulkier oxygen atoms. To test these two approximations, we tested two kinds of FFs for CH4 in siliceous zeolites as illustrated in Figure 2. One uses a UA model for CH4 and includes only vdW interactions between CH4 7 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 28

and framework oxygen atoms (Figure 2a). The other uses a rigid all-atom (AA) model with atomic charges for each atom of CH4 molecule and includes vdW and Coulombic interactions with both silicon and oxygen atoms between CH4 and zeolite (Figure 2b). The geometry of CH4 in the AA model was optimized using PBE-D2.

Figure 2: Two models used in force field development for CH4 in siliceous zeolites: (a) the united-atom (UA) model without atomic charges is for CH4, with only oxygen framework atoms being explicitly included for vdW interactions between CH4 and zeolite, and (b) the all-atom (AA) model with atomic charges is for CH4, and both silicon and oxygen framework atoms are explicitly included for van der Waals (vdW) and Coulomb interactions between CH4 and zeolite

When the AA model is used, we assume the interactions between each atom in CH4 and a zeolite atom are represented by pairwise vdW and Coulombic interactions,

EFF ( Rij ) = EvdW + ECoul = s12

C6ij qi q j C12ij − s + 6 Rij12 Rij6 Rij

(1)

where Rij is the distance between atoms i and j, C6ij and C12ij are the attractive and repulsive coefficients, qi and qj represent the charges for atoms i and j, respectively, and s6 and s12 are global scaling factors that are fitted to allow the closest correspondence between the first-principles results and the classical FF in a least-squares sense. C6ij values were adopted from Grimme’s work,48 and C12ij were determined by a simple relation we used before10

C12ij ( R0i + R0j )6 = C6ij 2

(2)

8 ACS Paragon Plus Environment

Page 9 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

where R0i and R0j are the vdW radii of atoms i and j from Grimme et al.48 Table 2 gives the C12ij and C6ij coefficients for each cross-species interaction between CH4 and Si-CHA. The Coulomb interactions were treated by a point charge model employing density derived electrostatic and chemical (DDEC3) charges for the atoms of zeolite framework and CH4 molecule.49-51 When the UA model is used, the Coulomb term in eq. (1) is not included. In the FF fitting, the residual standard deviation (RSD) is minimized, where n

∑(E

k FF

k − EDFT / CC )

2

k

RSD =

(3)

n−2

k k Here EFF is the interaction energy calculated at the FF level for CH4 configuration k, EDFT / CC is the

interaction energy calculated at the DFT/CC level, and n is the number of configurations in FF fitting dataset. The mean deviation (MD) is also calculated after the parameterization, n

∑( E MD =

k FF

k − EDFT / CC )

k

n

(4)

The RSD and MD can give an overall evaluation of the performance of the fitted FF in reproducing the QM data. We used the same primitive Si-CHA unit cell as the benchmark section for force field development. The unit cell was first fully optimized at the PBE-D2 level, which gives lattice constants of a = b = c = 9.333 Å and α = β = γ = 94.34°, close to the experimentally determined data, a = b = c = 9.229 Å and α = β = γ = 94.3°.16 DFT/CC calculations were performed with a loading of one CH4 per unit cell. To generate CH4 configurations in FF fitting dataset, we applied NVT MC (N = 1, T = 300 K, and AA model for CH4) simulations based on the 3×3×3 extended supercell. We observed that the ordinary NVT MC technique only leads to energetically favorable configurations for CH4 in the zeolite framework that are relevant to prediction of adsorption properties. As mentioned above, the energetically unfavorable configurations where CH4 sits inside 8-ring windows are critical for describing diffusion behaviors. To effectively generate configurations of CH4 in 8-ring windows, we firstly found the center of each supercage within the supercell of Si-CHA, and during the NVT MC simulations we blocked these centers using spheres of radius of 4, 3, 2, 1, and 0 Å, respectively. This way helps to make CH4 towards 8-ring windows, and configurations of CH4 in both super-cages and 8-ring windows can be obtained (see Figure S3 in Supporting Information). For each blockage radius, the same number of CH4 configurations was randomly chosen from the NVT MC snapshots. We generated CH4 configurations and derived FF parameters through two iterations. In the first iteration, the ClayFF parameters52 were used for the zeolite atoms (O and Si) and the TraPPE parameters46 9 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 28

and DDEC3 charges for the atoms of CH4. Lorentz-Berthelot mixing rules were applied for the LennardJones for cross interacting species between CH4 and zeolite. A total of 300 configurations were chosen from the NVT MC snapshots and used in the first iteration of force field fitting. The fitted FF parameters were then used in in the second iteration following the same procedure. Then a total of 1,000 configurations were included to refine the force field parameters. Further iterations did not influence the final FF parameters. The final FF parameters using AA and UA models are summarized in Table 3. For simplicity, we will refer to the two FFs derived in this way as CCFF_AA and CCFF_UA. It is useful to note that the number of independent fitting parameters in the two FFs is the same (s6 and s12 in eq. (1)), despite the difference in the functional form of the FFs. We have tested another set of blocking spheres (4.0, 3.5, 3.0, 2.5, 2.0, 1.5, 1.0, 0.5, and 0.0 Å). The interval of adjacent spheres is smaller, and found the final fitted force field parameters change by within 1% and 2% for for UA and AA models, respectively, and these changes would not significantly influence the adsorption and diffusion properties. Table 2: Attractive and repulsive coefficients (C6ij and C12ij) and the sum of van der Waals radii (R0i + R0j) for each cross species between CH4 and siliceous zeolites. Cross species C12ij C6ij R 0i + R 0j 12 -1 6 -1 (Jnm mol ) (Jnm mol ) (Å) AA model UA model Si–C Si–H Oz–C Oz–H

Oz– CH4

2.031×10-3 2.287×10-4 2.637×10-4 2.590×10-5

Table 3: CCFF_AA and CCFF_UA parameters for CH4 in siliceous zeolites. Cross species ε/kB (K) σ (Å) CCFF_AA Si–C 33.98 3.845 Si–H 24.15 3.298 Oz–C 19.88 3.391 Oz–H 16.17 2.844 CCFF_UA Oz–CH4 109.76 3.416

4.019 1.137 1.107 0.313

3.168 2.717 2.794 2.343

Charge (e) Si (2.2124) Oz (-1.1052) C (-0.516) H (0.129)

2.3. Force Field Fitting A detailed comparison of the interaction energies predicted using CCFF_AA and CCFF_UA and the corresponding energies at the DFT/CC level is shown in Figure 3. When the AA model is used, the 1000 DFT/CC calculations span a range of adsorption energies, from -17 to 7 kJ/mol, and the CCFF_AA interaction energies are reasonably consistent with the DFT/CC results (Figures 3a and 3b). Some

10 ACS Paragon Plus Environment

Page 11 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

deviations between CCFF_AA and DFT/CC energies exist, especially when the CH4 molecule is located near the center of 8-ring windows of Si-CHA (Figure 3c). The CCFF_AA interaction energies for 99.5% (100%) of the CH4 configurations have deviations within ±4 (±8 kJ/mol) of the DFT/CC energies. The RSD and MD between the FF and the DFT/CC results are 0.9 and 0.04 kJ/mol for the AA model, respectively. When the UA model is used, a larger deviation between CCFF_UA and DFT/CC energies is observed than that from the AA model (Figure 3d). This is not surprising since stronger approximations are made in the UA model to describe CH4 than that of the AA model. For the most energetically favorable configurations (Figure 3e), CCFF_UA underestimates the interaction energies by about 2 kJ/mol compared to the DFT/CC results, and deviation increases to 7 or -8 kJ/mol when the CH4 molecule is located near the center of 8-ring windows (Figure 3f). The CCFF_UA interaction energies for 94.7% (99.8%) of the CH4 configurations have deviations within ±4 (±8 kJ/mol) of the DFT/CC energies. Overall, RSD and MD between the FF and the DFT/CC results are 1.4 and 0.01 kJ/mol for the UA model, respectively. We also checked CCFF_AA and CCFF_UA in calculating the adsorption energies of the two configurations of CH4 in Si-CHA as shown in the benchmark calculations (section 2.2). During the calculation, the position of CH4 was optimized, while the framework of Si-CHA was fixed to the PBE-D2 minimized geometry. The results are summarized in Table 4. Both the CCFF_AA and CCFF_UA perform well in predicting the adsorption energies and diffusion barriers compared with the results using DFT/CC (Table 1). Geometrical parameters are also predicted well using the two FFs when compared to the results at the different levels of DFT.

Table 4: Adsorption energies (in kJ/mol), ∆Eads, of one CH4 molecule in siliceous CHA calculated by the developed force fields, and the nearest distance (in Å), Rmin (C–Si), between the carbon atom of CH4 and silicon atoms of zeolite. The adsorption configurations are shown in Figure 1. The energy difference from the super-cage to the 8-ring, ∆∆Eads, is also listed. Rmin (C–Si) ∆Eads ∆∆Eads Super-cage 8-Ring Super-cage 8-Ring CCFF_AA 4.13 3.92 -17.4 -5.7 11.7 CCFF_UA 4.26 3.94 -15.6 -1.7 13.9

11 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47

Page 12 of 28

Figure 3: Force field fitting results for CH4 in siliceous CHA: (a, d) Comparison of the interaction energies of CH4 in siliceous CHA zeolite for the CCFF and DFT/CC, (b, e) the difference in interaction energies (ECCFF –EDFT/CC) as a function of EDFT/CC, (c, f) ECCFF –EDFT/CC as a function of the nearest interatomic distance between the carbon atom of CH4 and the center of 8-membered ring window. (a, b, and c) are for the all-atom (AA) model, and (d, e, and f) for the united-atom (UA) model.

12 ACS Paragon Plus Environment

Page 13 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

3. Results and Discussion 3.1 Prediction of Adsorption Properties To examine the applicability of the FFs developed above, we compared adsorption properties predicted using CCFF_AA and CCFF_UA with experimental adsorption data for several siliceous zeolites. Adsorption isotherms and isosteric heats of adsorption were predicted computationally using standard GCMC methods (see Supporting Information). We emphasize that our calculated results come from our first-principles derived FFs and are not fitted to experimental data in any way. We first discuss results for CH4 adsorption in siliceous CHA, the zeolite on which the FF development is based. The simulated and experimental isotherms are shown in Figure 4a. The simulation results using CCFF_AA and CCFF_UA at 301 K are close each other, and in good agreement with the experimental isotherms from three different groups.53-55 The isosteric heats of adsorption at zero loading are predicted to be 17.7 and 17.0 kJ/mol for CCFF_AA and CCFF_UA (Figure 4b), respectively, which are accordance with the experimental result from Pham and Lobo.55 Our simulation results are also close to the experimental result (16.8 kJ/mol) obtained by Rubeš et al.35 To examine the transferability of our force fields, we performed GCMC simulations of CH4 in several other zeolites. Figures 4c and 4d show the isotherms and heats of adsorption for ITQ-29 (siliceous LTA). The agreement between the simulations and experiments is good for isotherms in Si-LTA, with some slight deviations compared to the work from Hedin et al.53 Considering the variations between experimental results from different groups, the predictions based on our method are good. The heats of adsorption of CH4 in ITQ-29 at all loadings are similar to each other with the two FFs and they are generally lower by about 2 kJ/mol compared with the experimental data from Palomino et al.56 Figure 5 shows the results for CH4 in siliceous MFI, which is an extensively studied system in experiment. The agreement between the simulations and experiments is excellent for isotherms at all the temperatures and pressures studied in this work (Figures 5a, 5b, and 5c). The deviation for heats of adsorption is within 2 kJ/mol compared with the experimental data from Sun et al.57 and Dunne et al (Figure 5d).58 DDR is another commonly studied zeolite. The comparison of our simulated and experimental data for CH4 in DDR is shown in Figure 6. CCFF_AA and CCFF_UA perform well for predicting isotherms at 195 K (Figure 6a), but slightly underestimate loadings at higher temperatures if compared with the experimental results from Bergh et al.59 When we include the experimental data from Himenno et al. in comparison,60 the agreement is excellent for both isotherms and heats of adsorption (Figures 6b and 6c). Based on the validation results for CH4 in these zeolites, the FF transferability with respect to adsorption isotherms is good.

13 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 28

Figure 4. Comparison of simulated (CCFF_AA and CCFF_UA) and experimental adsorption isotherms and heats of adsorption of CH4 in siliceous CHA (a and b), and ITQ-29 (siliceous LTA) (c and d). The experimental data are from Hedin et al.,53 Maghsoudi et al.,54 Pham et al.,55 and Palomino et al.56 Lines are drawn to guide the eye.

14 ACS Paragon Plus Environment

Page 15 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 5. Comparison of simulated (CCFF_AA and CCFF_UA) and experimental adsorption isotherms and heats of adsorption of CH4 in siliceous MFI. The experimental data are from Zhu et al.,61 Sun et al.,57 and Dunne et al.58 Lines are drawn to guide the eye.

15 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 28

Figure 6. Comparison of simulated (CCFF_AA and CCFF_UA) and experimental adsorption isotherms and heats of adsorption of CH4 in siliceous DDR. The experimental data are from Himeno et al.60 and Bergh et al.59 Lines are drawn to guide the eye.

16 ACS Paragon Plus Environment

Page 17 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

We have also checked the influence of the first-principles method on the performance of FFs. The effect of QM method is significant. For instance, while the FFs derived from DFT/CC work well for CH4 in CHA and ITQ-29, the ones derived from DFT-D2 overestimate the isotherms and heats of adsorption (Figure S4 in Supporting Information). This is different from the case of CO2, where the FFs derived from both DFT/CC and DFT-D2 give similar results.10,11 The difference in performance of DFT/CC and DFT-D2 for CH4 and CO2 in zeolites observed from our studies is consistent with that from Bludský and co-workers.30,35,62 This is probably because CH4 in zeolites is dispersion dominated and far more sensitive to force field van der Waals parameters compared to CO2 in zeolites, where both dispersion and electrostatic interactions dominate.

3.2 Prediction of Diffusion Properties To examine the applicability of the FFs for predicting diffusion properties, we performed standard MD simulations for CH4 in Si-CHA and Si-LTA zeolites, calculated self-diffusion coefficients,63 and compared them with experimental data obtained by Hedin et al. using pulsed field gradient (PFG) NMR.53 These two zeolites include 8-ring windows, through which CH4 molecules diffuse from cage to cage. It has been shown that for 8-ring zeolites, the flexibility of framework would have influence on molecular diffusion, so we tested both rigid and flexible frameworks. The modified Hill-Sauer FF from our recent work was used to account for flexibility of the zeolites,64 since it improves the description of pore sizes and framework flexibility in small pore zeolites relative to the original Hill-Suer FF.65 Because the atomic charges of the framework in our CCFF_AA are not consistent with those from Hill and Sauer FF,65 we only investigate flexible framework when using the CCFF_UA. Our MD simulations were performed using LAMMPS,66 and further details are given in the Supporting Information. The results for ITQ-29 are reported in Table 5. The experimental diffusion coefficient for CH4 in ITQ-29 was measured at the loading of 0.4 molecules per cage.53 To get better statistics, in simulations we calculated the diffusion coefficients at infinite dilution and the loading of 1 molecule per cage, and we expect that the simulated results should provide the lower and upper bounds for the result with the actual fractional loading (see section S3.2 in Supporting Information). As is seen, both CCFF_AA and CCFF_UA with rigid framework do a good job in predicting the diffusion coefficients when compared with the experimental data. When flexibility is taken into account, slightly better agreement is achieved, although the influence of flexibility is not significant for this case. Table 6 summarizes the results for siliceous CHA. CCFF_AA and CCFF_UA with rigid framework underestimate the diffusion coefficients compared to experimental data, but the values are in the same order of magnitude. When flexibility is taken into account, CCFF_UA performs better compared with the result using rigid framework. The 8-ring window sizes of Si-CHA and Si-LTA are 3.70 × 4.17 Å 17 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 28

and 4.00 × 4.22 Å, respectively,53 and the kinetic diameter of CH4 is 3.8 Å.67 Our previous studies show that the framework flexibility plays a crucial role when the size of the diffusing molecule approaches the size of the window.42,63 The window size of Si-CHA is more comparable to the kinetic diameter of CH4 than that of Si-LTA, and therefore the influence of framework flexibility is more significant for CH4 in Si-CHA than in Si-LTA. Multiple previous FFs have been developed for adsorption of CH4 in zeolites, and their performance for predicting adsorption properties is generally satisfactory.13,14,47,68-83 Here we checked their performance for predicting diffusion properties by calculating self-diffusion coefficients for CH4 in ITQ-29 and CHA at 300 K. The details about the FFs we used are summarized in Table S10 in Supporting Information. Since some of these FFs predict slow diffusion (less than 10-12 m2/s) that is beyond the time scale of standard MD, dynamically-corrected transition state theory (dcTST) method was used for these simulations.84-86 In dcTST simulations the loading of CH4 was at infinite dilution, and flexibility of framework was taken into account as described above. More details about dcTST are given in the section S3.3 in Supporting Information. Figure 7 summarizes the results of simulated self-diffusion coefficients as a function of the FF parameter, σLJ (Å) of CH4-Oz, for the various FFs we tested. It can be seen that the predicted selfdiffusion coefficients using different FFs differ by 8 and 10 orders of magnitude for ITQ-29 and CHA, respectively. The numerical values are reported in Table S11 in Supporting Information. By comparing with experimental data from Hedin et al.,53 most of earlier FFs that were developed to reproduce experimental adsorption do not perform reliably for predicting diffusivities. Our CCFF_AA and CCFF_UA stand out in the comparison. This is a striking example of how fitting a FF to only adsorption data can give predictions that are highly inaccurate for predicting molecular diffusion.

18 ACS Paragon Plus Environment

Page 19 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Table 5: Self-diffusion coefficients (m2/s) for CH4 in ITQ-29 (siliceous LTA) at 301 K, calculated using CCFF_AA and CCFF_UA with rigid and flexible frameworks. Standard Molecular Dynamics Exp. dataa Loading Infinite dilution 1.0 0.4 (molecule/cage) -10 -10 CCFF_AA (1.24 ± 0.25) × 10 (2.09 ± 0.53) × 10 (rigid framework) -10 -10 -10 CCFF_UA 1.42 × 10 (1.45 ± 0.34) × 10 (1.48 ± 0.54) × 10 (rigid framework) -10 -10 CCFF_UA (1.44 ± 0.29) × 10 (1.49 ± 0.26) × 10 (flexible framework) a Hedin et al. PFG-NMR.53 Table 6: Self-diffusion coefficients (m2/s) for CH4 in siliceous CHA at 301 K, calculated using CCFF_AA and CCFF_UA with rigid and flexible frameworks. Standard Molecular Dynamics Exp. dataa Loading Infinite dilution 1.0 0.4 (molecule/cage) -11 -11 CCFF_AA (0.27 ± 0.04) × 10 (0.57 ± 0.13) × 10 (rigid framework) -11 -11 -11 CCFF_UA 1.07 × 10 (0.35 ± 0.06) × 10 (0.63 ± 0.15) × 10 (rigid framework) -11 -11 CCFF_UA (0.49 ± 0.07) × 10 (0.96 ± 0.12) × 10 (flexible framework) a Hedin et al. PFG-NMR.53

19 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 28

Figure 7: CH4 self-diffusion coefficients at infinite dilution and 300 K in flexible zeolite ITQ-29 (black circles) and CHA (red squares) computed using CCFF (UA: filled symbols, AA: hashed symbols) and compared to 15 other force fields from literature (see Tables S10 and S11 in Supporting Information for details). Hedin et al.53 experimental self-diffusion coefficients are shown as dashed lines.

4. Conclusions We have developed first-principles-derived force fields for CH4 in siliceous zeolites that can accurately predict both adsorption and diffusion properties. Two models are tested for CH4 and zeolites. One uses a united-atom (UA) model for CH4 and includes only vdW interactions between CH4 and framework oxygen atoms. The other uses a rigid all-atom (AA) model with atomic charges for each atom of CH4 molecule and includes vdW and Coulombic interactions with both silicon and oxygen atoms between CH4 and zeolite. The results show that the first model that has stronger approximation performs similarly as the second model in predicting both adsorption and diffusion properties. The developed force fields show good transferability among different topologies. The force fields derived from DFT/CC method work well for CH4 in siliceous zeolites, while the ones derived from DFT-D2 method overestimate the isotherms and heats of adsorption. This observation is different from that of CO2 in 20 ACS Paragon Plus Environment

Page 21 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

siliceous zeolites, where the force fields derived from DFT/CC and DFT-D2 give similar results. Finally, a comparison between our force fields and those in literature indicate that most of earlier force fields that were developed to reproduce experimental adsorption do not perform reliably for predicting diffusivities of CH4 in silica zeolites. The extension of the methodology to higher hydrocarbon molecules and cationic zeolites is undergoing.

21 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 28

ASSOCIATED CONTENT Supporting Information Detailed description of DFT/CC method for CH4 in zeolite, periodic DFT calculations, and classical simulations.

AUTHOR INFORMATION Corresponding Authors [email protected] Author Contributions §

These authors contribute equally to this work

ORCID iD Hanjun Fang: 0000-0002-7125-6404 Salah E. Boulfelfel: 0000-0003-1412-1826 Rohan Awati: 0000-0001-9361-0622 David S. Sholl: 0000-0002-2771-9168

Notes The authors declare no competing financial interests.

ACKNOWLEGMENTS HF, RA, SB, and DSS acknowledge funding from ExxonMobil Research and Engineering. We thank Preeti Kamakoti and Lucas Koziol for helpful discussions.

22 ACS Paragon Plus Environment

Page 23 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

REFERENCES (1) Allen, M. P.; Tildesley, D. J. Computer Simulation of Liquids; Clarendon Press: Oxford, U.K., 1987. (2) Frenkel, D.; Smit, B. Understanding Molecular Simulation: From Algorithms to Applications 2nd ed.; Academic Press: San Diego, CA, 2002. (3) Fang, H. J.; Demir, H.; Kamakoti, P.; Sholl, D. S. Recent Developments in First-Principles Force Fields for Molecules in Nanoporous Materials. J. Mater. Chem. A 2014, 2, 274-291. (4) Dzubak, A. L.; Lin, L.-C.; Kim, J.; Swisher, J. A.; Poloni, R.; Maximoff, S. N.; Smit, B.; Gagliardi, L. Ab Initio Carbon Capture in Open-Site Metal–Organic Frameworks. Nat. Chem. 2012, 4, 810-816. (5) Chen, L.; Morrison, C. A.; Düren, T. Improving Predictions of Gas Adsorption in Metal–Organic Frameworks with Coordinatively Unsaturated Metal Sites: Model Potentials, Ab Initio Parameterization, and Gcmc Simulations. J. Phys. Chem. C 2012, 116, 18899-18909. (6) Lin, L. C.; Lee, K.; Gagliardi, L.; Neaton, J. B.; Smit, B. Force-Field Development from Electronic Structure Calculations with Periodic Boundary Conditions: Applications to Gaseous Adsorption and Transport in Metal-Organic Frameworks. J. Chem. Theory Comput. 2014, 10, 1477-1488. (7) Schmidt, J. R.; Yu, K.; McDaniel, J. G. Transferable Next-Generation Force Fields from Simple Liquids to Complex Materials. Accounts. Chem. Res. 2015, 48, 548-556. (8) Odoh, S. O.; Cramer, C. J.; Truhlar, D. G.; Gagliardi, L. Quantum-Chemical Characterization of the Properties and Reactivities of Metal-Organic Frameworks. Chem. Rev. 2015, 115, 6051-6111. (9) Sun, L.; Deng, W. Q. Recent Developments of First-Principles Force Fields. WIREs Comput. Mol. Sci. 2017, 7. (10) Fang, H.; Kamakoti, P.; Zang, J.; Cundy, S.; Paur, C.; Ravikovitch, P. I.; Sholl, D. S. Prediction of Co2 Adsorption Properties in Zeolites Using Force Fields Derived from Periodic Dispersion-Corrected Dft Calculations. J. Phys. Chem. C 2012, 116, 10692-10701. (11) Fang, H. J.; Kamakoti, P.; Ravikovitch, P. I.; Aronson, M.; Paur, C.; Sholl, D. S. First Principles Derived, Transferable Force Fields for Co2 Adsorption in Na-Exchanged Cationic Zeolites. Phys. Chem. Chem. Phys. 2013, 15, 12882-12894. (12) Fang, H. J.; Kulkarni, A.; Kamakoti, P.; Awati, R.; Ravikovitch, P. I.; Sholl, D. S. Identification of HighCo2-Capacity Cationic Zeolites by Accurate Computational Screening. Chem. Mater. 2016, 28, 3887-3896. (13) Dubbeldam, D.; Calero, S.; Vlugt, T. J. H.; Krishna, R.; Maesen, T. L. M.; Beerdsen, E.; Smit, B. Force Field Parametrization through Fitting on Inflection Points in Isotherms. Phys. Rev. Lett. 2004, 93, 088302. (14) Dubbeldam, D.; Calero, S.; Vlugt, T. J. H.; Krishna, R.; Maesen, T. L. M.; Smit, B. United Atom Force Field for Alkanes in Nanoporous Materials. J. Phys. Chem. B 2004, 108, 12301-12313. (15) Jee, S. E.; Sholl, D. S. Carbon Dioxide and Methane Transport in Ddr Zeolite: Insights from Molecular Simulations into Carbon Dioxide Separations in Small Pore Zeolites. J. Am. Chem. Soc. 2009, 131, 78967904. (16) Diaz-Cabanas, M. J.; Barrett, P. A.; Camblor, M. A. Synthesis and Structure of Pure Sio2 Chabazite: The Sio2 Polymorph with the Lowest Framework Density. Chem. Commun. 1998, 1881-1882. (17) Goltl, F.; Hafner, J. Alkane Adsorption in Na-Exchanged Chabazite: The Influence of Dispersion Forces. J. Chem. Phys. 2011, 134, 064102. (18) Goltl, F.; Gruneis, A.; Bucko, T.; Hafner, J. Van Der Waals Interactions between Hydrocarbon Molecules and Zeolites: Periodic Calculations at Different Levels of Theory, from Density Functional Theory to the Random Phase Approximation and Moller-Plesset Perturbation Theory. J. Chem. Phys. 2012, 137. (19) Goltl, F.; Hafner, J. Modelling the Adsorption of Short Alkanes in Protonated Chabazite: The Impact of Dispersion Forces and Temperature. Micropor. Mesopor. Mat. 2013, 166, 176-184. (20) Goltl, F.; Sautet, P. Modeling the Adsorption of Short Alkanes in the Zeolite Ssz-13 Using "Van Der Waals" Dft Exchange Correlation Functionals: Understanding the Advantages and Limitations of Such Functionals. J. Chem. Phys. 2014, 140. 23 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 28

(21) Shang, J.; Li, G.; Singh, R.; Xiao, P.; Danaci, D.; Liu, J. Z.; Webley, P. A. Adsorption of Co2, N-2, and Ch4 in Cs-Exchanged Chabazite: A Combination of Van Der Waals Density Functional Theory Calculations and Experiment Study. J. Chem. Phys. 2014, 140. (22) Piccini, G.; Alessio, M.; Sauer, J.; Zhi, Y. C.; Liu, Y.; Kolvenbach, R.; Jentys, A.; Lercher, J. A. Accurate Adsorption Thermodynamics of Small Alkanes in Zeolites. Ab Initio Theory and Experiment for HChabazite. J. Phys. Chem. C. 2015, 119, 6128-6137. (23) Fischer, M.; Bell, R. G. A Dft-D Study of the Interaction of Methane, Carbon Monoxide, and Nitrogen with Cation-Exchanged Sapo-34. Z. Kristallog. – Cryst. Mater. 2015, 230, 311-323. (24) Chiu, C. C.; Vayssilov, G. N.; Genest, A.; Borgna, A.; Rosch, N. Predicting Adsorption Enthalpies on Silicalite and Hzsm-5: A Benchmark Study on Dft Strategies Addressing Dispersion Interactions. J. Comput. Chem. 2014, 35, 809-819. (25) Zhang, Y. W.; Yu, J. Y.; Yeh, Y. H.; Gorte, R. J.; Rangarajan, S.; Mavrikakis, M. An Adsorption Study of Ch4 on Zsm-5, Mor, and Zsm-12 Zeolites. J. Phys. Chem. C. 2015, 119, 28970-28978. (26) Bludsky, O.; Rubes, M.; Soldan, P.; Nachtigall, P. Investigation of the Benzene-Dimer Potential Energy Surface: Dft/Ccsd(T) Correction Scheme. J. Chem. Phys. 2008, 128, 114102. (27) Rubes, M.; Nachtigall, P.; Vondrasek, J.; Bludsky, O. Structure and Stability of the Water-Graphite Complexes. J. Phys. Chem. C. 2009, 113, 8412-8419. (28) Rubes, M.; Kysilka, J.; Nachtigall, P.; Bludsky, O. Dft/Cc Investigation of Physical Adsorption on a Graphite (0001) Surface. Phys. Chem. Chem. Phys. 2010, 12, 6438-6444. (29) Grajciar, L.; Bludsky, O.; Nachtigall, P. Water Adsorption on Coordinatively Unsaturated Sites in Cubtc Mof. J. Phys. Chem. Lett. 2010, 1, 3354-3359. (30) Pulido, A.; Delgado, M. R.; Bludsky, O.; Rubes, M.; Nachtigall, P.; Arean, C. O. Combined Dft/Cc and Ir Spectroscopic Studies on Carbon Dioxide Adsorption on the Zeolite H-Fer. Energy Environ. Sci. 2009, 2, 1187-1195. (31) Cejka, J.; Zukal, A.; Pulido, A.; Gil, B.; Nachtigall, P.; Bludsky, O.; Rubes, M. Experimental and Theoretical Determination of Adsorption Heats of Co(2) over Alkali Metal Exchanged Ferrierites with Different Si/Al Ratio. Phys. Chem. Chem. Phys. 2010, 12, 6413-6422. (32) Arean, C. O.; Delgado, M. R.; Bibiloni, G. F.; Bludsky, O.; Nachtigall, P. Variable-Temperature Ir Spectroscopic and Theoretical Studies on Co(2) Adsorbed in Zeolite K-Fer. ChemPhysChem 2011, 12, 1435-1443. (33) Zukal, A.; Arean, C. O.; Delgado, M. R.; Nachtigall, P.; Pulido, A.; Mayerova, J.; Cejka, J. Combined Volumetric, Infrared Spectroscopic and Theoretical Investigation of Co(2) Adsorption on Na-a Zeolite. Micropor. Mesopor. Mat. 2011, 146, 97-105. (34) Chen, L. J.; Grajciar, L.; Nachtigall, P.; Duren, T. Accurate Prediction of Methane Adsorption in a Metal-Organic Framework with Unsaturated Metal Sites by Direct Implementation of an Ab Initio Derived Potential Energy Surface in Gcmc Simulation. J. Phys. Chem. C. 2011, 115, 23074-23080. (35) Rubes, M.; Trachta, M.; Koudelkova, E.; Bulanek, R.; Kasneryk, V.; Bludsky, O. Methane Adsorption in Ador Zeolites: A Combined Experimental and Dft/Cc Study. Phys. Chem. Chem. Phys. 2017, 19, 1653316540. (36) Tkatchenko, A.; Scheffler, M. Accurate Molecular Van Der Waals Interactions from Ground-State Electron Density and Free-Atom Reference Data. Phys. Rev. Lett. 2009, 102. (37) Bucko, T.; Lebegue, S.; Angyan, J. G.; Hafner, J. Extending the Applicability of the TkatchenkoScheffler Dispersion Correction Via Iterative Hirshfeld Partitioning. J. Chem. Phys. 2014, 141. (38) Tkatchenko, A.; DiStasio, R. A.; Car, R.; Scheffler, M. Accurate and Efficient Method for Many-Body Van Der Waals Interactions. Phys. Rev. Lett. 2012, 108. (39) Ambrosetti, A.; Reilly, A. M.; DiStasio, R. A.; Tkatchenko, A. Long-Range Correlation Energy Calculated from Coupled Atomic Response Functions. J. Chem. Phys. 2014, 140.

24 ACS Paragon Plus Environment

Page 25 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(40) Steinmann, S. N.; Corminboeuf, C. A Generalized-Gradient Approximation Exchange Hole Model for Dispersion Coefficients. J. Chem. Phys. 2011, 134. (41) Steinmann, S. N.; Corminboeuf, C. Comprehensive Bench Marking of a Density-Dependent Dispersion Correction. J. Chem. Theory Comput. 2011, 7, 3567-3577. (42) Awati, R. V.; Ravikovitch, P. I.; Sholl, D. S. Efficient and Accurate Methods for Characterizing Effects of Framework Flexibility on Molecular Diffusion in Zeolites: Ch4 Diffusion in Eight Member Ring Zeolites. J. Phys. Chem. C. 2013, 117, 13462-13473. (43) Kresse, G.; Hafner, J. Ab-Initio Molecular-Dynamics Simulation of the Liquid-Metal AmorphousSemiconductor Transition in Germanium. Phys. Rev. B 1994, 49, 14251-14269. (44) Kresse, G.; Furthmuller, J. Efficient Iterative Schemes for Ab Initio Total-Energy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B 1996, 54, 11169-11186. (45) Bucko, T.; Benco, L.; Hafner, J.; Angyan, J. G. Monomolecular Cracking of Propane over Acidic Chabazite: An Ab Initio Molecular Dynamics and Transition Path Sampling Study. J. Catal. 2011, 279, 220228. (46) Martin, M. G.; Siepmann, J. I. Transferable Potentials for Phase Equilibria. 1. United-Atom Description of N-Alkanes. J. Phys. Chem. B 1998, 102, 2569-2577. (47) Vlugt, T. J. H.; Krishna, R.; Smit, B. Molecular Simulations of Adsorption Isotherms for Linear and Branched Alkanes and Their Mixtures in Silicalite. J. Phys. Chem. B 1999, 103, 1102-1118. (48) Grimme, S. Semiempirical Gga-Type Density Functional Constructed with a Long-Range Dispersion Correction. J. Comput. Chem. 2006, 27, 1787-1799. (49) Manz, T. A.; Sholl, D. S. Chemically Meaningful Atomic Charges That Reproduce the Electrostatic Potential in Periodic and Nonperiodic Materials. J. Chem. Theory Comput. 2010, 6, 2455-2468. (50) Manz, T. A.; Sholl, D. S. Methods for Computing Accurate Atomic Spin Moments for Collinear and Noncollinear Magnetism in Periodic and Nonperiodic Materials. J. Chem. Theory Comput. 2011, 7, 41464164. (51) Manz, T. A.; Sholl, D. S. Improved Atoms-in-Molecule Charge Partitioning Functional for Simultaneously Reproducing the Electrostatic Potential and Chemical States in Periodic and Nonperiodic Materials. J. Chem. Theory Comput. 2012, 8, 2844-2867. (52) Cygan, R. T.; Liang, J.-J.; Kalinichev, A. G. Molecular Models of Hydroxide, Oxyhydroxide, and Clay Phases and the Development of a General Force Field. J. Phys. Chem. B 2004, 108, 1255-1266. (53) Hedin, N.; DeMartin, G. J.; Roth, W. J.; Strohmaier, K. G.; Reyes, S. C. Pfg Nmr Self-Diffusion of Small Hydrocarbons in High Silica Ddr, Cha and Lta Structures. Micropor. Mesopor. Mat. 2008, 109, 327-334. (54) Maghsoudi, H.; Soltanieh, M.; Bozorgzadeh, H.; Mohamadalizadeh, A. Adsorption Isotherms and Ideal Selectivities of Hydrogen Sulfide and Carbon Dioxide over Methane for the Si-Cha Zeolite: Comparison of Carbon Dioxide and Methane Adsorption with the All-Silica Dd3r Zeolite. Adsorption 2013, 19, 1045-1053. (55) Pham, T. D.; Lobo, R. F. Adsorption Equilibria of Co2 and Small Hydrocarbons in Aei-, Cha-, Sit-, and Rro-Type Siliceous Zeolites. Micropor. Mesopor. Mat. 2016, 236, 100-108. (56) Palomino, M.; Corma, A.; Rey, F.; Valencia, S. New Insights on Co2-Methane Separation Using Lta Zeolites with Different Si/Al Ratios and a First Comparison with Mofs. Langmuir 2010, 26, 1910-1917. (57) Sun, M. S.; Shah, D. B.; Xu, H. H.; Talu, O. Adsorption Equilibria of C(1) to C(4) Alkanes, Co(2), and Sf(6) on Silicalite. J. Phys. Chem. B 1998, 102, 1466-1473. (58) Dunne, J. A.; Mariwals, R.; Rao, M.; Sircar, S.; Gorte, R. J.; Myers, A. L. Calorimetric Heats of Adsorption and Adsorption Isotherms .1. O-2, N-2, Ar, Co2, Ch4, C2h6 and Sf6 on Silicalite. Langmuir 1996, 12, 5888-5895. (59) van den Bergh, J.; Zhu, W.; Gascon, J.; Moulijn, J. A.; Kapteijn, F. Separation and Permeation Characteristics of a Dd3r Zeolite Membrane. J. Membrane Sci. 2008, 316, 35-45.

25 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 28

(60) Himeno, S.; Tomita, T.; Suzuki, K.; Yoshida, S. Characterization and Selectivity for Methane and Carbon Dioxide Adsorption on the All-Silica Dd3r Zeolite. Micropor. Mesopor. Mat. 2007, 98, 62-69. (61) Zhu, W. D.; Hrabanek, P.; Gora, L.; Kapteijn, F.; Moulijn, J. A. Role of Adsorption in the Permeation of Ch4 and Co2 through a Silicalite-1 Membrane. Ind. Eng. Chem. Res. 2006, 45, 767-776. (62) Hermann, J.; Bludsky, O. A Novel Correction Scheme for Dft: A Combined Vdw-Df/Ccsd(T) Approach. J. Chem. Phys. 2013, 139. (63) Boulfelfel, S. E.; Ravikovitch, P. I.; Sholl, D. S. Modeling Diffusion of Linear Hydrocarbons in Silica Zeolite Lta Using Transition Path Sampling. J. Phys. Chem. C. 2015, 119, 15643-15653. (64) Boulfelfel, S. E.; Ravikovitch, P. I.; Koziol, L.; Sholl, D. S. Improved Hill-Sauer Force Field for Accurate Description of Pores in 8-Ring Zeolites. J. Phys. Chem. C. 2016, 120, 14140-14148. (65) Hill, J. R.; Sauer, J. Molecular Mechanics Potential for Silica and Zeolite Catalysts Based on Ab-Initio Calculations .2. Aluminosilicates. J. Phys. Chem. 1995, 99, 9536-9550. (66) Plimpton, S. Fast Parallel Algorithms for Short-Range Molecular-Dynamics. J. Comput. Phys. 1995, 117, 1-19. (67) Ismail, A. F.; Khulbe, K. C.; Matsuura, T. Gas Separation Membranes : Polymeric and Inorganic; Springer: Cham, 2015. (68) Yashonath, S.; Bandyopadhyay, S. Surprising Diffusion Behavior in the Restricted Regions of Silicalite. Chem. Phys. Lett. 1994, 228, 284-288. (69) Fritzsche, S.; Haberlandt, R.; Karger, J.; Pfeifer, H.; Heinzinger, K. On the Diffusion Mechanism of Methane in a Cation-Free Zeolite of Type Zk4. Chem. Phys. 1993, 174, 229-236. (70) Haberlandt, R. Transport Processes in Porous Media: Diffusion in Zeolites. Thin Solid Films 1998, 330, 34-45. (71) Ruthven, D. M.; Derrah, R. I. Transition-State Theory of Zeolitic Diffusion Diffusion of Ch4 and Cff Ch4 and Cf4 in 5a Zeolite. J. Chem. Soc., Faraday Trans. 1 1972, 68, 2332-2343. (72) Fritzsche, S.; Haberlandt, R.; Wolfsberg, M. Equilibration of the Kinetic Energy in Small Zeolite Cavities - the Thermalization Effect of Lattice Vibrations and of Mutual Interaction in the Diffusion of Methane in a Cation-Free Lta Zeolite. Chem. Phys. 2000, 253, 283-294. (73) Jost, S.; Fritzsche, S.; Haberlandt, R. An Md Study on the Diffusion of a Mixture of Methane and Xenon in Silicalite. Chem. Phys. Lett. 1997, 279, 385-388. (74) Hussain, I.; Titiloye, J. O. Molecular Dynamics Simulations of the Adsorption and Diffusion Behavior of Pure and Mixed Alkanes in Silicalite. Micropor. Mesopor. Mat. 2005, 85, 143-156. (75) Pascual, P.; Ungerer, P.; Tavitian, B.; Boutin, A. Development of a Transferable Guest-Host Force Field for Adsorption of Hydrocarbons in Zeolites. Ii. Prediction of Alkenes Adsorption and Alkane/Alkene Selectivity in Silicalite. J. Phys. Chem. B 2004, 108, 393-398. (76) Leroy, F.; Rousseau, B.; Fuchs, A. H. Self-Diffusion of N-Alkanes in Silicalite Using Molecular Dynamics Simulation: A Comparison between Rigid and Flexible Frameworks. Phys. Chem. Chem. Phys. 2004, 6, 775-783. (77) Goodbody, S. J.; Watanabe, K.; Macgowan, D.; Walton, J. P. R. B.; Quirke, N. Molecular Simulation of Methane and Butane in Silicalite. J. Chem. Soc., Faraday Trans. 1991, 87, 1951-1958. (78) Fritzsche, S.; Karger, J. Tracing Memory Effects in Correlated Diffusion Anisotropy in Mfi-Type Zeolites by Md Simulation. J. Phys. Chem. B 2003, 107, 3515-3521. (79) Kamat, M.; Dang, W. J.; Keffer, D. Agreement between Analytical Theory and Molecular Dynamics Simulation for Adsorption and Diffusion in Crystalline Nanoporous Materials. J. Phys. Chem. B 2004, 108, 376-386. (80) Bai, P.; Tsapatsis, M.; Siepmann, J. I. Trappe-Zeo: Transferable Potentials for Phase Equilibria Force Field for All-Silica Zeolites. J. Phys. Chem. C. 2013, 117, 24375-24387.

26 ACS Paragon Plus Environment

Page 27 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(81) Martin, M. G.; Thompson, A. P.; Nenoff, T. M. Effect of Pressure, Membrane Thickness, and Placement of Control Volumes on the Flux of Methane through Thin Silicalite Membranes: A Dual Control Volume Grand Canonical Molecular Dynamics Study. J. Chem. Phys. 2001, 114, 7174-7181. (82) Demontis, P.; Suffritti, G. B.; Fois, E. S.; Quartieri, S. Molecular-Dynamics Studies on Zeolites .6. Temperature-Dependence of Diffusion of Methane in Silicalite. J. Phys. Chem. 1992, 96, 1482-1490. (83) Talu, O.; Myers, A. L. Reference Potentials for Adsorption of Helium, Argon, Methane, and Krypton in High-Silica Zeolites. Colloid Surface A 2001, 187, 83-93. (84) Dubbeldam, D.; Beerdsen, E.; Vlugt, T. J. H.; Smit, B. Molecular Simulation of Loading-Dependent Diffusion in Nanoporous Materials Using Extended Dynamically Corrected Transition State Theory. J. Chem. Phys. 2005, 122. (85) Camp, J. S.; Sholl, D. S. Transition State Theory Methods to Measure Diffusion in Flexible Nanoporous Materials: Application to a Porous Organic Cage Crystal. J. Phys. Chem. C. 2016, 120, 11101120. (86) Verploegh, R. J.; Wu, Y.; Sholl, D. S. Lattice-Gas Modeling of Adsorbate Diffusion in Mixed-Linker Zeolitic Imidazolate Frameworks: Effect of Local Imidazolate Ordering. Langmuir 2017, 33, 6481-6491.

27 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 28

TOC Graphic

28 ACS Paragon Plus Environment