Chemical Shift Tensors of Cimetidine Form A Modeled with Density Functional Theory Calculations: Implications for NMR Crystallography

Abstract

The principal components of the ¹³C chemical shift tensors for the ten crystallographically distinct carbon atoms of the active pharmaceutical ingredient cimetidine Form A have been measured using the FIREMAT technique. Density functional theory (DFT) calculations of ¹³C and ¹⁵N magnetic shielding tensors are used to assign the ¹³C and ¹⁵N peaks. DFT calculations were performed on cimetidine and a training set of organic crystals using both plane-wave and cluster-based approaches. The former set of calculations allowed several structural refinement strategies to be employed, including calculations utilizing a dispersion-corrected force field that was parametrized using ¹³C and ¹⁵N magnetic shielding tensors. The latter set of calculations featured the use of resource-intensive hybrid-DFT methods for the calculation of magnetic shielding tensors. Calculations on structures refined using the new force-field correction result in improved values of ¹⁵N magnetic shielding tensors (as gauged by agreement with experimental chemical shift tensors), although little improvement is seen in the prediction of ¹³C shielding tensors. Calculations of ¹³C and ¹⁵N magnetic shielding tensors using hybrid functionals show better agreement with experimental values in comparison to those using GGA functionals, independent of the method of structural refinement; the shielding of carbon atoms bonded to nitrogen are especially improved using hybrid DFT methods.