Robust Computation and Analysis of Vibrational Spectra of Layered Framework Materials including Host-Guest Interactions

Layered framework materials, a rapidly advancing class of porous materials, are composed of molecular components stitched together via covalent bonds and are usually synthesized through wet-chemical methods. Computational infrared (IR) and Raman spectra are among the most important characterization tools for this materials class. Besides the a priori known spectra of the molecular building blocks and the solvent, they allow for in situ monitoring of the framework formation during synthesis. Therefore, they need to capture the additional peaks from host-guest interactions and the bands from emerging bonds between the molecular building blocks, verifying the successful synthesis of the desired material. In this work, we propose a robust computational framework based on ab initio molecular dynamics (AIMD), where we compute IR and Raman spectra from the time-correlation functions of dipole moments and polarizability tensors, respectively. As a case study, we apply our methodology to a covalent organic framework (COF) material, COF-1, and present its AIMD-computed IR and Raman spectra with and without 1,4-dioxane solvent molecules in its pores. To determine robust settings, we meticulously validate our model and explore how stacking disorder and different methods for computing dipole moments and polarizabilities affect IR and Raman intensities. Using our robust computational protocol, we achieve excellent agreement with experimental data. Furthermore, we illustrate how the computed spectra can be dissected into individual contributions from the solvent molecules, the molecular building blocks of COF-1, and the bonds connecting them.