J. Autschbach et al.
potential is included in V’ are treated relativistically. The approach in es-
sence amounts to a position-dependent scaling of the speed of light in
order to suppress relativistic effects locally, with the atomic potentials in
k providing a suitable and convenient set of switching functions. Fur-
ther, one may use a partial V’ only in the SO term, along with the full V’
in the scalar ZORA operator. This then creates a relativistic operator in
which scalar relativistic effects are fully present, but SO coupling is sup-
pressed for selected atoms. It needs to be stressed that the suppressed
“atom in molecule” contributions are specific to the molecule under con-
sideration and not predetermined for a given atom type. A similar idea
has previously been applied by van Wꢀllen in ZORA calculations of the
magnetic anisotropy of a molecular magnet.[28]
Results and Discussion
Experimental structures: All molecules [X-Si-(m-mt)4-M-Y]
were found to reflect the shape of a fourfold propeller with
the methimazolyl moieties tilted in the same direction. Both
Si and M are located within a square planar environment of
four N or four S atoms, respectively, almost perfectly in
plane. Flexibility of the molecular architecture is reflected
by the variable bond lengths along the XSiMY axis and by
the variable tilt angles of opposing methimazolyl groups for
one compound. These features are summarized in Table 1.
All computations were carried out with a locally modified 2011 developer
version of ADF.[29] A new keyword, partialzora, reads a list of atoms or
atom types for which the atomic potentials are subsequently excluded in
the construction of k. Since the focus of this work is on SO effects on
NMR shielding, the suppression of relativistic effects has been imple-
mented separately for the scalar and the SO ZORA operator for the
ground-state DFT calculations, as outlined in the previous paragraph. In
the computations discussed below, the partial ZORA functionality has
been applied only to suppress SO effects, while keeping scalar relativistic
effects for all atoms. Modifications have been made accordingly in the
NMR nuclear magnetic shielding module in ADF. Here, the partial
ZORA functionality can be chosen separately for magnetic perturbation
operators originating from the scalar and the SO terms of the ZORA
Hamiltonian, respectively, for consistency with the separation of scalar
and SO effects in the ground state calculations.
ꢀ
With respect to M Si bond lengths, one can find an over-
all trend of decreasing interatomic distance when replacing
Ni by Pd by Pt. This is rationalized by the different orbital
contributions to the metal–silicon bond and the resultant in-
crease in covalent bond character when going to Pt (vide
infra) as well as typical variations in the metal atomic radii.
ꢀ
This M Si bond shortening has direct influence on the
ꢀ
trans-disposed Si X bond, which responds with bond length-
ening, to pronounced extent in the case of M=Pt.
As to the influence of the halide, replacing the Si-bound
ꢀ
Cl by F generally results in M Si bond lengthening. We can
qualitatively attribute this behavior to the thermodynamical-
The ADF package supports scaled-ZORA calculations in which the orbi-
tal energies and the total energy are scaled to improve them towards
fully relativistic results.[21] A similar concept has been suggested for NMR
calculations.[30] The partial ZORA implementation is presently incompat-
ible with scaled ZORA. Data reported in the Results and Discussion for
partial-ZORA shielding contributions are therefore based on unscaled
ZORA calculations. For consistency with our previous work on metallasi-
latrane complexes, all other 29Si shielding data reported herein were ob-
tained with scaled ZORA. For the 29Si shielding, the effects from ZORA
scaling are small (below 2 ppm) and almost completely cancel in the
chemical shifts.
Further details of the computations are as follows: In the 29Si NMR
shielding calculations, a combination of the Vosko–Wilk–Nusair (VWN)
local density functional[31] with revised Perdew–Burke–Ernzerhof[32]
(revPBE) generalized gradient approximation (GGA) terms was utilized
for compatibility with previous work. A hybrid variant of PBE with 25%
Hartree–Fock exchange, PBE0,[33] was used in order to investigate the
sensitivity of the results to approximations in the functional. To better
describe the metal (M) and silicon valence shells in the complexes, a flex-
ible and highly polarized Slater-type orbital (STO) basis set of quadru-
ple-z quality (QZ4P) was used for M and Si. For the ligand atoms in the
complexes, a triple-z doubly polarized STO basis set (TZ2P) designed for
ZORA computations was utilized. NMR chemical shifts were calculated
with the ZORA NMR shielding code of the ADF package.[30,34] All cal-
culations were performed with the X-ray diffraction (XRD) structures, in
selected cases after optimization of the hydrogen positions (HOP).
ꢀ
ꢀ
ly weaker Si Cl bond (compared to Si F), which can re-
spond more easily to trans-disposed lone-pair donation by
ꢀ
bond lengthening, thus strengthening the trans-Si M bond.
ꢀ
Similar features were observed for N!Si X moieties (X=
Cl, F), that is, a shorter N!Si bond for X=Cl.[36] A similar
trend, although to marginal extent, is found for the metal-
bound halide when Cl is replaced by Br by I.
ꢀ
In some crystal structures with several independent M Y
bonds one can see notable differences in M Y bond lengths,
in particular for entries 4, 6, 11, 12, and 16 in Table 1, which
reflect weak M Y bonding and great flexibility of these
ꢀ
ꢀ
bonds resulting therefrom.
Finally, the tilt angle of opposing methimazolyl ring
planes appears to be also very flexible. Although one can
discern the influence of the Si-bound halide (pronounced
tilt angle for X=Cl, resulting from intramolecular H-X-re-
pulsion, Figure 3), for the com-
pounds with the smaller F atom
at silicon crystal packing effects
seem to determine this tilt
angle to a pronounced extent.
This is particularly evident from
In order to investigate the impact of a polarized environment on the cal-
culated NMR parameters, additional computations were performed with
the conductor-like screening model (COSMO) as implemented in
ADF.[35] Whereas the dielectric constants (e) of standard insulating mate-
rials are typically below 10, we chose to investigate the polarization
effect of water (e=78), assuming that this gives an estimation of the
upper limit of the crystal packing or long-range electrostatic interactions
in the solids, as well as the influence of aqueous solvent in the case of
liquid-state NMR. We found that the application of COSMO led to
mean variations of the silicon shielding of about ꢃ5 ppm depending on
halides and the metal centers. Results are reported in Table S7 of the
Supporting Information. The overall effect is not dramatic. The observed
trends for the 29Si NMR parameters discussed herein can therefore be
considered as “intrinsic” to the complexes and not dominantly caused by
crystal environments in the solids.
entry 2 in Table 1, with tilt
angles ranging between 12.7(6)
and 39.0(4)8. The space-filling
Figure 3. Propeller tilt of op-
posing methimazolyl moieties
enforced by repulsion between
the Si-bound halide and methi-
mazolyl-4-H atoms.
ꢀ
models in Figure 4 show the H
X contacts for [F-Si-(m-mt)4-Pt-
Cl] (in I4) and for [Cl-Si-(m-
mt)4-Pd-I] (in Cc).
Overall performance of the calculations: Figure 5 shows the
correlation between ZORA calculations of Si chemical
shifts and experiment for representative members of the [X-
Si-(m-mt)4-M-Y] complex family that are analyzed in detail
&
4
&
ꢂ 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim
Chem. Eur. J. 0000, 00, 0 – 0
ÝÝ
These are not the final page numbers!