Spectroscopic and Quantum Mechanical Calculation Study of the Effect of Isotopic Substitution on NIR Spectra of Methanol

In silico NIR spectroscopy - A review. Molecular fingerprint, interpretation of calibration models, understanding of matrix effects and instrumental difference

Quantum mechanical calculations are routinely used as a major support in mid-infrared (MIR) and Raman spectroscopy. In contrast, practical limitations for long time formed a barrier to developing a similar synergy between near-infrared (NIR) spectroscopy and computational chemistry. Recent advances in theoretical methods suitable for calculation of NIR spectra opened the pathway to modeling NIR spectra of various molecules. Accurate theoretical reproduction of NIR spectra of molecules reaching the size of long-chain fatty acids was accomplished so far. In silico NIR spectroscopy, where the spectra are calculated ab initio, provides substantial improvement in our understanding of the overtones and combination bands that overlap in staggering numbers and create complex lineshape typical for NIR spectra. This improves the comprehension of the spectral information enabling access to rich and detail molecular footprint, essential for fundamental research and useful in routine analysis by NIR spectroscopy and chemometrics. This review article summarizes the most recent accomplishments in the emerging field with examples of simulated NIR spectra of molecules reaching long-chain fatty acids and polymers. In addition to detailed NIR band assignments and new physical insights, simulated spectra enable innovative support in applications. Understanding of the difference in the performance observed between miniaturized NIR spectrometers and chemical interpretation of the chemometric models are noteworthy here. These new elements integrated into NIR spectroscopy framework enable a knowledge-based design of the analysis with comprehension of the processed chemical information.

Spectra-structure correlations in NIR region of polymers from quantum chemical calculations. The cases of aromatic ring, C=O, C≡N and C-Cl functionalities

Near-infrared (NIR) spectroscopy is a valued analytical tool in various applications involving polymers. However, complex nature of NIR spectra imposes difficulties in their direct interpretation. Here, anharmonic quantum chemical calculations are used to simulate NIR spectra of nine polymers; acrylonitrile butadiene styrene (ABS), ethylene-vinyl acetate (EVAC), polycarbonate (PC), polyethylene terephthalate (PET), polylactide or polylactic acid (PLA), polymethylmethacrylate (PMMA), polyoxymethylene (POM), polystyrene (PS) and polyvinylchloride (PVC). The generalized spectra-structure correlations are derived for these systems with focus given to the manifestation in NIR spectra of aromatic ring, C=O, C≡N and C-Cl functionalities. It is concluded that the nature of NIR polymer bands is only moderately sensitive to the remote chemical neighborhood. The majority of NIR absorption of polymers originates from binary combination bands, while the first overtones are meaningful only in ca. 6200-5500 cm^-1 region. The contribution of the overtone bands is relatively higher for the polymers bearing aromatic rings because of higher intensity of C-H stretching overtones. Highly characteristic combination bands of the modes localized in aromatic ring (ring deformation and CH stretching) are relatively independent on the remaining structure of the polymer. The combination bands originating from C=O group are more sensitive to the chemical neighborhood in near proximity, forming a useful fingerprint for a specific polymer. In contrast, the vibrational bands of C≡N functionality are far less useful in NIR region than in infrared (IR) region. With aid of the calculated absorption bands, structural specificity of NIR spectroscopy of polymers can be markedly improved.

Anharmonicity and Spectra-Structure Correlations in MIR and NIR Spectra of Crystalline Menadione (Vitamin K3)

Mid-infrared (MIR) and near-infrared (NIR) spectra of crystalline menadione (vitamin K3) were measured and analyzed with aid of quantum chemical calculations. The calculations were carried out using the harmonic approach for the periodic model of crystal lattice and the anharmonic DVPT2 calculations applied for the single molecule model. The theoretical spectra accurately reconstructed the experimental ones permitting for reliable assignment of the MIR and NIR bands. For the first time, a detailed analysis of the NIR spectrum of a molecular system based on a naphthoquinone moiety was performed to elucidate the relationship between the chemical structure of menadione and the origin of the overtones and combination bands. In addition, the importance of these bands during interpretation of the MIR spectrum was demonstrated. The overtones and combination bands contribute to 46.4% of the total intensity of menadione in the range of 3600-2600 cm-1. Evidently, these bands play a key role in shaping of the C-H stretching region of MIR spectrum. We have shown also that the spectral regions without fundamentals may provide valuable structural information. For example, the theoretical calculations reliably reconstructed numerous overtones and combination bands in the 4000-3600 and 2800-1800 cm-1 ranges. These results, provide a comprehensive origin of the fundamentals, overtones and combination bands in the NIR and MIR spectra of menadione, and the relationship of these spectral features with the molecular structure.

