Selective calculation of high-intensity vibrations in molecular resonance Raman spectra

Harmonic Infrared and Raman Spectra in Molecular Environments Using the Polarizable Embedding Model

We present a fully analytic approach to calculate infrared (IR) and Raman spectra of molecules embedded in complex molecular environments modeled using the fragment-based polarizable embedding (PE) model. We provide the theory for the calculation of analytic second-order geometric derivatives of molecular energies and first-order geometric derivatives of electric dipole moments and dipole-dipole polarizabilities within the PE model. The derivatives are implemented using a general open-ended response theory framework, thus allowing for an extension to higher-order derivatives. The embedding-potential parameters used to describe the environment in the PE model are derived through first-principles calculations, thus allowing a wide variety of systems to be modeled, including solvents, proteins, and other large and complex molecular environments. Here, we present proof-of-principle calculations of IR and Raman spectra of acetone in different solvents. This work is an important step toward calculating accurate vibrational spectra of molecules embedded in realistic environments.

Comprehensive Assessment of GFN Tight-Binding and Composite Density Functional Theory Methods for Calculating Gas-Phase Infrared Spectra

Vibrational spectroscopy is a valuable and widely used analytical tool for the characterization of chemical substances. We investigate the performance of semiempirical quantum mechanical GFN tight-binding and force-field methods for the calculation of gas-phase infrared spectra in comparison to experiment and low-cost (B3LYP-3c) density functional theory. A data set of 7247 experimental references was used to evaluate method performance based on automatic spectra comparison. Various quantitative spectral similarity measures were employed for the comparison between theory and experiment and for determining empirical scaling factors. It is shown that the scaling of atomic masses provides an accurate yet simple alternative to standard global frequency scaling in density functional theory (DFT) and semiempirical calculations. Furthermore, the method performance for 58 exemplary transition metal complexes was investigated. The efficient DFT composite method B3LYP-3c, that was introduced in the course of this work, was found to be excellently suited for general IR spectra calculations. The GFN1- and GFN2-xTB tight-binding methods clearly outperformed the PMx competitors. Conformational changes were investigated for a subset of the data and are found to have a mediocre strong influence on the simulated spectra suggesting that the corresponding elaborate sampling steps may be neglected in automated compound identification workflows.

Computing Bulk Phase Resonance Raman Spectra from ab Initio Molecular Dynamics and Real-Time TDDFT

We present our novel approach for computing resonance Raman (RR) spectra of periodic bulk phase systems from ab initio molecular dynamics, including solvent influence and some anharmonic effects. Based on real-time time-dependent density functional theory, we obtain the RR spectra for all laser wavelengths in one pass. We compute the RR spectrum of uracil in aqueous solution, which is in good agreement with experiment. This is the first simulation of a bulk phase RR spectrum.

Frequency Range Selection Method for Vibrational Spectra

Theoretical calculations of vibrational properties are widely used to explain and predict experimental spectra. However, with standard quantum chemical methods all molecular motions are considered, which is rather time-consuming for large molecules. Because typically only a specific spectral region is of experimental interest, we propose here an efficient method that allows calculation of only a selected frequency interval. After a computationally cheap low-level estimate of the molecular motions, the computational time is proportional to the number of normal modes needed to describe this frequency range. Results for a medium-sized molecule show a reduction in computational time of up to 1 order of magnitude with negligible loss in accuracy. We also show that still larger computational savings are possible by using an additional intensity-selection procedure.

Multiconfigurational Effects in Theoretical Resonance Raman Spectra

We analyze resonance Raman spectra of the nucleobase uracil in the short‐time approximation calculated with multiconfigurational methods. We discuss the importance of static electron correlation by means of density‐matrix renormalization group self‐consistent field (DMRG‐SCF) calculations. Our DMRG‐SCF results reveal that a minimal active orbital space that leads to a qualitatively correct description of the resonance Raman spectrum of uracil should encompass parts of the σ/σ* bonding/anti‐bonding orbitals of the pyrimidine ring. We trace these findings back to the considerable entanglement between the σ/σ* bonding/anti‐bonding as well as valence π/π* orbitals in the excited‐state electronic structure of uracil, which indicates non‐negligible non‐dynamical correlation effects that are less pronounced in the electronic ground state.

