Citation for published item:hr § ¡ %nsk¡ yD w rtin nd fourD etr nd rodgkinsonD ul @PHITA 9 emper ture dependen e of xw p r meters l ul ted from p th integr l mole ul r dyn mi s simul tionsF9D tourn l of hemi l theory nd omput tionFD IP @QAF ppF WTVEWUQF Further information on publisher's website: Use policyThe full-text may be used and/or reproduced, and given to third parties in any format or medium, without prior permission or charge, for personal research or study, educational, or not-for-pro t purposes provided that:• a full bibliographic reference is made to the original source • a link is made to the metadata record in DRO • the full-text is not changed in any way The full-text must not be sold in any format or medium without the formal permission of the copyright holders.Please consult the full DRO policy for further details. The influence of temperature on NMR chemical shifts and quadrupolar couplings in model molecular organic solids is explored using path integral molecular dynamics (PIMD) and density functional theory (DFT) calculations of shielding and electric field gradient (EFG) tensors. An approach based on convoluting calculated shielding or EFG tensor components with probability distributions of selected bond distances and valence angles obtained from DFT-PIMD simulations at several temperatures is used to calculate the temperature effects. The probability distributions obtained from the quantum PIMD simulations, which includes nuclear quantum effects, are significantly broader and less temperature dependent than those obtained with 2 conventional DFT molecular dynamics or with 1D scans through potential energy surface.Predicted NMR observables for the model systems were in excellent agreement with experiment.