Citation for published item:ohD ynd § rej nd rodgkinsonD ul nd iddi(eldD gory wF nd tesD tonthn F nd hr § ¡ %nsk¡ yD wrtin @PHIUA 9ixploring systemti disrepnies in hp lultions of hlorine nuler qudrupole ouplingsF9D tournl of physil hemistry eFD IPI @PIAF ppF RIHQERIIQF 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-prot 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. ABSTRACT: Previous studies revealed significant discrepancies between DFT-calculated and experimental nuclear quadrupolar coupling constants (C Q ) at chlorine atoms, particularly in ionic solids. Various aspects of the computations are systematically investigated here, including the choice of the DFT functional, basis set convergence, and geometry optimization protocol. The effect of fast (fs) time-scale dynamics are probed using molecular dynamics (MD) and nuclear quantum effects (NQEs) are considered using path-integral MD calculations. It is shown that the functional choice is the most important factor related to improving the accuracy of the 2 quadrupolar coupling calculations and functionals beyond the generalized gradient approximation (GGA) level, such as hybrid and meta-GGA functionals, are required for good correlations with experiment. The influence of molecular dynamics and NQEs is less important than the functional choice in the studied systems. A method which involves scaling the calculated quadrupolar coupling constant is proposed here; its application leads to good agreement with experimental data.