3.1. X-ray Diffraction Study
Molecular structure of the gold(I) complex [AuCl(CNC
6H
3-4-Cl-2-I)] (
1) (
Figure 1) is analogous to those reported for other [AuCl(CNR)] complexes [
50,
51,
52,
53]. The gold(I) ion is in a linear coordination environment, which is formed by one isocyanide and one chlorido ligands. Bond distances and angles are typical for this type of complexes [
50,
51,
52].
In the structure of
1, each gold(I) center forms several intermolecular contacts. The shortest separation, which involve a metal atom, is Au···Au (3.3050(8) Å; vs Bondi [
54,
55] Σ
vdW 3.32 Å); this metallophilic bond is characteristic for (RNC)[Au
I] species [
20,
50,
52]. The Au···Au contact is involved in the system of noncovalent interactions which includes also two antiparallel C
CN···Cl short separations. The linear C–Au–Cl group is subject to the solid-state head-to-tail pairing accordingly to the association mode predicted and studied computationally [
20] for similar gold(I) species. Another Au···Au contact in the system of noncovalent interactions includes two C
CN···Au and one Au···Au contact (
Figure 1). In the latter aurophilic contact, the Au···Au distance is 3.5672(7) Å. It is still larger than the sum of the corresponding Bondi Σ
vdWAu + Au = 3.32 Å, but lower than Alvarez [
54,
55] Σ
vdWAu + Au = 4.64 Å.
Other unconventional contacts in the crystal structure are C
CN···Cl (3.428(6) Å) and C
CN···Au (3.495(5) Å). For the C
CN···Cl separation, the C···Cl distance is smaller than Bondi Σ
vdWC + Au 3.45 Å. Geometry considerations based on the data collected in
Table 1 indicate that the Cl nucleophilic site is directed to the isocyano carbon. The availability of the π-hole
CN···LP
Cl interaction was confirmed in appropriate theoretical study (see later).
ϕ is a torsion angle ClAuAuC.
To find our similarities with the crystal structure of our gold complex – which based on one Au···Au and two C
CN···Cl contacts – we performed a Cambridge Structural Database (CSD) search using the following parameters: the structures of [AuX(CNR)] (X = Cl, Br, I) exhibiting one Au···Au contact with
d(Au + Au)<Bondi Σ
vdW(Au + Au) and two antiparallel C
CN···X with
d(C + X)<Bondi Σ
vdW(C + X). By applying these criteria, we retrieved six structures of [AuCl(CNAr)] species, but only 4 out of 6 structures were of acceptable quality (
Rw <5%). Examination of the geometrical parameters of the four structures (
Table 1) allows the assignment of these contacts to aurophilic Au···Au and π-hole
CN···LP
Cl interactions.
Another unconventional contact is C
CN···Au, in which the distances between the C-atom of the isocyanide group and the gold(I) site are slightly larger than Bondi Σ
vdWC + Au 3.36 Å, but are significantly less than Alvarez Σ
vdWC + Au = 4.09 Å [
56]. The C
CN···Au contact exhibits a comparable length with those in the earlier reported structures of [AuCl(CNR)] (R = C
6H
4OMe-4; BUVCAX [
57]; 3.662(5) Å and R = 2-naphthyl; TAHVIK;[
58] 3.605(8) Å) and even in other (RNC)Au
I species (see CSD search with brief analysis for C
CN···M contacts in refs. [
27,
59]). In these reports, π-hole···[Au
I] nature of the C···Au contact has not been analyzed, while analogous π-hole···[M
II] (M = Pd, Pt) contacts were thoroughly investigated, including DFT studies of their nature [
27,
59].
Based on the angular parameters analysis (both ∠N≡C···Au and ∠C···Au–C angles are
ca. 100°) one can conclude that the C and Au atoms are in contact but parallel displaced. Our theoretical calculations revealed no bond critical point and bond path and only a small value of energetic stabilization due to electron transfer from Au(5
dx2-y2) to the π*(C≡N) orbital. As further explained below, the NCIplot analysis also suggests the presence of a weak interaction and confirms that both the CN group and Au-atoms are in contact. The existence of such π-hole
CN···[Au
I] interaction seems quite logical as gold(I) site, despite its positive charge, exhibits
d-orbital nucleophilicity and can function as a nucleophilic coformer of noncovalent interactions, for instance, halogen bonding [
60,
61,
62]. However, the results of our geometry considerations and the computations data give only a hint that such interaction likely exists and it is very weak.
The system of Au···Au and C
CN···Cl contacts accomplishes 2D extended ladder-type architecture (
Figure 2). Molecules from different 2D-ladder arrays are linked to each other via I···I contacts (3.7932(7) Å vs. Bondi Σ
vdWI + I = 3.96 Å) (
Figure 3). These contacts, as follows from consideration of the ∠C–I–I values (150.10(15) and 128.76(14)°), are close to Type-II halogen···halogen interactions (idealized angles
θ1 ≈ 180 and
θ2 ≈ 90°) [
63,
64]. In addition, molecules from the same 2D-ladder array exhibit short I···I contacts (4.0486(1) Å vs Bondi Σ
vdWI + I = 3.96 Å) of Type-I (∠C–I–I values are 103.04(14) and 76.96(14)° that are close to 90°) halogen···halogen interactions [
64]. To be conclusive regarding the nature of such contacts, an elaborated theoretical study is detailed below, where different computational tools were used to investigate the donor/acceptor role of the I-atoms.
3.2. Theoretical Study
First, to study the electron rich and electron poor parts of compound
1, the MEP surface was computed and corresponding surface is represented in
Figure 4. It can be observed that the MEP minimum, as expected, is located at the chlorido ligand (–33.8 kcal/mol). The MEP is also negative at the Au and negative belts around the iodine and chlorine substituents (–5.0 kcal/mol). The MEP maximum is located at the iodine’s σ-hole (33.9 kcal/mol) and it is also large and positive at the aromatic H-atoms, ranging from 30 to 34 kcal/mol. The MEP is also positive over the center of the aromatic ring (16.9 kcal/mol), σ-hole of chlorine (18.8 kcal/mol), and over the isocyano group (11.9 over C and 18.5 over N, measured perpendicular to the molecular plane). This distribution of the electron density explains the formation of the 2D-ladder type architecture with the antiparallel arrangement of the molecules.
To characterize the 2D-ladder packing and to gain an insight into the nature of the noncovalent interactions, DFT calculations were performed using two independent dimers, denoted as “A” and “B” that exemplify the two binding modes observed in the infinite assembly depicted in
Figure 2. The energetic features of the dimers and the QTAIM/NCIplot analyses are given in
Figure 5. As can be inferred from consideration of our computational results, both binding modes exhibit very similar interaction energies (–13.6 kcal/mol and –13.4 kcal/mol for dimers A and B, respectively). Both dimers A and B are interconnected by five bond critical points (CPs, represented as small red spheres) and bond paths (orange lines). One bond CP interconnects both Au-atoms, thus disclosing the existence of aurophilic interactions. Moreover, in both dimers the Cl is connected to the isocyano group by a bond CP and bond path, in the case of dimer A to the C-atom and in the case of B to the N-atom. Finally, in dimer A, two symmetrically equivalent CPs and bond paths connect two aromatic CH bonds to the chlorido ligands. For dimer B, the additional bond CPs and bond paths connect the Cl to one C-atom of the aromatic ring, thus evidencing the occurrence of π···LP interactions.
Such combination of interactions explains the rather strong dimerization energies. The main difference between both dimers is observed in the NCIplot analysis and shape of the reduced density gradient (RDG) isosurfaces. That is, dimer A exhibits well defined disk-shape isosurfaces for each contact, revealing that the Au···Au interaction is the strongest one (blue RDG isosurface). In contrast, dimer B presents a continuous and green RDG isosurface, thus suggesting a more complicated (cooperative) binding and supporting the possibility of a partial contribution of π-holeCN···[AuI] interaction in combination with the aurophilic one and also π-holeCN···Cl contacts.
To shed light into the physical nature interactions in both dimers, we performed the natural bond orbital (NBO) analysis since it is useful to study charge-transfer effects from an orbital donor-acceptor viewpoint. The most important contributions (>1 kcal/mol) obtained from the NBO analysis of both dimers are provided in
Figure 6. For dimer A, two orbital donor-acceptor contributions were found: one corresponds to the aurophilic interaction that is composed by an electron donation from the LP located at the filled 5
dx2-y2 atomic orbital of gold to the antibonding C–Au orbital of the adjacent molecule, and vice versa. The associated stabilization energy is E
(2) = 3.7 kcal/mol due to 5
dx2-y2[Au]→σ*(C–Au) is the largest donor-acceptor orbital contribution. In addition, an electron donation from one LP of the chlorido ligand to one of the π*(C≡N) orbitals of the isocyano group is also observed, with a concomitant second order stabilization of E
(2) = 1.4 kcal/mol.
For dimer B it is interesting to highlight that the most important charge transfer is from the LP of gold (5dx2-y2) to the π*(C≡N) orbital, thus confirming the existence of the interaction between the π-hole at the isocyano group and gold.
Echeverria[
20] has also analyzed C
CN···Cl contacts in Au
I dimers similar to those observed in Dimer A (characterized by bond CPs and disk shaped RDG isosurfaces). The reported NBO stabilization energies are ranged from 0.05 to 0.2 kcal/mol for the LP(Cl)→π*(C≡N) donor-acceptor interactions, thus revealing that the C
CN···Cl contacts studied herein are stronger than those previously studied.
Finally, the noncovalent three-center I···I contacts in the structure of
1 (
Figure 3) were also studied using QTAIM analysis and the electron localization function (ELF) method to verify halogen bonding (for reviews on halogen bonding see refs. [
1,
2,
3,
4,
5,
65]) nature of these contacts. The results are shown in
Figure 7 and it can be observed that the QTAIM analysis reveals the existence of three I···I contacts characterized by the corresponding bond CPs and bond paths. The trimerization energy is moderately strong (–10.3 kcal/mol) likely due to the contribution of the π-stacking of two out of the three molecules that form the I···I triangle. The ELF 2D-map is represented in
Figure 7b, using the plane defined by the three I-atoms for the orientation. This ELF analysis reveals that the I···I contacts characterized by the bond CPs denoted as “a” and “b” can be defined as halogen bonds since the bond path that connects the I-atoms crosses LP of one I (electron donor) and the σ-hole of the other one (electron acceptor). In contrast, the bond CP denoted as “c” corresponds to a I···I interaction that cannot classified as halogen bond, since the bond path crosses the LPs of both iodine atoms that are, in fact, pointing to each other. The size of these electron rich regions is significantly smaller than other similar regions (LPs) that are not pointing to each other, indicating some polarization and/or conjugation of the LPs to the aromatic ring, likely reducing the electrostatic repulsion. This I···I van der Waals contact can be rationalized considering that other forces like dispersion or polarization compensate the electrostatics or that is a consequence of packing effects.