Anharmonic DFT Study of Near-Infrared Spectra of Caffeine: Vibrational Analysis of the Second Overtones and Ternary Combinations

Anharmonic quantum chemical calculations were employed to simulate and interpret a near-infrared (NIR) spectrum of caffeine. First and second overtones, as well as binary and ternary combination bands, were obtained, accurately reproducing the lineshape of the experimental spectrum in the region of 10,000–4000 cm⁻¹ (1000–2500 nm). The calculations enabled performing a detailed analysis of NIR spectra of caffeine, including weak bands due to the second overtones and ternary combinations. A highly convoluted nature of NIR spectrum of caffeine was unveiled, with numerous overlapping bands found beneath the observed spectral lineshape. To properly reflect that intrinsic complexity, the band assignments were provided in the form of heat maps presenting the contributions to the NIR spectrum from various kinds of vibrational transitions. These contributions were also quantitatively assessed in terms of the integral intensities. It was found that the combination bands provide the decisively dominant contributions to the NIR spectrum of caffeine. The first overtones gain significant importance between 6500–5500 cm⁻¹, while the second overtones are meaningful in the higher wavenumber regions, particularly in the 10,000–7000 cm⁻¹ region. The obtained detailed band assignments enabled deep interpretation of the absorption regions of caffeine identified in the literature as meaningful for analytical applications of NIR spectroscopy focused on quantitative analysis of caffeine content in drugs and natural products.

Advances, challenges and perspectives of quantum chemical approaches in molecular spectroscopy of the condensed phase

