Spectra-structure correlations of saturated and unsaturated medium-chain fatty acids. Near-infrared and anharmonic DFT study of hexanoic acid and sorbic acid

Development of a CH2-dependent analytical method using near-infrared spectroscopy via the integration of two algorithms: non-dominated sorting genetic-II and competitive adaptive reweighted sampling (NSGAII-CARS)

In most of the near-infrared studies, near-infrared spectra (NIRS) were often mathematically treated. However, these algorithms selected a large number of variables and latent variables, and they caused the over-fitting phenomenon, which became very common. The large number of variables made it impossible to extract the "chemical information" directly from the NIRS. To build robust and interpretable mathematical models, the non-dominated sorting genetic-II-competitive adaptive reweighted sampling (NSGAII-CARS) algorithm was proposed to determine influential functional groups for quantitative analysis. In this research, data on a primary mixture of two amino acids (AAs), namely NH₂(CH₂)₃COOH and HOOC(NH₂)CH(CH₂)₂COOH, was used to illustrate the algorithm. The principle of the algorithm was first to find out the different characteristic spectral regions of two amino acids by extreme points according to Non-dominated Sorting Genetic-II (NSGAII). Second, based on the absolute value of the regression coefficient, we found out [ν(CH₂) + 2δ(CH₂)] and [2ν(CH₂)], where the wavenumber ranged from 6165 to 5683 cm^-1, were the influential functional groups for quantitative analysis. Finally, the CARS (competitive adaptive reweighted sampling) algorithm was combined with NSGAII to find the specific fingerprint points for the determination of two AAs. Compared with the previous results, the NSGAII-CARS algorithm not only pointed out the influential quantitative functional groups but also used only 6 points for HOOC(NH₂)CH(CH₂)₂COOH and 18 points for NH₂(CH₂)₃COOH to achieve the full-spectrum quantitative effect. The results proposed a general algorithm for the quantitative analysis of NIRS obtained in the binary or ternary mixed systems. The MATLAB codes of the NSGAII-CARS algorithm are available on the website: https://github.com/Mark1988NK/NSGAII-CARS-Algorithm.git.

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.

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.

Current and future research directions in computer-aided near-infrared spectroscopy: A perspective

The present review aims to draw a perspective on the vibrational spectroscopy combined with the tools of computational chemistry. This includes an overview of the accomplishments made so far, the assessment of the present development trends and the prospects for continuing these advances. State-of-the-art methods, current challenges and the expected future advances are evaluated from the point-of-view of the practical application in vibrational spectroscopy. A special attention is given to near-infrared (NIR) spectroscopy, which occupies a distinct position among the techniques of vibrational spectroscopy. As the result of intrinsically complex spectra, reliance on the anharmonicity as well as keen interest given to complex materials, NIR spectroscopy may particularly benefit from computational chemistry. The present key limitations hindering development of NIR spectroscopy are identified; these constitute primarily the limit in the treatable system size and the inability to effectively include chemical matrix effects. Given the expanding role of NIR spectroscopy in science and industry, lifting these limitations would directly enhance the general potential of this technique.

Use of Near-Infrared-Mid-Infrared Dual-Wavelength Spectrometry to Obtain Two-Dimensional Difference Spectra of Sesame Oil as Inactive Drug Ingredient

The present study has investigated the transformation of sesame oil kept at low temperature during a definite period of time for refinement (called winterization) as an inactive drug ingredient by using two-dimensional difference spectra (2D-DS) analysis of spectra collected using a near-infrared (NIR) and mid-infrared (MIR) dual-wavelength spectrometer (NIR-MIR-DWS). The NIR and MIR spectra were measured nearly simultaneously from samples of sesame oil before and after winterization. The difference spectrum analysis of the obtained NIR-MIR data elucidated that, after the winterization process, the absorbances at peaks attributed to C=O, C=C, and OH groups decrease while the absorbances arising from the main chain (CH₂) increase. The result indicated the removal of lignan and the fatty acids with relatively short main chains. Moreover, sesame oil unwinterized was cooled from room temperature to near 1 ℃ and subsequently warmed to room temperature. And the cycle was repeated two times. Real-time monitoring during the cooling and warming processes were carried out using the NIR-MIR-DWS. The prediction results obtained from partial least square calibration model for the temperature suggests that there are subtle differences in the oil composition between the first cooling process and after the warming and cooling cycle. For the more detailed analysis, the 2D-DS method is proposed. The results of the analyses using 2D-DS revealed that the starting point of the transformation is around 15 ℃. It can be estimated that sesame oil is mainly transformed by the first cooling down. Moreover, it was implied that the structure of methylene (CH₂) was significantly related to the modifications in sesame oil with temperature change. A series of experimental results elucidated that the winterization of sesame oil removed its impurities and stabilized its conditions. These results are probably the first report on the effect of the winterization process on sesame oil.

Near-Infrared Spectroscopy in Bio-Applications

Near-infrared (NIR) spectroscopy occupies a specific spot across the field of bioscience and related disciplines. Its characteristics and application potential differs from infrared (IR) or Raman spectroscopy. This vibrational spectroscopy technique elucidates molecular information from the examined sample by measuring absorption bands resulting from overtones and combination excitations. Recent decades brought significant progress in the instrumentation (e.g., miniaturized spectrometers) and spectral analysis methods (e.g., spectral image processing and analysis, quantum chemical calculation of NIR spectra), which made notable impact on its applicability. This review aims to present NIR spectroscopy as a matured technique, yet with great potential for further advances in several directions throughout broadly understood bio-applications. Its practical value is critically assessed and compared with competing techniques. Attention is given to link the bio-application potential of NIR spectroscopy with its fundamental characteristics and principal features of NIR spectra.

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.

