2.2. DFT Calculations
In order to elucidate the molecular structures of fragment ions detected in mass spectrometry experiments, a series of DFT calculations has been performed, starting from the computation of the minimum energy structures of the three neutral complexes [Ni(L)
2TMEDA] (
Figure 3). Additionally, all fragments discussed in the previous section have been subjected to structural optimizations. The resulting optimized geometries are depicted in
Figure 3,
Figure 4,
Figures S1 and S2, while relevant geometrical parameters are reported in
Table 2 2,
Table 4 and
Table S1.
A graphical representation of the optimized structures of the neutral complexes obtained for L= tfa, fod, and thd is shown in
Figure 3a,c,e, respectively. Since all ESI-HRMS spectra showed a base peak due to the loss of a β-diketonate ligand L, a search for the minimum energy structures of [M-L]
+ species has been performed starting from the optimized structures of neutral complexes and removing a β-diketonate ligand. The resulting optimized geometries of [Ni(L)TMEDA]
+ with L = tfa, fod, and thd are depicted in
Figure 3b,d,f, respectively.
As far as the neutral complexes are concerned, their molecular geometries are similar among each other and close to the structural parameters computed with another density functional approximation [
32]. In particular, in all the three neutral complexes Ni exhibits a slightly distorted octahedral coordination environment, characterized by Ni-O distances shorter than Ni-N ones (see
Table 2), suggesting that the interactions between metal center and oxygen atoms should be stronger than those with nitrogen ones. This hypothesis is substantiated by the bond order (BO) trend, as the Ni-O BOs are systematically higher than the Ni-N ones for all three compounds (see
Table 2). Interestingly, whereas BO values for Ni-O bonds are quite similar among the three complexes, the Ni-N bond order for the non-fluorinated compound
3 - bearing two tert-butyl group
per diketonate ligand - is slightly lower than that found for fluorinated complexes
1 and
2. These data indicate that, whereas F-containing groups cope to increase the Ni-N bond strength, the presence of tert-butyl substituents has the opposite effect –
i.e., weakening Ni-N interactions. Nonetheless, beside these slight differences, both the molecular geometry and the electronic structure underlying the metal-ligand bonding scheme show a strong similarity in the three compounds, as evidenced by
Figure 3 and the data in
Table 2. This result might rationalize the similar fragmentation behavior of the target complexes as deduced from
Figure 1, in particular the presence of [M-L]
+ as the base peak in all ESI-HRMS spectra.
Indeed, similar considerations hold for the minimum energy structures calculated for [M-L]
+ fragments, all characterized by very similar geometric arrangements and metal-ligand bonding schemes. Specifically, as illustrated by
Figure 3b,d,f, all three ions are characterized by a square-planar geometry, typical of tetra-coordinated Ni(II) complexes in the singlet spin state. Moreover, data in
Table 2 indicate that both geometrical parameters and BOs values show only very slight differences in passing from compound
1 to
2 to
3. This finding is further supported by the Natural Bond Orbital (NBO) analysis performed on the electronic structure of the three [M-L]
+ fragments. In particular, the data reported in
Table S2 indicate a very strong similarity of the NBOs localized on the metal-ligand bonds. Irrespective of the diketonate nature, all NBOs have σ-character. While the bonding components (BD(σ)) are mostly localized on ligand’s nitrogen and oxygen atoms (≈90%), the antibonding components (BD(σ*)) are predominantly localized on Ni. Taken together, NBO data indicate that electron charge is transferred by the ligands towards Ni. This finding is further clarified by the graphical representation of the aforementioned bonding NBO’s - depicted in
Figures S3, S4, and S5 for L = tfa, fod, and thd, respectively –, highlighting the net σ-character of the orbital, as well as Ni participation to the bonding scheme. Such a close similarity in the electronic structure might suggest that all three ions are characterized by a similar stability and may all easily form in ESI-HRMS conditions.
Notably, for all the target complexes, Ni-O and Ni-N distances undergo a significant shortening upon passing from neutral compounds to [M-L]+ fragments. This trend is paralleled by the considerable increase of the corresponding bond orders, which are nearly doubled in the fragment ions. Hence, in terms of the metal-ligand bonding, the loss of a diketonate ligand results in a significant strengthening of Ni-TMEDA interactions and of Ni-O bonds of the remaining diketonate moiety.
A useful feature of NBO analysis is the opportunity of partitioning the electronic density between ‘subsystems’ –
i.e., specific portions of a given molecule, allowing thus to estimate the total charge on each subsystem. In this case, it is particularly instructive to inspect the total NBO charges localized on Ni, diketonate ligand L, and TMEDA, reported in
Table 3. Interestingly, the data show that the positive charge on Ni decreases by ≈ 25-30% on passing from neutral complex to [M-L]
+ cations. This result, that may appear counterintuitive given the fragments’ positive charge, may be explained by considering that the donation of electronic charge from both L and TMEDA to Ni increases significantly upon the loss of a diketonate ligand. Indeed, the total negative charge on L decreases by 50%, whereas the total (positive) charge on TMEDA is more than doubled (
Table 3). Hence, these data clearly indicate that all [M-L]
+ species are stabilized by a strong donation of electronic density from the ligands towards the metal center. This effect is slightly stronger in the fragment derived from compound
3, due to the higher electron-donor character of the two tert-butyl substituents on the diketonate moiety with respect to fluorine-containing groups, present in compounds
1-2.
Beside the main peak, which is common to the three complexes, ESI-HRMS spectra indicated the formation of other cations (see
Figure 1 and
Table 1). A careful computational search of the minimum energy structure of these species has been performed as well. Graphical representations of optimized geometries of [M+Na]
+, [M-TMEDA+Na]
+, [M-TMEDA+H]
+, [HL+H]
+ and [TMEDA+H]
+ are shown in
Figure S1, while selected geometrical parameters are reported in
Table S1. It should be pointed out that Na
+ is present only in ESI-HRMS spectra of the fluorine-richest compound
2. In [M+Na]
+, Ni maintains its octahedral environment in spite of the entrance of Na
+ in the second coordination shell (
Figure S1a). Specifically, three O atoms appear to be at coordination distances from Na (average Na-O distance = 2.416 Å), which is also interacting with a fluorine atom (Na-F distance = 2.441 Å). The consequent distortion of the octahedral Ni environment is evidenced by the lengthening of Ni-O and Ni-N bonds with respect to the neutral complex [Ni(fod)
2TMEDA]. This effect is particularly relevant for the bonds involving the three oxygen atoms also coordinated with Na
+ (
Table S1). Similarly, the species [M-TMEDA+Na]
+ - exhibiting a distorted square-planar Ni coordination geometry, is characterized by significant O-Na
+ and F-Na+ interactions (
Table S1,
Figure S1b). In a different way, the [M-TMEDA+H]
+ species - namely [Ni(thd)
2+H]
+ - is protonated at the CH1 carbon. For this ion, two structures very close in energy (ΔE = 0.32 kcal/mol) have been found (
Figure S1c,d). The most stable one exhibits a distorted tetrahedral arrangement around the Ni center, and has triplet spin multiplicity, whereas the less stable one shows a distorted square planar arrangement and singlet spin multiplicity. Given the modest energy difference, both structures might be formed during the ESI process. The presence of the protonated diketonate species [HL+H]
+ has been detected only for the case of the more electron-donor thd ligand. Its optimized structure is characterized by the protonation of both carbonylic groups (
Figure S1e). On the other hand, the other protonated ligand [TMEDA+H]
+ is detected in the fragmentation of all three complexes. In this case, the protonation occurs on one of the nitrogen atoms (see
Figure S1f).
As far as MS
2 spectra are concerned, attention was focused on the subsequent fragmentation of [M-L]
+ cations. Even in this case, the spectra of the three cations presented common features, in particular the formation of [M-L]
+ - NH(CH
3)
2 species. The theoretical search of the minimum energy structure of these moieties led to the molecular geometries depicted in
Figure 4a,c,f for L = tfa, fod, thd, respectively. Remarkably, these cations have very similar geometries and Ni-ligands bonding patterns. Whereas the diketonate remains coordinated to Ni with both its oxygen atoms, only one nitrogen atom (the one survived to the fragmentation) is linked to the metal center. Nevertheless, Ni maintains its nearly square-planar coordination geometry thanks to the interaction with a double C=C bond formed upon - NH(CH
3)
2 release. Specifically, the Ni-CHT and Ni-CH2T distances become ≈ 2 Å in these fragments (see
Table 4). Moreover, data in
Table 4 clearly indicate a net shortening of CHT-CH2T bond length (
Figure 4) compared to the same distance in [M-L]
+ species (
Table 2). This distance shortening is accompanied by a significant increase of the CHT-CH2T bond order (from values of ≈ 1 to values greater than 1.5, see
Table 2 and
Table 4). Correspondingly, Ni-CHT and Ni-CH2T BOs values become comparable to those of Ni-O/Ni-N bonds, indicating appreciable Ni-C interactions which are greater for the terminal C atom (CH2T), closer to Ni. These observations are common to the three cations, underlining once again the similarity of the metal-ligand bonding scheme among systems derived from compounds
1,
2,
3.
The loss of a -CH
3CH
2N(CH
3)
2 group is detected only for the fluorinated [M-L]
+ species derived from complexes
1 and
2. In these cations (see
Figure 4 b,d), the bonding pattern around Ni is practically the same as in [M-L]
+ - NH(CH
3)
2 species, discussed previously. The planar tetra-coordinated Ni structure is characterized by two Ni-O bonds, one Ni-N bond and a cation-π intramolecular interaction between the Ni center and the CHT=CH2T double bond.
Figures S6-S10 report a graphical representation of the relevant orbitals involved in this cation-π interaction, whereas
Tables S3 and S4 report the quantitative NBO analysis. As can be seen from both orbital representations and NBO analysis, Ni cation actively participates to such an interaction, with the π structure predominantly localized on the CHT=CH2T double bond.
Concerning complex
2, species corresponding to the loss of a CH
3C=CH
2 group from [Ni(fod)TMEDA]
+ have been detected ([Ni(fod)TMEDA]
+-(CH
3)
2C=CH
2) in ESI-MS spectra. Their stoichiometry indicates a rearrangement of the fod ligand (see
Figure 1 and
Table 1). In order to identify the structure of these cations, geometry optimizations have been performed on different hypothetical guess structures characterized by such a stoichiometry.
Figure S2 contains a graphical representation of the six most stable structures, along with their relative stability (including the zero-point contribution). The most stable structure (
Figure S2a) is characterized by a distribution of fluorine atoms on both diketonate ligand sides. A similar feature is also present in the other lower energy isomers.
Finally, the trend of the Bader charges on passing form the neutral complex to the positive ions further evidences common features in the fragmentation of the three complexes (see
Tables S5-S8). Indeed, it is possible to deduce a displacement of electrons from ligands to the metal center, whose positive charge is significantly depleted along the fragmentation. This finding is also in line with the arguments deduced from NBO analysis. Notably, such an electron-enrichment tendency of the metal center during [Ni(L)
2TMEDA] precursors’ fragmentation might be relevant for CVD processes, where an electron-richer Ni cation (nearly +1) could be more prone to be oxidized by the molecular oxygen present in the reaction atmosphere, facilitating thus the formation of NiO films.
Table 4.
Selected bond lengths (BL) (in Å) and corresponding bond orders (BO) from the optimized structures of fragments detected in ESI-MS
2 spectra of [M-L]
+ ions for compounds
1,
2, and
3. Different colors in the BO columns refer to different ligands (brown = diketonate, gray = TMEDA). Atom labels as in
Figure 4.
Table 4.
Selected bond lengths (BL) (in Å) and corresponding bond orders (BO) from the optimized structures of fragments detected in ESI-MS
2 spectra of [M-L]
+ ions for compounds
1,
2, and
3. Different colors in the BO columns refer to different ligands (brown = diketonate, gray = TMEDA). Atom labels as in
Figure 4.
Bond distance |
1 |
2 |
|
33 |
|
[M-tfa]+- NH(CH3)2
|
[M-tfa]+- CH3CH2N(CH3)2
|
[M-fod]+- NH(CH3)2
|
[M-fod]+- CH3CH2N(CH3)2
|
[M-fod]+- -(CH3)2C=CH2
|
[M-thd]+- NH(CH3)2
|
|
BL BO
|
BL BO
|
BL BO
|
BL BO
|
BL BO
|
BL BO
|
Ni-O1 |
1.811 0.523
|
1.799 0.557
|
1.814 0.506
|
1.794 0.554
|
1.837 0.450
|
1.808 0.535
|
Ni-O2 |
1.830 0.511
|
1.819 0.541
|
1.816 0.518
|
1.815 0.543
|
1.842 0.453
|
1.803 0.546
|
Ni-N1 |
- |
- |
1.929 0.391
|
1.946 0.402
|
1.938 0.415
|
1.941 0.368
|
Ni-N2 |
1.932 0.393
|
1.943 0.408
|
2.153 0.201 |
1.947 0.387 |
1.937 0.414
|
1.957 0.368 |
Ni-CHT |
1.972 0.248
|
1.987 0.264
|
1.976 0.239
|
1.986 0.264
|
2.165 0.178 |
1.974 0.243
|
Ni-CH2T |
2.096 0.394
|
2.128 0.368
|
2.098 0.393
|
2.128 0.363
|
2.165 0.178 |
2.108 0.375
|
O1-C5 |
1.277 1.288
|
1.281 1.267
|
1.280 1.258
|
1.281 1.263
|
1.265 1.342
|
1.288 1.264
|
O2-C7 |
1.281 1.297
|
1.279 1.308
|
1.279 1.314
|
1.280 1.299
|
1.270 1.318
|
1.287 1.264
|
C5-CH1 |
1.384 1.444
|
1.381 1.463
|
1.380 1.468
|
1.381 1.458
|
1.400 1.345
|
1.399 1.381
|
C7-CH1 |
1.409 1.315
|
1.412 1.300
|
1.417 1.281
|
1.413 1.292
|
1.389 1.404
|
1.403 1.361
|
CHT-CH2T |
1.379 1.566
|
1.373 1.588
|
1.380 1.567
|
1.373 1.590
|
1.508 1.027
|
1.378 1.579
|