Error bounds on goodness of fit indicators in vibrational circular dichroism spectroscopy

Significance 

Prior information is very useful in selecting the right and accurate models for particular experimental tasks. This is the same case in vibrational circular dichroism analysis where few alternatives to the correct model are available. In choosing the best alternative model, researchers have employed various methods. For instance, researchers have assigned absolute configuration to a stereoisomer with a spectrum that is consistent with the experimental data. This can be is, however, a complicated process, especially in cases involving more than one stereocenters and occurrence of more conformations, where computations are more demanding, or in case of low concentrations leading to spectra with low signal-to-noise ratio. Alternatively, a combination of several techniques and software has been employed to handle vibrational circular dichroism spectra with non-trivial situations. As such, robust goodness-of-fit indicators have been identified as a powerful tool. Unlike the mean square error, they are less biased by outliers and guarantee accurate results.

To this note, a group of Università di Salerno scientists: Professor Guglielmo Monaco, Dr. Giovanni Procida, Dr. Antonia Di Mola, and Professor Antonio Massa together with Prof. Wouter Herrebout from the University of Antwerp used the bootstrap method to estimate the error bound on five different goodness-of-fit indicators for potential use in vibrational spectroscopy analysis. The indicators were specifically employed to assess and select the suitable and best stereochemical model from possible stereoisomer alternatives. The work is currently published in the research journal Chemical Physics Letters.

In brief, the research team applied the bootstrap method together with the model-averaging method. First, the methods were tested on a 3-chloro-1-butyne of known absolute configuration and then adopted to assign the unknown absolute configuration of synthesized 3-methyl-3-nitromethyl-isoindolinone. The selection criterion was based on comparison with mean square errors taking into consideration the error-normalized dimensionless values. A total of five different goodness-of-fit indicators were applied with and without using the model-averaging method.

Results showed that the stereoisomer turned out to have the known absolute configuration as indicated by the obtained error bounds on the goodness-of-fit indicators. The correctness of the assignment of the absolute configuration could be validated through inspection. For all the utilized goodness-of-fit indicators, the difference between the two models was observed to be significant since the bootstrap estimate for standard deviation of the difference of the goodness-of-fit indicators was lower than its corresponding mean. All the goodness-of-fit indicators were capable of selecting the correct absolute configuration. For the unknown configuration a similar agreement between all goodness-of-fit indicators was observed. Furthermore, restricting the analysis to the few intense peaks whose frequencies were computed accurately, also led to same best model. As such, the difference of goodness-of-fit indicator and its corresponding error is a key tool in indicating the correctness of the absolute configuration.

In a nutshell, the Italian study employed goodness-of-fit indicators, together with error bounds obtained by the bootstrap method, to determine the stereochemical model which best matches the experimental vibrational circular dichroism spectrum. The method was applied to a known case as well as for the assignment of an unknown absolute configuration. The study provides useful insights that will be of great help in vibrational circular dichroism spectroscopy analysis in key fields like pharmaceutical industries.

Error bounds on goodness of fit indicators in vibrational circular dichroism spectroscopy - Advances in Engineering
Model selection requires error bounds on target data as well as on the model itself. Error bounds on goodness-of-fit indicators are used effectively in selecting the best stereochemical model in Vibrational Circular Dichroism spectroscopy

About the author

Guglielmo Monaco graduated in Chemistry in 1996 with 110/110 cum laude at the University of Naples “Federico II”, defending a dissertation on high-resolution NMR spectra of stereoregular polymers, under the leadership of Profs. V. Busico and G. Del Re. In 1999 he defended his doctorate thesis in heterogeneous catalysis completed in the group of Prof. P. Corradini. He has therefore obtained a research grant in Salerno to work with Prof. A. Immirzi on the X-ray diffraction by polymer samples. Since 2002 he is an assistant professor in Physical Chemistry at the Department of Chemistry and Biology at the University of Salerno. Starting from 2006, he has collaborated with Prof. R. Zanasi on:

1) calculations of molecular magnetic properties via the quantum-mechanical current density, and 2) electronic and vibrational circular dichroism spectra.

His major achievements in the first field of investigation are the introduction of the altan-molecules, designed to have magnetically decoupled molecular moieties, the generalization of the method to identify the magnetic character of virtual transitions within the ipsocentric Steiner-Fowler method, the introduction of the anisotropy of the asymmetric current density tensor (AACID), the investigation of stabilizing yet debated H—H interaction in congested molecules.

As for vibrational circular dichroism spectra, he introduced in 2017 the model-averaging technique, which considers the uncertainties of the model itself before setting up comparisons with experimental spectra also affected by the experimental error. The technique has been successfully applied to vibrational circular dichroism spectra in several scientific papers.

Reference

Monaco, G., Procida, G., Di Mola, A., Herrebout, W., & Massa, A. (2020). Error bounds on goodness of fit indicators in vibrational circular dichroism spectroscopy. Chemical Physics Letters, 739, 137000.

Go To Chemical Physics Letters

Check Also

A DFT-based simulated annealing method for the optimization of global energy in zeolite framework systems: Application to natrolite, chabazite and clinoptilolite - Advances in Engineering

A DFT-based simulated annealing method for the optimization of global energy in zeolite framework systems