Towards the Accurate Calculation of 183W NMR Chemical Shifts in Polyoxometalates: The Relevance of the Structure

DFT calculations were carried out to study ¹⁸³W NMR chemical shifts in the family of the Keggin anions with formula α‐XW₁₂O₄₀]^q? (X=B, Al, Si, P, Ga, Ge, As, Zn), in the β‐ and γ‐SiW₁₂O₄₀]^4? geometric isomers, in the derivative Dawson anion P₂W₁₈O₆₂]^6?, and in the most symmetrical Lindqvist W₆O₁₉]^2? anion and its derivative W₁₀O₃₂]^4?. In this article, we show that the geometry employed in the calculation of NMR chemical shifts in polyoxotungstates is extremely important if we want to be quantitative. Using very large basis sets of QZ4P quality and taking into account the conductor‐like screening model (COSMO) to account for solvent effects (aqueous and organic solutions), good geometries were found for the polyoxoanions. From these optimal geometries the ¹⁸³W NMR chemical shifts were computed with the more standard basis sets of TZP quality and including spin–orbit corrections inside the zero‐order regular approximation (ZORA) to describe the relativistic effects of the internal electrons. With this strategy the mean absolute error between experimental and theoretical values was found to be less than 10 ppm, which is similar to the experimental error. We also discuss how the geometry of the polyoxoanion influences on the shielding.