Analytic calculations of anharmonic infrared and Raman vibrational spectra

Harmonic and anharmonic vibrational computations for biomolecular building blocks: Benchmarking DFT and basis sets by theoretical and experimental IR spectrum of glycine conformers

Advanced vibrational spectroscopic experiments have reached a level of sophistication that can only be matched by numerical simulations in order to provide an unequivocal analysis, a crucial step to understand the structure-function relationship of biomolecules. While density functional theory (DFT) has become the standard method when targeting medium-size or larger systems, the problem of its reliability and accuracy are well-known and have been abundantly documented. To establish a reliable computational protocol, especially when accuracy is critical, a tailored benchmark is usually required. This is generally done over a short list of known candidates, with the basis set often fixed a priori. In this work, we present a systematic study of the performance of DFT-based hybrid and double-hybrid functionals in the prediction of vibrational energies and infrared intensities at the harmonic level and beyond, considering anharmonic effects through vibrational perturbation theory at the second order. The study is performed for the six-lowest energy glycine conformers, utilizing available "state-of-the-art" accurate theoretical and experimental data as reference. Focusing on the most intense fundamental vibrations in the mid-infrared range of glycine conformers, the role of the basis sets is also investigated considering the balance between computational cost and accuracy. Targeting larger systems, a broad range of hybrid schemes with different computational costs is also tested.

Spectroscopic and optical properties of 1,2,4-triazolo[4,3-a]pyridin-3(2H)-one as a component of herbicides

The herbicides azafenidin [(2-(2,4-dichloro-5-prop-2-ynoxyphenyl)-5,6,7,8-tetrahydro-[1,2,4]triazolo[4,3-a]pyridin-3-one)] and flumetsulam [(N-(2,6-difluorophenyl)-5-methyl-[1,2,4]triazolo[1,5-a]pyrimidine-2-sulfonamide)] were subjected to IR, Raman, UV-Vis and emission studies. As triazolopyridine is the most prominent and active component of these herbicides, this molecule was characterised by XRD studies, FTIR, Raman, UV-Vis and emission spectra. The experimental data were compared to the results of the DFT quantum chemical calculations carried out for its optimised structure, IR intensities and Raman activities, HOMO-LUMO transitions, and energies of the singlet and triplet states. The characteristics for triazolopyridine quantities were used in the analysis of the studied herbicides.

Electronic State Spectroscopy of Nitromethane and Nitroethane

High-resolution photoabsorption cross-sections in the 3.7-10.8 eV energy range are reinvestigated for nitromethane (CH₃NO₂), while for nitroethane (C₂H₅NO₂), they are reported for the first time. New absorption features are observed for both molecules which have been assigned to vibronic excitations of valence, Rydberg, and mixed valence-Rydberg characters. In comparison with nitromethane, nitroethane shows mainly broad absorption bands with diffuse structures, which can be interpreted as a result of the side-chain effect contributing to an increased number of internal degrees of freedom. New theoretical quantum chemical calculations performed at the time-dependent density functional theory (TD-DFT) level were used to qualitatively help interpret the recorded photoabsorption spectra. From the photoabsorption cross-sections, photolysis lifetimes in the terrestrial atmosphere have been obtained for both compounds. Relevant internal conversion from Rydberg to valence character is noted for both molecules, while the nuclear dynamics of CH₃NO₂ and C₂H₅NO₂ along the C-N reaction coordinate have been evaluated through potential energy curves at the TD-DFT level of theory, showing that the pre-dissociative character is more prevalent in nitromethane than in nitroethane.

Calculation of absolute Raman scattering cross-sections using vibrational self-consistent field/vibrational configuration interaction wave functions

In the present study, the differential scattering cross-sections, depolarization ratios and Raman shifts of small molecular systems are obtained from configuration iteration wave functions of vibrational self-consistent field (VSCF) states. The transition polarizabilities were modeled using the Placzek approximation, neglecting those contributions not arising from the electric dipole mechanism. This theoretical approach is considered a good approximation for samples that absorb in the UV range if the excitation radiation falls in the visible region, as is the case of the molecules selected for the present study, namely: water, methane, and acetylene. Potential energy and electronic polarizability surfaces are calculated by the CCSD(T) and CC3 methods with aug-cc-p(C)V(T,Q,5)Z basis sets. The vibrational Hamiltonian includes the vibrational angular momentum contribution of the Watson kinetic energy operator. As expected, due to the variational nature of the VSCF and vibrational configuration interaction (VCI) methods, the Raman transition wavenumbers are substantially improved over the harmonic predictions. Surprisingly, the scattering cross-sections obtained using the harmonic approximation or the VSCF method better agrees with the experimental values than those cross-sections predicted using VCI wave functions. The more significant deviations of the VCI results from the experimental reference may be related to the significant uncertainties of the measured cross-sections. Still, it may also indicate that the VCI Raman transition moments may require a more accurate description of the electronic polarizability surface. Finally, the depolarization ratios calculated for H₂ O and C₂ D₂ using harmonic and VCI wave functions have similar accuracy, whereas, for C₂ H₂ and C₂ HD, the VCI results are more accurate.

Electromagnetic bioeffects: a multiscale molecular simulation perspective

