4
I.D. Gorokh et al. / Journal of Molecular Structure 1199 (2020) 126955
Table 3
Values of the density of all electrons (r(r)), Laplacian of electron density (V [2]r(r)) and appropriate l2 eigenvalue, energy density (Hb), potential energy density (V(r)), and
Lagrangian kinetic energy (G(r)) (a.u.) at the bond critical point (3, ꢀ1), corresponding to non-covalent interactions Br/Br (type I halogen-halogen short contacts) in su-
pramolecular structure of 1, bond length (l) (Å), Wiberg bond index (WI), as well as energy for this contact Eint (kcal/mol), defined by different approaches.
a
b
c
d
r
(r)
V [2]
0.018
r
(r)
l2
Hb
V(r)
G(r)
Eint
1.3
Eint
1.1
Eint
1.5
Eint
1.4
le
WI
0.007
ꢀ0.008
0.000
ꢀ0.004
0.004
3.647
0.01
a
Eint ¼ eV(r)/253
.
b
Eint ¼ 0.429G(r) [54].
c
d
e
Eint ¼ 0.58(ꢀV(r)) (correlation developed exclusively for non-covalent interactions involving bromine atoms) [55].
Eint ¼ 0.57G(r) (correlation developed exclusively for non-covalent interactions involving bromine atoms) [54].
The shortest van der Waals radius for Br atom is 1.83 Å29
.
contacts are 2.710 Å (with H atoms of phenyl rings), and the Br/Br
distances exceed the sum of corresponding van der Waals radii,
indicating lack of such interactions. Both complexes were isolated
as pure phases, according to the PXRD data (see SI).
contacts [41,42]) takes place as well.
For estimation of the energies of these type I halogen-halogen
short contacts which, in our opinion, are related to the phenome-
non of halogen bonding [41,42], we performed DFTcalculations and
topological analysis of the electron density distribution (QTAIM
theory) [43] for 1. Previously, this approach was successfully uti-
lized for the studies of various non-covalent interactions in coor-
dination compounds [44e51]. Results are summarized in Table 3,
For estimation of relative contributions of different contacts to
the crystal packing, we performed Hirshfeld surface analysis (HSA)
for 1 and 2, see Table 2 and SI for details. For the visualization, we
have used a mapping of the normalized contact distance (dnorm); its
negative value enables identification of molecular regions of sub-
stantial importance for detection of short contacts. Fig. 4 depicts
the Hirshfeld surfaces for 1 and 2 (anionic and cationic parts). In
these Hirshfeld surfaces, the regions of shortest intermolecular
contacts visualized by red circle areas. In both cases, the domi-
nating contacts are Br/H and/or H/H (this situation is common
for halometalates [39,40]); however, in the case of 1, contribution of
abovementioned Br/Br interactions (type I halogen-halogen short
the contour line diagram of the Laplacian distribution V [2]r(r),
bond paths, and selected zero-flux surfaces as well as reduced
density gradient (RDG) isosurface [52] for non-covalent in-
teractions Br/Br in supramolecular structure of 1 are shown in
Fig. 5.
The Laplacian of electron density is typically decomposed into
the sum of contributions along the three principal axes of maximal
variation, giving the three eigen values (l1, l2 and l3) of the Hessian
matrix. In line with common practice, the sign of l2 can be utilized
to distinguish bonding (attractive, l2 < 0) weak interactions from
non-bonding ones (repulsive, l2 > 0). Thus, the discussed above
non-covalent interactions Br/Br (type I halogen-halogen short
contacts) in supramolecular structure of 1 are attractive (Table 3).
The QTAIM analysis demonstrates the presence of appropriate
bond critical point (3, ꢀ1) (BCP) for non-covalent interactions
Br/Br in supramolecular structure of 1 (Table 3). The low magni-
tude of the electron density, positive value of the Laplacian, and
very low energy density in this BCP are typical for non-covalent
interactions involving halogen atoms [56,57]. We have defined
energy for this weak contact Br/Br according to several proced-
ures [53e55]; depending on the approach used, the strength of this
contact vary in the ranges 1.1e1.5 kcal/mol. The balance between
the Lagrangian kinetic energy G(r) and potential energy density
V(r) at the BCP reveals that a covalent contribution is absent in
Br/Br supramolecular contacts in 1 [58], which correlates well
with negligible Wiberg bond indices for these interactions. The
NBO charges on the bromine atoms involving in non-covalent in-
teractions Br/Br in supramolecular structure of 1 are ꢀ0.44
(anionic part) and 0.12 (cationic part).
Overall, analysis of both geometries and energies of Br/Br in-
teractions reveals that those should be rather classified as Type I
halogen$$$halogen short contacts according to the classification by
Metrangolo et al.
According to the TGA data (see SI), both 1 and 2 reveal high
thermal stability (decomposition begins at z 200 ꢃC for
1
and z220 ꢃC for 2). These observations agree well with the in-
vestigations of thermal stability of bromobismuthates(III) made by
us earlier [59].
Fig. 5. Contour line diagram of the Laplacian distribution V [2]
selected zero-flux surfaces (top) and RDG isosurface (bottom) referring to non-
covalent interactions Br/Br in (packing induced contact aka type halogen-
halogen short contact). Bond critical points (3, ꢀ1) are shown in blue, nuclear crit-
ical points (3, ꢀ3) e in pale brown, length units e Å, isovalues on contour lines and
RDG isosurface values on the color scheme are given in a.u.
r(r), bond paths and
4. Conclusions
1
I
With 4-bromobenzyl-substituted pyridinium salts as pre-
cursors, we prepared two new bromobismuthate complexes.
Combination of Hirshfeld Surface Analysis (HSA) and DFT show that