To what extent do bond length and angle govern the 13C and 1H NMR response to weak CH⋯O hydrogen bonds? A case study of caffeine and theophylline cocrystals

1H and 13C chemical shift-structure effects in anhydrous β-caffeine and four caffeine-diacid cocrystals probed by solid-state NMR experiments and DFT calculations

By using density functional theory (DFT) calculations, we refined the H atom positions in the structures of β-caffeine (C), α-oxalic acid (OA; (COOH)2), α-(COOH)2·2H2O, β-malonic acid (MA), β-glutaric acid (GA), and I-maleic acid (ME), along with their corresponding cocrystals of 2 : 1 (2C-OA, 2C-MA) or 1 : 1 (C-GA, C-ME) stoichiometry. The corresponding 13C/1H chemical shifts obtained by gauge including projector augmented wave (GIPAW) calculations agreed overall very well with results from magic-angle-spinning (MAS) nuclear magnetic resonance (NMR) spectroscopy experiments. Chemical-shift/structure trends of the precursors and cocrystals were examined, where good linear correlations resulted for all COO1H sites against the H⋯O and/or H⋯N H-bond distance, whereas a general correlation was neither found for the aliphatic/caffeine-stemming 1H sites nor any 13C chemical shift against either the intermolecular hydrogen- or tetrel-bond distance, except for the 13COOH sites of the 2C-OA, 2C-MA, and C-GA cocrystals, which are involved in a strong COOH⋯N bond with caffeine that is responsible for the main supramolecular stabilization of the cocrystal. We provide the first complete 13C NMR spectral assignment of the structurally disordered anhydrous β-caffeine polymorph. The results are discussed in relation to previous literature on the disordered α-caffeine polymorph and the ordered hydrated counterpart, along with recommendations for NMR experimentation that will secure sufficient 13C signal-resolution for reliable resonance/site assignments.

Quantifying the Intrinsic Strength of C-H⋯O Intermolecular Interactions

It has been recognized that the C-H⋯O structural motif can be present in destabilizing as well as highly stabilizing intermolecular environments. Thus, it should be of interest to describe the strength of the C-H⋯O hydrogen bond for constant structural factors so that this intrinsic strength can be quantified and compared to other types of interactions. This description is provided here for C_2h-symmetric dimers of acrylic acid by means of the calculations that employ the coupled-cluster theory with singles, doubles, and perturbative triples [CCSD(T)] together with an extrapolation to the complete basis set (CBS) limit. Dimers featuring the C-H⋯O and O-H⋯O hydrogens bonds are carefully investigated in a wide range of intermolecular separations by the CCSD(T)/CBS approach, and also by the symmetry-adapted perturbation theory (SAPT) method, which is based on the density-functional theory (DFT) treatment of monomers. While the nature of these two types of hydrogen bonding is very similar according to the SAPT-DFT/CBS calculations and on the basis of a comparison of the intermolecular potential curves, the intrinsic strength of the C-H⋯O interaction is found to be about a quarter of its O-H⋯O counterpart that is less than one might anticipate.