ABSTRACT
2H solid-state NMR and atomistic molecular dynamics (MD) simulations are used to understand the disorder of guest solvent molecules in two cocrystal solvates of the pharmaceutical furosemide. Traditional approaches to interpreting the NMR data fail to provide a coherent model of molecular behavior and indeed give misleading kinetic data. In contrast, the direct prediction of the NMR properties from MD simulation trajectories allows the NMR data to be correctly interpreted in terms of combined jump-type and libration-type motions. Time-independent component analysis of the MD trajectories provides additional insights, particularly for motions that are invisible to NMR. This allows a coherent picture of the dynamics of molecules restricted in molecular-sized cavities to be determined.
ABSTRACT
Van Vleck's classic theory of the second moment of lineshapes in 1H nuclear magnetic resonance (NMR) is reworked in a form that allows the effect of rapid molecular motion on second moments to be calculated in a semi-analytical fashion. This is much more efficient than existing approaches and also extends previous analyses of (non-dynamic) dipolar networks in terms of site-specific root-sum-square dipolar couplings. The non-local nature of the second moment means that it can discriminate between overall motions that are difficult to discriminate using alternative approaches, such as measurements of NMR relaxation. The value of reviving second moment studies is illustrated on the plastic solids diamantane and triamantane. In the case of triamantane, straightforward measurements of 1H lineshapes on milligram samples show that the molecules in the higher temperature phase undergo multi-axis jumps, information that is not accessible either to diffraction studies or to alternative NMR approaches. The efficiency of the computational methods means that the second moments can be calculated using a readily extensible and open-source Python code.
