Contribution of high-energy conformations to NMR chemical shifts, a DFT-BOMD study

Abstract

This paper highlights the relevance of including the high-energy conformational states sampled by Born–Oppenheimer molecular dynamics (BOMD) in the calculation of time-averaged NMR chemical shifts. Our case study is the very flexible glycerol molecule that undergoes interconversion between conformers in a nonrandom way. Along the sequence of structures from one backbone conformer to another, transition states have been identified. The three ¹³C NMR chemical shifts of the molecule were estimated by averaging their calculated values over a large set of BOMD snapshots. The simulation time needed to obtain a good agreement with the two signals present in the experimental spectrum is shown to be dependent on the atomic orbital basis set used for the dynamics, with a necessary longer trajectory for the most extended basis sets. The large structural deformations with respect to the optimized conformer geometries that occur along the dynamics are related to a kinetically driven conformer distribution. Calculated conformer type populations are in good agreement with experimental gas phase microwave results.