To cite this article: Jochen Autschbach (2013) The role of the exchange-correlation response kernel and scaling corrections in relativistic density functional nuclear magnetic shielding calculations with the zeroth-order regular approximation,The relativistic NMR module of the Amsterdam Density Functional (ADF) package, which is frequently utilised in studies of heavy atom NMR chemical shifts, is extended to calculate a hitherto neglected term from the response of the exchangecorrelation (XC) potential. The term vanishes in the absence of spin-orbit coupling. Further, corrections to the shielding arising from scaling factors in the zeroth-order regular approximation (ZORA) relativistic framework are investigated. The XC response markedly improves calculated proton chemical shifts for hydrogen halides. Mercury chemical shifts for mercury dihalides are also noticeably altered. Contributions from density-gradient dependent terms in the response kernel contribute about 30-40%. New fully relativistic density functional theory (DFT) benchmark data are compared with ZORA and literature reference values. In line with previous work, it is found that absolute shielding constants for Hg are not accurately predicted with ZORA. However, chemical shifts agree well with fully relativistic calculations. The application of 'scaled-ZORA' scaling factors deteriorates the shielding constants and is therefore not recommended. The scaling hardly affects chemical shifts. ZORA calculations are not suitable for absolute shielding of heavy atoms but they can be used safely for chemical shifts in most application scenarios.
IntroductionFirst-principles calculations of structure, thermochemistry, and spectroscopic properties of molecules with heavy elements require a relativistic quantum chemistry framework [1][2][3][4][5][6][7]. This is particularly true for molecular properties for which the relevant quantum mechanical operators predominantly sample the electronic structure close to a nucleus. NMR parameters are prime examples of such properties. The zeroth-order regular approximation (ZORA) [8-10] is a variationally stable approximate (quasirelativistic) twocomponent (2c) all-electron relativistic framework. In conjunction with density functional theory (DFT), ZORA is widely applied to calculations of NMR parameters. For a small selection of recent applications to a variety of problems in chemistry see, for instance [11][12][13][14][15][16]. A ZORA module for NMR shielding tensors has been part of the Amsterdam Density Functional (ADF) [17] package since 1999 [18,19], followed by an implementation of a module for calculations of indirect nuclear spin-spin coupling (J-coupling) tensors in the year 2000 [20,21]. Among the reasons for the early development of ZORA-based magnetic property modules have been the ease by which magnetic perturbation operators can be derived and the rather straightforward way of calculating matrix elements with the help of numerical integration techniques. Extensions for hybrid DFT calculations have been