{'fraction_non_alphanumeric': 0.025322997416020673, 'fraction_numerical': 0.0031007751937984496, 'mean_word_length': 5.585034013605442, 'pattern_counts': {'":': 0, '<': 1, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Calculation of zero-field splitting parameters: comparison of a two-component noncolinear spin-density-functional method and a one-component perturbational approach.

Two different sets of approaches for the density-functional calculation of the spin-orbit contributions to zero-field splitting (ZFS) parameters of high-spin systems have been implemented within the same quantum chemistry code ReSpect and have been validated and compared for a series of model systems. The first approach includes spin-orbit coupling variationally in a two-component calculation, using either an all-electron Douglas-Kroll-Hess ansatz or two-component relativistic pseudopotentials. The ZFS parameters are computed directly from energy differences between different relativistic states. Additionally, an approximate second-order perturbation theory approach has been implemented, based on nonrelativistic or scalar relativistic wave functions. For a series of group 16 triplet diatomics and for the octet GdH3 molecules, two-component density functional calculations underestimate the zero-field splitting D systematically by a factor of 2. This may be rationalized readily by the incomplete description of states with absolute value MJ < J by a single-determinantal wave function built from two-component spinors. In the case of two 3d transition metal complexes and for GdH3, the results depend furthermore sensitively on exchange-correlation functional. Results of the alternative one-component approach agree strikingly with the two-component data for systems with small spin-orbit effects and start to deviate from them only for heavier systems with large spin-orbit effects. These results have fundamental implications for the achievable accuracy of one-component density-functional approaches used widely to compute ZFS parameters in the field of molecular magnetism. Possible refinements of both one-and two-component approaches are discussed.