The solvent effects on dimethyl phthalate investigated by FTIR characterization, solvent parameter correlation and DFT computation

This study set out with the aim of investigating the solvent effects on dimethyl phthalate (DMP) using FTIR characterization, solvent parameter correlation and DFT calculation. DMP exposed to 17 organic solvents manifested varying shift in the carbonyl stretching vibration frequency (ν_CO). Non-alkanols induced Band I and alkanols produced Band I and Band II. Through correlating the ν_CO with the empirical solvent scales including acceptor parameter (AN), Schleyer's linear free energy parameter (G), and linear free salvation energy relationships (LSER), Band I was mainly ascribed to non-specific effects from either non-alkanols or alkanol polymers ((alkanol)_n). ν_CO of the latter indicated minor red shift and less variability compared to the former. An assumption was made and validated about the sequestering of hydroxyl group by the bulky hydrophobic chain in (alkanol)_n, creating what we refer to as "screening effects". Ab initio calculation, on the other hand, provided insights for possible hydrogen binding between DMP and (ethanol)_n or between ethanol monomers. The two components of Band I observed in inert solvents were assigned to the two CO groups adopting differentiated conformations. This in turn prompted our consideration that hydrogen binding was highly selective in favor of lowly associated (alkanol)_n and the particular CO group having relatively less steric hindrance and stronger electron-donating capacity. Band II was therefore believed to derive from hydrogen-bond interactions mainly in manner of 1:1 and 1:2 DMP-(alkanol)_n complexes.

Non-linearity correction in NIR absorption spectra by grouping modeling according to the content of analyte

To correct the non-linearity caused by light scattering in quantitative analysis with near infrared absorption spectra, a new modeling analysis method was proposed: grouping modeling according to the content of analyte. In this study, we tested the proposed method for non-invasive detection of human hemoglobin (Hb) based on dynamic spectrum (DS). We compared the prediction performance of the proposed method with non-grouping modeling method. Experimental results showed that the root mean square error of the prediction set (RMSEP) by the proposed method was reduced by 9.96% and relative standard deviation of the prediction set (RSDP) was reduced by 4.73%. The results demonstrated that the proposed method could reduce the effects of non-linearity on the composition analysis by spectroscopy. This research provides a new method for correcting the non-linearity stemming from light scattering. And the proposed method will accelerate the pace of non-invasive detection of blood components into clinical application.

Anharmonic calculations of frequencies and intensities of OH stretching vibrations of (R)-1,3-butanediol conformers in the fundamentals and first overtones by density functional theory

The frequencies and absorption intensities of the five kinds of conformers of 1,3-butanediol with the same carbon skeleton (GG') were calculated by anharmonic calculation for the fundamentals and first overtones of OH stretching vibrations. The four kinds of conformers form intramolecular hydrogen bonds and one conformer did not. Intramolecular hydrogen bond formation shifted the frequency of fundamental and first overtone of H-bonding OH stretching vibration to the lower frequency. The absorption intensities of the fundamentals as well as the vibrational anharmonicities increased upon hydrogen bond formation, while the intensities of first overtones decreased. The differences of conformers were clearly seen in the frequencies of the first overtones of free OH.

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.

Numerical modeling of polymorphic transformation of oleic acid via near-infrared spectroscopy and factor analysis

Near-infrared (NIR) spectroscopy as a tool for direct and quantitatively screening the minute polymorphic transitions of bioactive fatty acids was assessed basing on a thermal heating process of oleic acid. Temperature-dependent NIR spectral profiles indicate that dynamical variances of COOH group dominate its γ → α phase transition, while the transition from active α to β phase mainly relates to the conformational transfer of acyl chain. Through operating multivariate curve resolution-alternating least squares with factor analysis, instantaneous contribution of each active polymorph during the transition process was illustrated for displaying the progressive evolutions of functional groups. Calculated contributions reveal that the α phase of oleic acid initially is present at around -18 °C, but sharply grows up around -2.2 °C from the transformation of γ phase and finally disappears at the melting point. On the other hand, the β phase of oleic acid is sole self-generation after melt even it embryonically appears at -2.2 °C. Such mathematical approach based on NIR spectroscopy and factor analysis calculation provides a volatile strategy in quantitatively exploring the transition processes of bioactive fatty acids; meanwhile, it maintains promising possibility for instantaneous quantifying each active polymorph of lipid materials.

Inferring anatomical therapeutic chemical (ATC) class of drugs using shortest path and random walk with restart algorithms

The anatomical therapeutic chemical (ATC) classification system is a widely accepted drug classification scheme. This system comprises five levels and includes several classes in each level. Drugs are classified into classes according to their therapeutic effects and characteristics. The first level includes 14 main classes. In this study, we proposed two network-based models to infer novel potential chemicals deemed to belong in the first level of ATC classification. To build these models, two large chemical networks were constructed using the chemical-chemical interaction information retrieved from the Search Tool for Interactions of Chemicals (STITCH). Two classic network algorithms, shortest path (SP) and random walk with restart (RWR) algorithms, were executed on the corresponding network to mine novel chemicals for each ATC class using the validated drugs in a class as seed nodes. Then, the obtained chemicals yielded by these two algorithms were further evaluated by a permutation test and an association test. The former can exclude chemicals produced by the structure of the network, i.e., false positive discoveries. By contrast, the latter identifies the most important chemicals that have strong associations with the ATC class. Comparisons indicated that the two models can provide quite dissimilar results, suggesting that the results yielded by one model can be essential supplements for those obtained by the other model. In addition, several representative inferred chemicals were analyzed to confirm the reliability of the results generated by the two models. This article is part of a Special Issue entitled: Accelerating Precision Medicine through Genetic and Genomic Big Data Analysis edited by Yudong Cai & Tao Huang.