Benchmark Data for Interactions in Zeolite Model Complexes and Their Use for Assessment and Validation of Electronic Structure Methods
The research described in this product was performed in part in the Environmental Molecular Sciences Laboratory, a national scientific user facility sponsored by the Department of Energy's Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory. We present benchmark binding energies for five zeolite model complexes, with four of the adsorbates bound noncovalently and one covalently. The binding energies were determined as the sum of the infinite-basis-set limit of Møller-Plesset second-order perturbation theory (MP2) energies and a CCSD(T) correction term evaluated with the aug-cc-pVDZ basis set. The basis set limit of MP2 energies was determined by two-point extrapolation using the aug-cc-pVXZ basis sets for X ) D and T and separate extrapolation of the Hartree- Fock and correlation energies. We found that correlation contributions beyond MP2 to the final binding energies are small; their magnitude is in the range of 0.02-1.0 kcal/mol. For the MP2 method to describe the interactions in these zeolite model systems accurately, one needs to use a basis set at least the size of aug-cc-pVTZ in conjunction with counterpoise corrections. Final binding energies (kcal/mol) of the model complexes are in the range of 3.5-19.5 kcal/mol, and they were used as reference data to test 6 wave function methods and 41 density functionals. Among the tested density functional methods, M06-L/6-31+G(d,p) gives a mean unsigned error (MUE) without counterpoise correction of 0.87 kcal/mol. With counterpoise corrections, the M06-2X and M05-2X functionals give the best performance. The MUE with counterpoise corrections for the M06-2X/6-311+G(2df,2p)//MP2/6-311+G(2df,2p) level of theory is 0.39 kcal/mol. With the DFT/6-31+G- (d,p) geometries and the 6-311+G(2df,2p) basis set, M05-2X and M06-2X give MUEs with counterpoise corrections of 0.40 and 0.52 kcal/mol, respectively. Tests against the binding energies of four complexes (two noncovalent and two covalent) of the adsoption of isobutene on a large 16T zeolite model cluster confirmed that M06-L, M06, M05-2X, and M06-2X are very promising quantum mechanical methods for hybrid quantum mechanical/molecular mechanical (QM/MM) simulations of zeolites. In fact the performance of these four Minnesota functionals, as compared to other high-quality functionals, is relatively even better for the larger 16T clusters than for the smaller 3T ones.
- Research Organization:
- Pacific Northwest National Lab. (PNNL), Richland, WA (United States). Environmental Molecular Sciences Lab. (EMSL)
- Sponsoring Organization:
- USDOE
- DOE Contract Number:
- AC05-76RL01830
- OSTI ID:
- 959203
- Journal Information:
- Journal of Physical Chemistry C, 112(17):6860-6868, Vol. 112, Issue 17; ISSN 1932-7447
- Country of Publication:
- United States
- Language:
- English
Similar Records
On the Structure and Geometry of Biomolecular Binding Motifs (Hydrogen-Bonding, Stacking, X-H· · ·π): WFT and DFT Calculations
Density functional theory study of the interaction of vinyl radical, ethyne, and ethene with benzene, aimed to define an affordable computational level to investigate stability trends in large van der Waals complexes