Hoffman Alexander E J, Temmerman Wim, Campbell Emma, Damin Alessandro Ali, Lezcano-Gonzalez Ines, Beale Andrew M, Bordiga Silvia, Hofkens Johan, Van Speybroeck Veronique
Center for Molecular Modeling, Ghent University, 9000 Ghent, Belgium.
Cardiff Catalysis Institute, Cardiff University, Cardiff CF10 3AT, United Kingdom.
J Chem Theory Comput. 2024 Jan 23;20(2):513-531. doi: 10.1021/acs.jctc.3c00942. Epub 2023 Dec 29.
Vibrational spectroscopy is an omnipresent spectroscopic technique to characterize functional nanostructured materials such as zeolites, metal-organic frameworks (MOFs), and metal-halide perovskites (MHPs). The resulting experimental spectra are usually complex, with both low-frequency framework modes and high-frequency functional group vibrations. Therefore, theoretically calculated spectra are often an essential element to elucidate the vibrational fingerprint. In principle, there are two possible approaches to calculate vibrational spectra: (i) a static approach that approximates the potential energy surface (PES) as a set of independent harmonic oscillators and (ii) a dynamic approach that explicitly samples the PES around equilibrium by integrating Newton's equations of motions. The dynamic approach considers anharmonic and temperature effects and provides a more genuine representation of materials at true operating conditions; however, such simulations come at a substantially increased computational cost. This is certainly true when forces and energy evaluations are performed at the quantum mechanical level. Molecular dynamics (MD) techniques have become more established within the field of computational chemistry. Yet, for the prediction of infrared (IR) and Raman spectra of nanostructured materials, their usage has been less explored and remain restricted to some isolated successes. Therefore, it is currently not clear which methodology should be used to accurately predict vibrational spectra for a given system. A comprehensive comparative study between various theoretical methods and experimental spectra for a broad set of nanostructured materials is so far lacking. To fill this gap, we herein present a concise overview on which methodology is suited to accurately predict vibrational spectra for a broad range of nanostructured materials and formulate a series of theoretical guidelines to this purpose. To this end, four different case studies are considered, each treating a particular material aspect, namely breathing in flexible MOFs, characterization of defects in the rigid MOF UiO-66, anharmonic vibrations in the metal-halide perovskite CsPbBr, and guest adsorption on the pores of the zeolite H-SSZ-13. For all four materials, in their guest- and defect-free state and at sufficiently low temperatures, both the static and dynamic approach yield qualitatively similar spectra in agreement with experimental results. When the temperature is increased, the harmonic approximation starts to fail for CsPbBr due to the presence of anharmonic phonon modes. Also, the spectroscopic fingerprints of defects and guest species are insufficiently well predicted by a simple harmonic model. Both phenomena flatten the potential energy surface (PES), which facilitates the transitions between metastable states, necessitating dynamic sampling. On the basis of the four case studies treated in this Review, we can propose the following theoretical guidelines to simulate accurate vibrational spectra of functional solid-state materials: (i) For nanostructured crystalline framework materials at low temperature, insights into the lattice dynamics can be obtained using a static approach relying on a few points on the PES and an independent set of harmonic oscillators. (ii) When the material is evaluated at higher temperatures or when additional complexity enters the system, , strong anharmonicity, defects, or guest species, the harmonic regime breaks down and dynamic sampling is required for a correct prediction of the phonon spectrum. These guidelines and their illustrations for prototype material classes can help experimental and theoretical researchers to enhance the knowledge obtained from a lattice dynamics study.
振动光谱是一种广泛应用的光谱技术,用于表征功能纳米结构材料,如沸石、金属有机框架(MOF)和金属卤化物钙钛矿(MHP)。所得的实验光谱通常很复杂,包含低频骨架模式和高频官能团振动。因此,理论计算光谱往往是阐明振动指纹的关键要素。原则上,有两种计算振动光谱的方法:(i)静态方法,将势能面(PES)近似为一组独立的简谐振子;(ii)动态方法,通过积分牛顿运动方程在平衡态附近显式采样PES。动态方法考虑了非谐性和温度效应,能更真实地反映材料在实际工作条件下的情况;然而,这种模拟的计算成本大幅增加。在量子力学层面进行力和能量评估时尤其如此。分子动力学(MD)技术在计算化学领域已更为成熟。然而,对于纳米结构材料的红外(IR)和拉曼光谱预测,其应用尚未得到充分探索,仍局限于一些孤立的成功案例。因此,目前尚不清楚对于给定系统应使用哪种方法来准确预测振动光谱。迄今为止,缺乏针对广泛的纳米结构材料进行各种理论方法与实验光谱的全面比较研究。为填补这一空白,我们在此简要概述哪种方法适用于准确预测多种纳米结构材料的振动光谱,并为此制定一系列理论指导原则。为此,考虑了四个不同的案例研究,每个案例研究一个特定的材料方面,即柔性MOF中的呼吸振动、刚性MOF UiO - 66中缺陷的表征、金属卤化物钙钛矿CsPbBr中的非谐振动以及沸石H - SSZ - 13孔道上的客体吸附。对于所有这四种材料,在无客体和无缺陷状态且温度足够低时,静态和动态方法都能产生与实验结果定性相似的光谱。当温度升高时,由于存在非谐声子模式,CsPbBr的简谐近似开始失效。此外,简单的简谐模型对缺陷和客体物种的光谱指纹预测不够准确。这两种现象都会使势能面(PES)变平,从而促进亚稳态之间的跃迁,因此需要动态采样。基于本综述中处理的四个案例研究,我们可以提出以下用于模拟功能固态材料准确振动光谱的理论指导原则:(i)对于低温下的纳米结构晶体骨架材料,可以使用基于PES上的几个点和一组独立简谐振子的静态方法来深入了解晶格动力学。(ii)当在较高温度下评估材料或系统中引入额外的复杂性,如强非谐性、缺陷或客体物种时,简谐区域会失效,需要进行动态采样才能正确预测声子光谱。这些指导原则及其对原型材料类别的说明有助于实验和理论研究人员增进从晶格动力学研究中获得的知识。