Publications

In this work we investigate anharmonic vibrational polaritons formed due to strong light-matter interactions in an optical cavity between radiation modes and anharmonic vibrations beyond the long-wavelength limit. We introduce a conceptually simple description of light-matter interactions, where spatially localized cavity radiation modes couple to localized vibrations. Within this theoretical framework, we employ self-consistent phonon theory and vibrational dynamical mean-field theory to efficiently simulate momentum-resolved vibrational-polariton spectra, including effects of anharmonicity. Numerical simulations in model systems demonstrate the accuracy and applicability of our approach.

The application of an external, oriented electric field has emerged as an attractive technique for manipulating chemical reactions. Because most applications occur in solution, a theory of electric field catalysis requires treatment of the solvent, whose interaction with both the external field and the reacting species modifies the reaction energetics and thus the reaction rate. Here, we formulate such a transition state theory using a dielectric continuum model, and we incorporate dynamical effects due to solvent motion via Grote–Hynes corrections. We apply our theory to the Menshutkin reaction between CH3I and pyridine, which is catalyzed by polar solvents, and to the symmetric SN2 reaction of F– with CH3F, which is inhibited by polar solvents. At low applied field strengths when the solvent responds linearly, our theory predicts near-complete quenching of electric field catalysis. However, a qualitative treatment of the nonlinear response (i.e., dielectric saturation) shows that catalysis can be recovered at appreciable field strengths as solvent molecules begin to align with the applied field direction. The dynamical correction to the rate constant is seen to vary nonmonotonically with increasing solvent polarity due to contrasting effects of the screening ability and the longitudinal relaxation time of the solvent.

The adsorption of CO on the surface of MgO has long been a model problem in surface chemistry. Here, we report periodic Gaussian-based calculations for this problem using second-order perturbation theory (MP2) and coupled-cluster theory with single and double excitations (CCSD) and perturbative triple excitations [CCSD(T)], with the latter two performed using a recently developed extension of the local natural orbital approximation to problems with periodic boundary conditions. The low cost of periodic local correlation calculations allows us to calculate the full CCSD(T) binding curve of CO approaching the surface of MgO (and thus the adsorption energy) and the two-dimensional potential energy surface (PES) as a function of the distance from the surface and the CO stretching coordinate. From the PES, we obtain the fundamental vibrational frequency of CO on MgO, whose shift from the gas phase value is a common experimental probe of surface adsorption. We find that CCSD(T) correctly predicts a positive frequency shift upon adsorption of $+14.7~\textrm{cm}^{-1}$, in excellent agreement with the experimental shift of $+14.3~\textrm{cm}^{-1}$. We use our CCSD(T) results to assess the accuracy of MP2, CCSD, and several density functional theory (DFT) approximations, including exchange correlation functionals and dispersion corrections. We find that MP2 and CCSD yield reasonable binding energies and frequency shifts, whereas many DFT calculations overestimate the magnitude of the adsorption energy by $5$ -- $15$~kJ/mol and predict a negative frequency shift of about $-20~\textrm{cm}^{-1}$, which we attribute to self-interaction-induced delocalization errors that are mildly ameliorated with hybrid functionals. Our findings highlight the accuracy and computational efficiency of the periodic local correlation for the simulation of surface chemistry with accurate wavefunction methods.