First principles, microkinetic, and experimental analysis of Lewis acid site speciation during ethanol dehydration on Sn-Beta zeolites

Density functional theory calculations are combined with kinetic measurements of ethanol dehydration to diethyl ether to identify the relative catalytic contributions of structurally distinct speciations of Sn sites in zeolite Beta frameworks. The structural complexities of the Beta framework require nonstandard techniques for entropy and energy calculations, including consideration of anharmonic effects in vibrational modes, employment of quasi-harmonic densities of states methods to evaluate entropies, and use of hybrid density functionals to evaluate binding energies. Calculated energies and entropies are used to construct a microkinetic model that is iteratively refined to identify all kinetically and thermodynamically sensitive reaction steps and intermediates which are subsequently treated with the higher-level methods. The rate and equilibrium constants obtained from this tiered approach agree well with measured reaction orders in ethanol and water. Site balances provide evidence for the interconversion of Sn sites between “closed” configurations that are tetra-coordinated to the framework, open configurations formed by hydrolysis that are tri-coordinated to the framework and contain a hydroxyl ligand and proximal silanol group (“hydroxy-open”), and open tri-coordinate configurations formed by reaction with ethanol, yielding an ethoxy ligand and proximal silanol group (“ethoxy-open”). Closed, hydroxy-open, and ethoxy-open Sn sites adsorb ethanol via distinct modes and react via distinct pathways, with the consequence that prevalent dehydration pathways depend on the speciation of Sn sites under reaction conditions. The kinetic modeling indicates that, under the conditions studied (404 K, 0.5–35 kPa ethanol, 0.1–50 kPa water), bimolecular dehydration on the closed Sn site is the sole kinetically-relevant step, and ethanol, ethanol-ethanol dimers, and ethanol-water dimers are the most abundant surface intermediates. These results highlight the importance of considering the distribution of heteroatom coordination modes to zeolite frameworks under reaction conditions, as well as pathways for their interconversion in the presence of reacting molecules, to obtain more robust mechanistic and kinetic interpretations of catalytic pathways in Lewis acid zeolites.