Accurate structures and energetics of neutral-framework zeotypes from dispersion-corrected DFT calculations.

Abstract

Density-functional theory (DFT) calculations incorporating a pairwise dispersion correction were employed to optimize the structures of various neutral-framework compounds with zeolite topologies. The calculations used the PBE functional for solids (PBEsol) in combination with two different dispersion correction schemes, the D2 correction devised by Grimme and the TS correction of Tkatchenko and Scheffler. In the first part of the study, a benchmarking of the DFT-optimized structures against experimental crystal structure data was carried out, considering a total of 14 structures (8 all-silica zeolites, 4 aluminophosphate zeotypes, and 2 dense phases). Both PBEsol-D2 and PBEsol-TS showed an excellent performance, improving significantly over the best-performing approach identified in a previous study (PBE-TS). The temperature dependence of lattice parameters and bond lengths was assessed for those zeotypes where the available experimental data permitted such an analysis. In most instances, the agreement between DFT and experiment improved when the experimental data were corrected for the effects of thermal motion and when low-temperature structure data rather than room-temperature structure data were used as a reference. In the second part, a benchmarking against experimental enthalpies of transition (with respect to α-quartz) was carried out for 16 all-silica zeolites. Excellent agreement was obtained with the PBEsol-D2 functional, with the overall error being in the same range as the experimental uncertainty. Altogether, PBEsol-D2 can be recommended as a computationally efficient DFT approach that simultaneously delivers accurate structures and energetics of neutral-framework zeotypes.