Electromagnetic bioeffects remain an enigma from both the experimental and theoretical perspectives despite the ubiquitous presence of related technologies in contemporary life. Multiscale computational modelling can provide valuable insights into biochemical systems and predict how they will be perturbed by external stimuli. At a microscopic level, it can be used to determine what (sub)molecular scale reactions various stimuli might induce; at a macroscopic level, it can be used to examine how these changes affect dynamic behaviour of essential molecules within the crowded biomolecular milieu in living tissues. In this review, we summarise and evaluate recent computational studies that examined the impact of externally applied electric and electromagnetic fields on biologically relevant molecular systems. First, we briefly outline the various methodological approaches that have been employed to study static and oscillating field effects across different time and length scales. The practical value of such modelling is then illustrated through representative case-studies that showcase the diverse effects of electric and electromagnetic field on the main physiological solvent - water, and the essential biomolecules - DNA, proteins, lipids, as well as some novel biomedically relevant nanomaterials. The implications and relevance of the theoretical multiscale modelling to practical applications in therapeutic medicine are also discussed. Finally, we summarise ongoing challenges and potential opportunities for theoretical modelling to advance the current understanding of electromagnetic bioeffects for their modulation and/or beneficial exploitation in biomedicine and industry.

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.

AC/DC Analysis: Broad and Comprehensive Approach to Analyze Infrared Intensities at the Atomic Level

We present a complete theoretical protocol to partition infrared intensities into terms owing to individual atoms by two different but related approaches: the atomic contributions (ACs) show how the entire molecular vibrational motion affects the electronic structure of a single atom and the total infrared intensity. On the other hand, the dynamic contributions (DCs) show how the displacement of a single atom alters the electronic structure of the entire molecule and the total intensity. The two analyses are complementary ways of partitioning the same total intensity and conserve most of the features of the total intensity itself. Combined, they are called the AC/DC analysis. These can be further partitioned following the CCTDP (or CCT) models according to the population analysis chosen by the researcher. The main conceptual features of the equations are highlighted, and representative numerical results are shown to support the interpretation of the equations. The results are invariant to rotation and translation and can readily be extended to molecules of any size, shape, or symmetry. Although the AC/DC analysis requires the choice of a charge model, all charge models that correctly reproduce the total molecular dipole moment can be used. A fully automated protocol managed by the Placzek program is made available, free of charge and with input examples.

Complete analytic anharmonic hyper-Raman scattering spectra

We present the first computational treatment of the complete second-order vibrational perturbation theory applied to hyper-Raman scattering spectroscopy. The required molecular properties are calculated in a fully analytic manner using a recently developed program [Ringholm, Jonsson and Ruud, J. Comp. Chem., 2014, 35, 622] that utilizes recursive routines. For some of the properties, these calculations are the first analytic calculations of their kind at their respective levels of theory. We apply this approach to the calculation of the hyper-Raman spectra of methane, ethane and ethylene and compare these to available experimental data. We show that the anharmonic corrections have a larger effect on the vibrational frequencies than on the spectral intensities, but that the inclusion of combination and overtone bands in the anharmonic treatment can improve the agreement with the experimental data, although the quality of available experimental data limits a detailed comparison.

Cyanoform and Its Isomers. Relative Stabilities, Spectroscopic Features, and Rearrangements by Coupled Cluster and MCSCF-Based Methods

Although an isolation of elusive tricyanomethane HC(CN)₃ was recently reported, the existence of other HC₄N₃ species has yet to be confirmed. In this work, the relative stabilities, spectroscopic features, and rearrangements of tricyanomethane and its four isomers are examined using single- (CCSD(T), CCSD(T)-F12) and multireference (MCSCF, MRPT2) methods. Tricyanomethane and dicyanoketenimine (NC)₂CCNH, which are found to be the two most stable HC₄N₃ isomers lying within 9 kcal/mol, can be discriminated by their spectroscopic parameters. The predicted stepwise interconversion path relating HC(CN)₃ and (NC)₂CCNH features the HC₄N₃ species comprising the C-C-N ring moiety, with the largest barrier being associated with the initial H migration to one of the CN carbons. Adding a water molecule reduces the H migration barrier strongly and makes it possible to interconvert tricyanomethane to dicyanoketenimine in a "concerted" way.

On the benefits of localized modes in anharmonic vibrational calculations for small molecules

Anharmonic vibrational calculations can already be computationally demanding for relatively small molecules. The main bottlenecks lie in the construction of the potential energy surface and in the size of the excitation space in the vibrational configuration interaction (VCI) calculations. To address these challenges, we use localized-mode coordinates to construct potential energy surfaces and perform vibrational self-consistent field and L-VCI calculations [P. T. Panek and C. R. Jacob, ChemPhysChem 15, 3365 (2014)] for all vibrational modes of two prototypical test cases, the ethene and furan molecules. We find that the mutual coupling between modes is reduced when switching from normal-mode coordinates to localized-mode coordinates. When using such localized-mode coordinates, we observe a faster convergence of the n-mode expansion of the potential energy surface. This makes it possible to neglect higher-order contributions in the n-mode expansion of the potential energy surface or to approximate higher-order contributions in hybrid potential energy surfaces, which reduced the computational effort for the construction of the anharmonic potential energy surface significantly. Moreover, we find that when using localized-mode coordinates, the convergence with respect to the VCI excitation space proceeds more smoothly and that the error at low orders is reduced significantly. This makes it possible to devise low-cost models for obtaining a first approximation of anharmonic corrections. This demonstrates that the use of localized-mode coordinates can be beneficial already in anharmonic vibrational calculations of small molecules and provides a possible avenue for enabling such accurate calculations also for larger molecules.