The purpose of this review is to demonstrate advances, challenges and perspectives of quantum chemical approaches in molecular spectroscopy of the condensed phase. Molecular spectroscopy, particularly vibrational spectroscopy and electronic spectroscopy, has been used extensively for a wide range of areas of chemical sciences and materials science as well as nano- and biosciences because it provides valuable information about structure, functions, and reactions of molecules. In the meantime, quantum chemical approaches play crucial roles in the spectral analysis. They also yield important knowledge about molecular and electronic structures as well as electronic transitions. The combination of spectroscopic approaches and quantum chemical calculations is a powerful tool for science, in general. Thus, our article, which treats various spectroscopy and quantum chemical approaches, should have strong implications in the wider scientific community. This review covers a wide area of molecular spectroscopy from far-ultraviolet (FUV, 120-200 nm) to far-infrared (FIR, 400-10 cm-1)/terahertz and Raman spectroscopy. As quantum chemical approaches, we introduce several anharmonic approaches such as vibrational self-consistent field (VSCF) and the combination of periodic harmonic calculations with anharmonic corrections based on finite models, grid-based techniques like the Numerov approach, the Cartesian coordinate tensor transfer (CCT) method, Symmetry-Adapted Cluster Configuration-Interaction (SAC-CI), and the ZINDO (Semi-empirical calculations at Zerner's Intermediate Neglect of Differential Overlap). One can use anharmonic approaches and grid-based approaches for both infrared (IR) and near-infrared (NIR) spectroscopy, while CCT methods are employed for Raman, Raman optical activity (ROA), FIR/terahertz and low-frequency Raman spectroscopy. Therefore, this review overviews cross relations between molecular spectroscopy and quantum chemical approaches, and provides various kinds of close-reality advanced spectral simulation for condensed phases.

Theoretical Simulation of Near-Infrared Spectrum of Piperine: Insight into Band Origins and the Features of Regression Models

We investigated the near-infrared spectrum of piperine using quantum mechanical calculations. We evaluated two efficient approaches, DVPT2//PM6 and DVPT2//ONIOM [PM6:B3LYP/6-311++G(2df, 2pd)] that yielded a simulated spectrum with varying accuracy versus computing time factor. We performed vibrational assignments and unveiled complex nature of the near-infrared spectrum of piperine, resulting from a high level of band convolution. The most meaningful contribution to the near-infrared absorption of piperine results from binary combination bands. With the available detailed near-infrared assignment of piperine, we interpreted the properties of partial least square regression models constructed in our earlier study to describe the piperine content in black pepper samples. Two models were compared with spectral data sets obtained with a benchtop and a miniaturized spectrometer. The two spectrometers implement distinct technology which leads to a profound instrumental difference and discrepancy in the predictive performance when analyzing piperine content. We concluded that the sensitivity of the two instruments to certain types of piperine vibrations is different and that the benchtop spectrometer unveiled higher selectivity. Such difference in obtaining chemical information from a sample can be one of the reasons why the benchtop spectrometer performs better in analyzing the piperine content of black pepper. This evidenced direct correspondence between the features critical for applied near-infrared spectroscopic routine and the underlying vibrational properties of the analyzed constituent in a complex sample.

How much anharmonicity is in vibrational spectra of CH3I and CD3I?

This work presents new experimental and theoretical insights on vibrational spectra of CH₃I and CD₃I in the liquid phase. For the first time, we provided the contributions from different vibrational modes to mid-infrared (MIR) and near-infrared (NIR) spectra and estimated the extent of anharmonicity in the MIR region. Direct comparison of the intensities from ATR-IR and NIR transmission spectra was possible due to normalization of ATR-IR spectra. As a reference for normalization, we applied the area of the ν_s(CH₃)/ν_s(CD₃) band recorded in transmission mode. Our results show that the corresponding vibrational modes of CH₃I and CD₃I have similar contributions to the total intensity (MIR + NIR), however, these contributions are distributed in a different way between MIR and NIR regions. As expected, most of intensity in MIR spectra originates from the fundamental transitions (>90%). The fundamental bands together with the first overtones and the binary combinations contribute to more than 99% of MIR intensity for both compounds. Therefore, reliable reconstruction of MIR spectra can be achieved by considering only these vibrational modes. On the other hand, accurate simulation of NIR spectra requires including the higher-order transitions. In the case of CD₃I, the fourth-order transitions contribute to 12.7% of NIR intensity. The contributions from NIR region are significantly smaller than those from MIR range and were estimated to be 6.7% for CH₃I and 2.3% for CD₃I. The theoretical calculations provide a reasonable estimation of the total contribution from the fundamental bands. Yet, the calculated contributions from the anharmonic transitions are different from those obtained from the experimental data. MIR spectra of CH₃I and CD₃I reveal an unexpected increase in the intensity of some overtones and combination bands indicating the presence of Fermi resonances. These resonances are responsible for differences in contributions from the first overtones and binary combinations between CH₃I and CD₃I.

The coupling of localised, vibrational modes - Probing OH-bands of organic molecules via a two dimensional Numerov approach

In this article the extension of the grid-based Numerov approach to probe two coupled, localised vibrational modes is assessed. The theoretically obtained wave numbers are compared to experimental results for five increasingly complex organic molecules carrying two OH groups measured in gas-phase as well as carbon tetrachloride. By using an appropriate spacing of the associated potential energy grid a deviation of the predicted wave numbers with experiment of ≤1% is achieved for both the fundamental and the first overtone bands. In particular the calculated wave numbers of aliphatic species in vacuum underline the versatility of this approach. In addition, it is demonstrated that bicubic interpolation is a viable strategy to greatly reduce the required data points and thus, the computational effort. Comparison of predicted wave numbers obtained for different conformers with experimental data enables the identification of the most relevant conformer present in solution. Since especially the accurate calculation of overtone vibrations is known to be challenging in case of strongly anharmonic potentials such as OH bonds, the presented approach provides a particularly efficient route to study the properties of the associated overtone contribution under the influence of inter-mode coupling. This is due to the fact that the Numerov approach requires no assumption about form and composition of the vibrational wave functions. In addition, the presented method also provides one of the simplest routes to access combined excitations of the considered vibrational modes.

Solvation effects on wavenumbers and absorption intensities of the OH-stretch vibration in phenolic compounds – electrical- and mechanical anharmonicity via a combined DFT/Numerov approach

A previously measured oscillatory intensity pattern in phenolic compounds between different solvents was successfully reproduced for the first time, employing modern grid-based methods to solve the time-independent vibrational Schrödinger equation.

Simulated NIR spectra as sensitive markers of the structure and interactions in nucleobases

Near-infrared (near-IR; NIR) spectroscopy is continuously advancing in biophysical and biochemical fields of investigation. For instance, recent progresses in NIR hyperspectral imaging of biological systems may be noted. However, interpretation of NIR bands for biological samples is difficult and creates a considerable barrier in exploring the full potential of NIR spectroscopy in bioscience. For this reason, we carried out a systematic study of NIR spectra of adenine, cytosine, guanine, and thymine in polycrystalline state. Interpretation of NIR spectra of these nucleobases was supported by anharmonic vibrational analysis using Deperturbed Vibrational Second-Order Perturbation Theory (DVPT2). A number of molecular models of nucleobases was applied to study the effect of the inter-molecular interactions on the NIR spectra. The accuracy of simulated NIR spectra appears to depend on the intra-layer interactions; in contrast, the inter-layer interactions are less influential. The best results were achieved by combining the simulated spectra of monomers and dimers. It is of particular note that in-plane deformation bands are far more populated than out-of-plane ones and the importance of ring modes is relatively small. This trend is in contrast to that observed in mid-IR region. As shown, the local, short-range chemical neighborhood of nucleobase molecules influence their NIR spectra more considerably. This suggests that NIR spectra are more sensitive probe of the nucleobase pairing than mid-IR ones. The obtained results allow, for the first time, to construct a frequency correlation table for NIR spectra of purines and pyrimidines.

Spectra-Structure Correlations in Isotopomers of Ethanol (CX3CX2OX; X = H, D): Combined Near-Infrared and Anharmonic Computational Study

The effect of isotopic substitution on near-infrared (NIR) spectra has not been studied in detail. With an exception of few major bands, it is difficult to follow the spectral changes due to complexity of NIR spectra. Recent progress in anharmonic quantum mechanical calculations allows for accurate reconstruction of NIR spectra. Taking this opportunity, we carried out a systematic study of NIR spectra of six isotopomers of ethanol (CX₃CX₂OX; X = H, D). Besides, we calculated the theoretical spectra of two other isotopomers (CH₃CD₂OD and CD₃CH₂OD) for which the experimental spectra are not available. The anharmonic calculations were based on generalized vibrational second-order perturbation theory (GVPT2) at DFT and MP2 levels with several basis sets. We compared the accuracy and efficiency of various computational methods. It appears that the best results were obtained with B2PLYP-GD3BJ/def2-TZVP//CPCM approach. Our simulations included the first and second overtones, as well as binary and ternary combinations bands. This way, we reliably reproduced even minor bands in the spectra of diluted samples (0.1 M in CCl₄). On this basis, the effect of isotopic substitution on NIR spectra of ethanol was accurately reproduced and comprehensively explained.

Distinct Difference in Sensitivity of NIR vs. IR Bands of Melamine to Inter-Molecular Interactions with Impact on Analytical Spectroscopy Explained by Anharmonic Quantum Mechanical Study

Melamine (IUPAC: 1,3,5-Triazine-2,4,6-triamine) attracts high attention in analytical vibrational spectroscopy due to its misuse as a food adulterant. Vibrational spectroscopy [infrared (IR) and Raman and near-infrared (NIR) spectroscopy] is a major quality control tool in the detection and quantification of melamine content. The physical background for the measured spectra is not interpreted in analytical spectroscopy using chemometrics. In contrast, quantum mechanical calculations are capable of providing deep and independent insights therein. So far, the NIR region of crystalline melamine has not been studied by quantum mechanical calculations, while the investigations of its IR spectra have remained limited. In the present work, we employed fully anharmonic calculation of the NIR spectrum of melamine based on finite models, and also performed IR spectral simulation by using an infinite crystal model—periodic in three dimensions. This yielded detailed and unambiguous NIR band assignments and revised the previously known IR band assignments. We found that the out-of-plane fundamental transitions, which are essential in the IR region, are markedly more sensitive to out-of-plane inter-molecular interactions of melamine than NIR transitions. Proper description of the chemical surrounding of the molecule of melamine is more important than the anharmonicity of its vibrations. In contrast, the NIR bands mostly arise from in-plane vibrations, and remain surprisingly insensitive to the chemical environment. These findings explain previous observations that were reported in IR and NIR analytical studies of melamine.

Breakthrough Potential in Near-Infrared Spectroscopy: Spectra Simulation. A Review of Recent Developments

Near-infrared (12,500–4,000 cm⁻¹; 800–2,500 nm) spectroscopy is the hallmark for one of the most rapidly advancing analytical techniques over the last few decades. Although it is mainly recognized as an analytical tool, near-infrared spectroscopy has also contributed significantly to physical chemistry, e.g., by delivering invaluable data on the anharmonic nature of molecular vibrations or peculiarities of intermolecular interactions. In all these contexts, a major barrier in the form of an intrinsic complexity of near-infrared spectra has been encountered. A large number of overlapping vibrational contributions influenced by anharmonic effects create complex patterns of spectral dependencies, in many cases hindering our comprehension of near-infrared spectra. Quantum mechanical calculations commonly serve as a major support to infrared and Raman studies; conversely, near-infrared spectroscopy has long been hindered in this regard due to practical limitations. Advances in anharmonic theories in hyphenation with ever-growing computer technology have enabled feasible theoretical near-infrared spectroscopy in recent times. Accordingly, a growing number of quantum mechanical investigations aimed at near-infrared region has been witnessed. The present review article summarizes these most recent accomplishments in the emerging field. Applications of generalized approaches, such as vibrational self-consistent field and vibrational second order perturbation theories as well as their derivatives, and dense grid-based studies of vibrational potential, are overviewed. Basic and applied studies are discussed, with special attention paid to the ones which aim at improving analytical spectroscopy. A remarkable potential arises from the growing applicability of anharmonic computations to solving the problems which arise in both basic and analytical near-infrared spectroscopy. This review highlights an increased value of quantum mechanical calculations to near-infrared spectroscopy in relation to other kinds of vibrational spectroscopy.

Probing vibrational coupling via a grid-based quantum approach-an efficient strategy for accurate calculations of localized normal modes in solid-state systems

In this work an approach to investigate the properties of strongly localized vibrational modes of functional groups in bulk material and on solid-state surfaces is presented. The associated normal mode vectors are approximated solely on the basis of structural information and obtained via diagonalization of a reduced Hessian. The grid-based Numerov procedure in one and two dimensions is then applied to an adequate scan of the respective potential surface yielding the associated vibrational wave functions and energy eigenvalues. This not only provides a detailed description of anharmonic effects but also an accurate inclusion of the coupling between the investigated vibrational states on a quantum mechanical level. All results obtained for the constructed normal modes are benchmarked against their analytical counterparts obtained from the diagonalization of the total Hessian of the entire system. Three increasingly complex systems treated at quantum chemical level of theory have been considered, namely the symmetric and asymmetric stretch vibrations of an isolated water molecule, hydroxyl groups bound to the surface of GeO2 (001), α-quartz(001) and Rutil (001) as well as crystalline Li2 NH serving as an example for a bulk material. While the data obtained for the individual systems verify the applicability of the proposed methodology, comparison to experimental data demonstrates the accuracy of this methodology despite the restriction to limit this methodology to a few selected vibrational modes. The possibility to investigate vibrational phenomena of localized normal modes without the requirement of executing costly harmonic frequency calculations of the entire system enables the application of this method to cases in which the determination of normal modes is prohibitively expensive or not available for a particular level of theory. © 2018 Wiley Periodicals, Inc.

Spectra-structure correlations in NIR region: Spectroscopic and anharmonic DFT study of n-hexanol, cyclohexanol and phenol

We investigated near-infrared (7500-4000 cm^-1) spectra of n-hexanol, cyclohexanol and phenol in CCl₄ (0.2 M) by using anharmonic quantum calculations. These molecules represent three major kinds of alcohols; linear and cyclic aliphatic, and aromatic ones. Vibrational second-order perturbation theory (VPT2) was employed to calculate the first overtones and binary combination modes and to reproduce the experimental NIR spectra. The level of conformational flexibility of these three alcohols varies from one stable conformer of phenol through four conformers of cyclohexanol to few hundreds conformers in the case of n-hexanol. To take into account the most relevant conformational population of n-hexanol, a systematic conformational search was performed. Accurate reproduction of the experimental NIR spectra was achieved and detailed spectra-structure correlations were obtained for these three alcohols. VPT2 approach provides less reliable description of highly anharmonic modes, i.e. OH stretching. In the present work this limitation was manifested in erroneous results yielded by VPT2 for 2νOH mode of cyclohexanol. To study the anharmonicity of this mode we solved the corresponding time-independent Schrödinger equation based on a dense-grid probing of the relevant vibrational potential. These results allowed for significant improvement of the agreement between the calculated and experimental 2νOH band of cyclohexanol. Various important biomolecules include similar structural units to the systems investigated here. A detailed knowledge on spectral properties of these three types of alcohols is therefore essential for advancing our understanding of NIR spectroscopy of biomolecules.