Anharmonic Theoretical Vibrational Spectroscopy of Polypeptides

Because of the size of polypeptides and proteins, the quantum-chemical prediction of their vibrational spectra presents an exceptionally challenging task. Here, we address one of these challenges, namely, the inclusion of anharmonicities. By performing the expansion of the potential energy surface in localized-mode coordinates instead of the normal-mode coordinates, it becomes possible to calculate anharmonic vibrational spectra of polypeptides efficiently and reliably. We apply this approach to calculate the infrared, Raman, and Raman optical activity spectra of helical alanine polypeptides consisting of up to 20 amino acids. We find that while anharmonicities do not alter the band shapes, simple scaling procedures cannot account for the different shifts found for the individual bands. This closes an important gap in theoretical vibrational spectroscopy by making it possible to quantify the anharmonic contributions and opens the door to a first-principles calculation of multidimensional vibrational spectra.

Vibronic-structure tracking: a shortcut for vibrationally resolved UV/Vis-spectra calculations

The vibrational coarse structure and the band shapes of electronic absorption spectra are often dominated by just a few molecular vibrations. By contrast, the simulation of the vibronic structure even in the simplest theoretical models usually requires the calculation of the entire set of normal modes of vibration. Here, we exploit the idea of the mode-tracking protocol [M. Reiher and J. Neugebauer, J. Chem. Phys. 118, 1634 (2003)] in order to directly target and selectively calculate those normal modes which have the largest effect on the vibronic band shape for a certain electronic excitation. This is achieved by defining a criterion for the importance of a normal mode to the vibrational progressions in the absorption band within the so-called "independent mode, displaced harmonic oscillator" (IMDHO) model. We use this approach for a vibronic-structure investigation for several small test molecules as well as for a comparison of the vibronic absorption spectra of a truncated chlorophyll a model and the full chlorophyll a molecule. We show that the method allows to go beyond the often-used strategy to simulate absorption spectra based on broadened vertical excitation peaks with just a minimum of computational effort, which in case of chlorophyll a corresponds to about 10% of the cost for a full simulation within the IMDHO approach.

Efficient calculation of anharmonic vibrational spectra of large molecules with localized modes

The analysis and interpretation of the vibrational spectra of complex (bio)molecular systems, such as polypeptides and proteins, requires support from quantum-chemical calculations. Such calculations are currently restricted to the harmonic approximation. Here, we show how one of the main bottlenecks in such calculations, the evaluation of the potential energy surface, can be overcome by using localized modes instead of the commonly employed normal modes. We apply such local vibrational self-consistent field (L-VSCF) and vibrational configuration interaction (L-VCI) calculations to a cyclic water tetramer and a helical hexa-alanine peptide. The results show that the use of localized modes is equivalent to the commonly used normal modes, but offers several advantages. First, a faster convergence with respect to the excitation level is observed in L-VCI calculations. Second, the localized modes provide a reduced representation of the couplings between modes that show a regular coupling pattern. This can be used to disregard a significant number of small two-mode potentials a priori. Several such reduced coupling approximations are explored, and we show that the number of single-point calculations required to evaluate the potential energy surface can be significantly reduced without introducing noticeable errors in the resulting vibrational spectra.

PICVib: an accurate, fast, and simple procedure to investigate selected vibrational modes at high theoretical levels

A new approach Procedure for Investigating Categories of Vibrations (PICVib) for estimating vibrational frequencies of selected modes using only the structure and energy calculations at a more demanding computational level is presented and explored. The PICVib has an excellent performance at only a small fraction of the computational demand required for a complete analytical calculation. The errors are smaller than ca. 0.5% when DFT functionals are combined with high level ab initio methods. The approach is general because it can use any quantum chemical program and electronic structure method. It is very robust because it was validated for a wide range of frequency values (ca. 20-4800 cm(-1)) and systems: XH(3) (D(3h) ) with X = B, Al, Ga, N, P, As, O, S, and Se, YH(4) (D(4h) ) with Y = C, Si, and Ge, conformers of RDX, S(N) 2 and E2 reactions, [W(dppe)(2)(NNC(5)H(10))] complex, carbon nanotubes, and hydrogen-bonded complexes including guanine-cytosine pair.