Next Article in Journal
Stannylenes and Germylenes Stabilized by Tetradentate Bis(amidine) Ligands with a Rigid Naphthalene Backbone
Previous Article in Journal
Interactions between Damaged Hair Keratin and Juglone as a Possible Restoring Agent: A Vibrational and Scanning Electron Microscopy Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Analysis of Nickel–Phosphine Complexes with Cumulated Double Bond Ligands: Structural Insights and Electronic Interactions via ETS-NOCV and QTAIM Approaches

1
National Laboratory of Renewable Energy, University of Pécs, H-7624 Pécs, Hungary
2
Department of General and Inorganic Chemistry and HUN-REN Research Group for Selective Chemical Syntheses, Hungary, University of Pécs, H-7624 Pécs, Hungary
*
Author to whom correspondence should be addressed.
Molecules 2024, 29(2), 324; https://doi.org/10.3390/molecules29020324
Submission received: 7 December 2023 / Revised: 29 December 2023 / Accepted: 3 January 2024 / Published: 9 January 2024
(This article belongs to the Section Computational and Theoretical Chemistry)

Abstract

:
This study presents a comprehensive analysis of nickel–phosphine complexes, specifically Ni(PH3)2(OCCH2), Ni(PH3)2(H2CCO), Ni(PH3)2(H2CCCH2), Ni(PH3)2(NNCH2), and Ni(PH3)2( η 1 -H2CNN). Utilizing ETS-NOCV analysis, we explored orbital energy decomposition and the Hirshfeld charges of the ligands, providing insights into the electronic structures and donor–acceptor interactions within these complexes. The interactions in the ketene and allene complexes exhibit similar deformation densities and NOCV orbital shapes to those calculated for Ni(PH3)2(NNCH2), indicating consistent interaction characteristics across these complexes. The total interaction energy for all η 2 complexes is observed to be over 60 kcal/mol, slightly exceeding that of the analogous carbon dioxide complex reported earlier. Furthermore, the study highlights the stronger back-donation as compared to donor interactions across all η 2 complexes. This is further corroborated by Hirshfeld analysis, revealing the charge distribution dynamics within the ligand fragments. The research offers new perspectives on the electron distribution and interaction energies in nickel–phosphine complexes, contributing to a deeper understanding of their catalytic and reactive behaviors.

Graphical Abstract

1. Introduction

Reactions of carbon dioxide complexes with transition metals deserve special attention, as CO2 can be considered an inexhaustible raw material source and is used as a C1 building block for the production of various organic compounds (formic acid, methanol, and cyclic carbonates) [1,2,3]. Gaining knowledge about the coordination of carbon dioxide to metals (Figure 1) can provide valuable information, which can help develop more efficient catalytic systems [4].
Among the catalytic reactions based on carbon dioxide substrates, hydrogenation is noteworthy [5,6], which results in the formation of formic acid. The reaction was discovered by Farlow and Adkins in 1935, using Raney nickel as a catalyst [7]. The first application in homogeneous catalysis was published by Inoue and colleagues in 1976 [8]. With suitable catalysts, such as the Pd-dppm complex, the reduction can proceed to the formation of methane [9,10]. Other noteworthy homogeneous catalytic reactions based on CO2 activation include rhodium-catalyzed hydrosilylation [11,12,13], as well as coupling reactions with alkenes [14,15,16] and alkynes [17].
The diazo group of diazoalkanes can be formally exchanged for the isoelectronic carbon monoxide. The resulting ketenes are versatile, generally highly reactive intermediates in organic syntheses [18,19,20,21,22,23,24]. The carbonylation of diazomethane was first published by Staudinger and Kupfer [25]. For the production of ketene, carbon monoxide was bubbled through a 1% ether solution of diazomethane, and the resulting vapor/gas mixture was heated in a quartz tube above 400 °C (Equation (1)). The exiting gas, when reacted with an ether solution of aniline, yielded acetanilide, which indirectly proved the formation of ketene (Equation (2)).
CH 2 N 2 + CO N 2 Δ CH 2 = C = O
CH 2 = C = O PhNH 2 PhNHC ( = O ) CH 3
Ketenes are generally produced from the corresponding carboxylic acids or carboxylic acid derivatives. The industrial production of unsubstituted ketene is performed via pyrolysis of acetic acid [26]:
CH 3 COOH Δ H 2 C = C = O + H 2 O
Previously, we studied the carbonylation of diazoalkanes in the presence of various cobalt catalysts as well as the mechanism of the reaction [27,28,29]. Moreover, the mechanism of diazocarbonylation was scrutinized for the nickel–carbonyl- [30] and the iron–phosphine–carbonyl- [31] catalyzed model reactions as well.
In connection with our previous studies, the goal of the present study is to scrutinize the electronic structure of complexes Ni(PH3)2( η 2 -(O,C)-OCCH2) (1), Ni(PH3)2( η 2 -(C,C)-H2CCO) (2), Ni(PH3)2( η 2 -(C,C)H2CCCH2) (3), Ni(PH3)2( η 2 -(N,N)NNCH2) (4), and complex Ni(PH3)2( κ -C-H2CNN) (5), as well as that of the standalone ketene, allene, and diazomethane ligands.

2. Computational Details

All the structures were optimized without symmetry constraints with tight convergence criteria using the Gaussian suite of programs [32], with the exchange and correlation functionals developed by Perdew, Burke, and Ernzherhof [33] and denoted as PBEPBE. For all the atoms, the def2-TZVP basis set [34] was employed. Natural bond orbital (NBO) analyses have been performed using the GENNBO 7.0 program [35]. QTAIM (Quantum Theory of Atoms In Molecules) analyses of the wave function [36] were carried out with the AIMAll software (Version 19.10.12) [37]. For the EDA-NOCV (Energy Decomposition Analysis in combination with the Natural Orbitals of Chemical Valence methodology) calculations, the ADF part of the AMS 2023 software was used [38], employing the PBEPBE functional in combination with the triple- ζ STO basis set for all atoms with one set of polarization functions (denoted as TZP) and a small frozen core. The EDA-NOCV method, also known as ETS-NOCV (where ETS stands for extended transition state), is an attractive method for studying metal–ligand interactions in a wide variety of complexes [39,40,41,42,43,44,45,46,47,48]. For the Domain-averaged Fermi hole analyses, the WinFermi program was used [49].

3. Results and Discussion

As carbon dioxide can be considered as the ’prototype’ of the other small molecules with a cumulative double bond, some of its properties will be discussed first followed by those of the other standalone ligands. The electronic structures of nickel(0) complexes with carbon dioxide, carbonyl sulfide, and carbon disulfide were reported previously [50,51,52]. The Cartesian geometries as well as internal coordinates of all species occuring in this study are in the Supporting Information.

3.1. The Electron Structure of Carbon Dioxide

The CO2 molecule is generally depicted with cumulated double bonds, in which the carbon atom is in an s p hybrid state. In the ground state of carbon dioxide, the 1 π g orbital is the degenerate HOMO, which is an antisymmetric linear combination of the p orbitals centered on oxygen (Figure 2), while the LUMO ( 2 π u ) is more significantly located on the carbon atom. The frontier orbitals provide a fairly good qualitative description of the coordination of CO2 to various metals, where the oxygen’s lone pair participates as a donor ligand, while the back-donation interaction is directed towards the carbon atom. However, in stark contrast to the cumulated representation mode (where the symmetry planes of the π orbitals are perpendicular to each other), CO2 follows D h symmetry, i.e., it is cylindrically symmetrical [53].
The distribution of electron density can be visually represented by its second derivative with respect to spatial coordinates, i.e., the Laplacian distribution ( 2 ρ ( r ) ). In Figure 2, the planar representation mode indicates regions with less electron density using dashed lines. Such is the immediate vicinity of the carbon atom, where the molecule shows Lewis acid character. The area enclosed by continuous lines indicates electron density accumulation, suggesting Lewis base character. The figure also shows lines of zero gradient, which separate atomic basins, i.e., the attraction spheres of the atomic nuclei. The surface of the zero Laplacian value also separates areas of electron accumulation and Lewis acid character regions. It is worth mentioning that, although not adjacent to each other, the exchange between the two oxygen atom basins (i.e., the δ (O,O) indicated and is somewhat related to the bond order delocalization index) is surprisingly high: 0.39. Bader and co-workers obtained a value of 0.38 at the HF level and 0.31 at the CISD level for the same delocalization index [54].
In the next step, using domain-averaged Fermi hole analysis, we examine what the valence electron pairs look like that ultimately result in a cylindrically symmetrical electron density distribution. For our studies, we consider one of the oxygen atoms as the domain (fragment) and depict the eigenvalues of the Fermi hole integrated over it in Figure 3. Each eigenvector corresponds to an electron pair, and the associated eigenvalues indicate the proportion that each pair can be attributed to the given fragment. The σ bond is almost evenly distributed between carbon and oxygen (54%:46%), but interestingly, about 20% of the oxygen’s σ -symmetry lone pair is located on the other (i.e. CO) fragment. The π CO bond is significantly more polarized, being 66.5% attributable to the oxygen fragment. We also find two of them, which are degenerate and, of course, perpendicular to each other.
Inspecting the distribution of electron pairs between fragments, the significant exchange between the basins of the non-adjacent oxygen atoms ( δ = 0.39) becomes more plausible. It is apparent that the π CO orbital, though to a small extent, extends to the second oxygen atom. In addition, the π bonds located on the CO fragment also extend 4% to the O fragment, but the π CO orbital, indeed the n O σ -symmetry lone pair, also appears on the domain we examined at 5.5% and 1%, respectively.
Observing the eigenvectors, it is noticeable that, contrary to the cumulative representation, there is only one σ -symmetry lone pair on the oxygen atoms, but two completely identical π bonds (which are perpendicular to each other) are found. Therefore, it seems logical to assume that the “disappearing” π -symmetry lone pairs are the sources of the “new” π bonds. To confirm this, we use the natural localized molecular orbitals (NLMOs), which are also suitable for visualizing individual electron pairs. The NLMO can be derived from its parent NBO orbital, whose occupancy of less than two can be supplemented to two via intramolecular donor–acceptor interactions. The 1.62 occupancy of the π -symmetry lone pair centered on the oxygen of CO2 is primarily the result of significant electron density transferred to the adjacent CO fragment’s π * antibonding orbital, as observed in Figure 4. The resulting NLMO from this transition, characterized by significant interaction energy (135.7 kcal/mol), appears as a highly polarized π orbital, showing similarity to the corresponding DAFH (Domain Averaged Fermi Hole) eigenvector, although the electron pair shape is not entirely identical to that derived purely from the π NBO ( π ˜ CO ). Similarity can also be found between the NLMO of the σ -symmetry lone pair ( ñ O ) and its DAFH eigenvector.
The π electron system of CO2 can also be described by the NBO method as a three-center, four-electron bond. A combination of one-center hybrid orbitals based on the p z orbitals found on the two oxygen atoms and the carbon atom can yield two resonance states, which are equivalent to each other and have equal probability.
O 1 : C O 2 O 1 C : O 2 ( 50 % : 50 % )
Both states contain a π -symmetry lone pair on the oxygen and a carbon–oxygen π bond. However, averaging out the two boundary structures, we can also arrive at a 3c/4e bond ( ω -bond), in which the four π electrons are delocalized across the three atoms of CO2.
In conclusion, it can be stated that the cylindrically symmetrical structure of CO2 can be interpreted using Bader analysis, the NBO method, and DAFH analysis, while also resolving the apparent contradiction in the appearance of the two π bonds assigned to carbon.

3.2. Electron Structure of Small Molecules Containing Cumulated Double Bonds

We continue our investigations with carbonyl sulfide and carbon disulfide and compare these to ketene, diazomethane, and allene (see Figure 5). Replacing the oxygen atoms in carbon dioxide with sulfur causes significant changes in the electron density distribution of the molecule. For carbonyl sulfide, the molecular symmetry naturally changes from D h to C v . Here, the partial charge of the carbon atom is less positive than in the case of carbon dioxide, but for carbon disulfide, the direction of bond polarity reverses, as the carbon atom’s charge is negative and the sulfur atoms are positive. In the case of ketene, diazomethane, and allene, the methylene carbon carries a partial negative charge, especially pronounced at the terminal carbon atom of ketene.
The C=C and C=N bonds involving the terminal methylene group function as true double bonds, as there is no lone pair available that could be the source of a second double bond. The bond ellipticity is smallest in the case of allene ( ε CC = 0.377 ), followed by ketene ( ε CC = 0.754 ), and then diazomethane ( ε CN = 0.810 ). In the ketene CO and diazomethane CN fragments, the electron density distribution is already very close to cylindrical symmetry. Regarding diazomethane, it is also worth mentioning that according to Ponec and Cooper’s DAFH calculations, there is a polar σ CN and a less polar π CN bond between the CH2 and N2 fragments [55].
All the examined molecules containing cumulated double bonds show considerable similarity in both HOMO and LUMO orbitals. Except for allene, a rather high delocalization index (0.26–0.41) exists between non-adjacent atoms in the cumulated system, whereas in allene δ (C,C’) = 0.08, meaning that the three-center π electron system characteristic of the other molecules does not form, which can otherwise be characterized as a consequence of an adjacent π π * donor–acceptor interaction u sing the NBO method.
In the case of diazomethane, the probabilities were examined that are assignable to each resonance structure using NRT (natural resonance theory). Although the Lewis resonance structure carrying a positive charge on the internal nitrogen atom is the most probable (44%), surprisingly, the structure indicating a bond between carbon and the terminal nitrogen atom is also assigned a probability of 18% (see Figure 6).

3.3. Electron Structure of Metal Complexes Containing Ligands with Cumulated Double Bonds

The ratio of potential to kinetic energy calculated at the bond critical point provides additional information for characterizing a given chemical bond. If 1 < | V(r) | / G(r) < 2, it indicates a donor–acceptor-type bond. A value less than one suggests ionic interaction, while |V(r) | / G(r) > 2 implies a classical ’shared-shell’ covalent interaction, i.e., when the fragments forming the bond each contribute one electron to the bond.
The bond paths, Laplacian distribution, and kinetic energy density map of ketene and allene complexes containing the Ni(PH3)2 fragment can be seen in Figure 7. For ketene (besides the less stable η 1 -O coordination), two types of η 2 coordination modes can be envisioned. Of the two, the η 2 -(C,C) coordination proved to be more stable, with a free energy difference of 5.1 kcal/mol at the PBEPBE/def2-TZVP level. It can be determined that, in terms of electron density distribution, the η 2 -(O,C) complex of the two ketene complexes is more similar to the analogous Ni(PH3)2(CO2) (6) carbon dioxide complex, while the η 2 -(C,C) structure (2) resembles the Ni(PH3)2(H2CCCH2) allene complex (3). Similar to the analogous carbon dioxide complex, the absence of a bond path in the η 2 -(O,C) tautomer (1) suggests high kinetic energy density in the Ni-O region, which is confirmed by the kinetic energy density path calculated on the G ( r ) distribution. However, in the η 2 -(C,C) complexes, a bond path forms between the central metal atom and both coordinating carbon atoms, and the kinetic energy density maximum extends from the metal to the middle atom of the ligand.
Among the diazomethane complexes, two proved to be stable; the “side-on” η 2 -(N,N) coordination 4 and the κ -C coordination 5. The difference in free energy here is 8.4 kcal/mol in favor of the former complex when the two nitrogen atoms are involved in the coordination to the metal center. In our previous studies on nickel–carbonyl [30] and iron–carbonyl–phosphine [31] complexes, no η 2 -(N,N) coordination was predicted; however, κ -N coordination was also found, in contrast with the nickel–phosphine system used in this study.
The calculated structures of the Ni(PH3)2 fragment’s ketene, allene, and diazomethane complexes are shown in Figure 8, while their Bader parameters and the Wiberg bond indices (WBIs) of some important bonds are summarized in Table 1. In all cases, the coordination sphere of the nickel center adopts a planar arrangement with no out-of-plane angles for any of the ligands.
The nickel–phosphorus distances are good indicators that the third ligand can be involved in donor–acceptor interactions, transferring electron density to the σ NiP * orbitals. While the η 2 -(N,N) diazomethane and the η 2 -(C,C) allene complexes show no significant difference in this regard, it is interesting to compare the Ni(PH3)2(ketene) complexes. Switching the coordination type from η 2 -(O,C) to η 2 -(C,C) results in a noticeable elongation in the Ni-P bonds. That might indicate that the higher thermodynamic stability of complex 2 is due to the more favored overlap between the π CC orbitals of the coordinated ketene with the σ NiP * antibonding orbitals.
Based on the electron density and WBI calculated at the bond critical point, the strength of the interaction between nickel and the internal carbon atom in ketene complexes is similar to that in the carbon dioxide complex, while it is weaker in the allene complex. In the η 2 -(N,N) diazomethane complex, there is almost no difference between the two nitrogen atoms and the metal in terms of ρ BCP , but according to the Wiberg bond indices, the bond between nickel and the terminal nitrogen is stronger. In the η 1 -C complex, the strength of the interaction between the ligand and the metal is significantly weaker, in line with the greater distance between the metal and the coordinating carbon atom compared to the nickel–nitrogen distances in the η 2 -(N,N) complex. Interestingly, despite the less favored ligand coordination, the metal–carbon distance in complex 1 is shorter than in complex 2.
The ellipticity of the bonds slightly exceeds the ε value calculated for the Ni-C bond in the carbon dioxide complex. In the η 2 -(C,C) complexes, an extreme degree of flattening is observed in the bond between the central atom and the terminal carbon atom. In the η 2 -(N,N) diazomethane complex, there is also a difference in the ellipticity of the interactions between the two nitrogen atoms and nickel, but here the bond involving the terminal nitrogen does not show such an outstanding ε value.
The dominant donor and back-donation interactions of the Ni(PH3)2 fragment’s diazomethane complexes are presented in Figure 9. Based on the ETS-NOCV method (Extended Transition State theory combined with the Natural Orbitals of Chemical Valence), in both cases, the back-donation is characterized by greater interaction energy; however, in the η 1 -C complex, this energy value is significantly lower (−17.6 kcal/mol) compared to that calculated in the η 2 -(N,N) complex (−99.3 kcal/mol). This significant difference is well interpretable based on the NOCV eigenvectors. While in the 4 complex, the electron transition largely originates from the nickel d x 2 y 2 orbital and primarily reaches the diazomethane’s π symmetry LUMO, in the 5 complex, the primary source of electron transition is the Ni atom’s d z 2 orbital, which allows for a much weaker overlap with the LUMO of the H2CNN fragment.
The primary source of the donor interaction in both diazomethane complexes is the diazomethane HOMO orbital, while the main component of the NOCV eigenvector, which represents the destination of the electron transition, is the nickel 4s orbital. In the case of the ketene and allene complexes, the nature of the interactions, i.e., the deformation densities and the shapes of the NOCV orbitals, do not significantly differ from those calculated for Ni(PH3)2(H2CNN); hence, they are not separately highlighted as figures.
The total interaction energy for all η 2 complexes is very similar, slightly more than 60 kcal/mol (see Table 2). For exception, the carbon dioxide complex was found previously [50], where Δ E i n t = 51.9 kcal/mol. The exceptionally high orbital energy value calculated for the η 2 -(N,N) complex (—123 kcal/mol, at least 10 kcal/mol higher than in all other cases) is compensated by the high steric interaction energy (59.6 kcal/mol), which can be described as the sum of the electrostatic and Pauli interactions.
Although not to the same extent as in the 4 complex, but in all η 2 cases, the back-donation is much stronger than the donor interaction. The Hirshfeld analysis of the fragments also confirms this, as the charge on the ligand fragment is least negative in the case of the allene complex 3 (—0.283), exceeding the —0.3 value in other complexes. In contrast, in the η 1 -C diazomethane complex (5), there is hardly any charge difference between the two fragments. The minimal negative charge on the diazomethane fragment (—0.038) is consistent with the small difference calculated between the coordination and back-donation interaction energies.

3.4. Summary and Conclusions

This research focused on the detailed analysis of nickel–phosphine complexes, specifically targeting Ni(PH3)2 with various ligands such as ketene, allene, and diazomethane. The investigation employed QTAIM and ETS-NOCV analyses to dissect the orbital energy components and explore the Hirshfeld charges of the ligands, thereby revealing the intricate electronic structures and interactions within these complexes. The results disclosed herein can be summarized as follows.
  • Orbital Energy Decomposition and Ligand Charges: The study revealed that the primary source of donor interactions in diazomethane complexes is the diazomethane HOMO orbital, with the nickel 4s orbital being integral to electron transitions. The ketene and allene complexes showed interaction patterns consistent with those in the Ni(PH3)2(H2CNN) complex, suggesting similar electronic behaviors.
  • Interaction Energy Analysis: The total interaction energy for all η 2 complexes was found to be slightly over 60 kcal/mol. The carbon dioxide complex emerged as an outlier with a slightly lower interaction energy. This finding enriches our understanding of the energy dynamics within these complexes and their potential reactivity profiles.
  • Back-donation versus Donor Interactions: A significant observation was the predominance of back-donation over donor interactions in all η 2 cases. This was supported by Hirshfeld analysis, which revealed the nuances of charge distributions within the ligand fragments and between different complexes.
  • Charge Distribution and Interaction Strengths: The study noted that in the η 1 -C diazomethane complex, the charge difference between the two fragments was minimal, aligning with the subtle differences in coordination and back-donation interaction energies.
In conclusion, this research provides a comprehensive electronic and structural analysis of nickel–phosphine complexes with various ligands. The findings contribute significantly to the understanding of the electronic properties of these complexes, potentially guiding future research and applications in catalysis and organic synthesis. The consistency of interaction patterns across different ligands and the insights into their energy dynamics offer valuable information for the development of more efficient catalytic systems and the exploration of new synthetic routes in organic chemistry.

Supplementary Materials

The following supporting information containing the internal energies and Cartesian coordinates of all species occuring in this study can be downloaded at: www.mdpi.com/xxx/s1.

Author Contributions

Conceptualization and methodology, T.K.; calculations, T.R.K. and T.K., writing, review, and editing, T.R.K. and T.K. All authors have read and agreed to the published version of the manuscript.

Funding

Project no. RRF-2.3.1-21-2022-00009, titled National Laboratory for Renewable Energy, has been implemented with the support provided by the Recovery and Resilience Facility of the European Union within the framework of Programme Széchenyi Plan Plus; Project number 014_2023_PTE_RK/14 of the University of Pécs.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The author declare no conflicts of interest.

References

  1. Wang, W.; Wang, S.; Ma, X.; Gong, J. Recent advances in catalytic hydrogenation of carbon dioxide. Chem. Soc. Rev. 2011, 40, 3703–3727. [Google Scholar] [CrossRef] [PubMed]
  2. Arakawa, H.; Aresta, M.; Armor, J.N.; Barteau, M.A.; Beckman, E.J.; Bell, A.T.; Bercaw, J.E.; Creutz, C.; Dinjus, E.; Dixon, D.A.; et al. Catalysis research of relevance to carbon management: Progress, challenges, and opportunities. Chem. Rev. 2001, 101, 953–996. [Google Scholar] [CrossRef] [PubMed]
  3. Sordakis, K.; Tang, C.; Vogt, L.K.; Junge, H.; Dyson, P.J.; Beller, M.; Laurenczy, G. Homogeneous catalysis for sustainable hydrogen storage in formic acid and alcohols. Chem. Rev. 2018, 118, 372–433. [Google Scholar] [CrossRef] [PubMed]
  4. Kégl, T. Computational aspects of carbon dioxide coordination to transition metals. Glob. J. Inorg. Chem. 2010, 2, 286–303. [Google Scholar]
  5. Jessop, P.G.; Ikariya, T.; Noyori, R. Homogeneous hydrogenation of carbon dioxide. Chem. Rev. 1995, 95, 259–272. [Google Scholar] [CrossRef]
  6. Leitner, W. Carbon dioxide as a raw material: The synthesis of formic acid and its derivatives from CO2. Angew. Chem. Int. Ed. Engl. 1995, 34, 2207–2221. [Google Scholar] [CrossRef]
  7. Farlow, M.W.; Adkins, H. The hydrogenation of carbon dioxide and a correction of the reported synthesis of urethans. J. Am. Chem. Soc. 1935, 57, 2222–2223. [Google Scholar] [CrossRef]
  8. Inoue, Y.; Izumida, H.; Sasaki, Y.; Hashimoto, H. Catalytic fixation of carbon dioxide to formic acid by transition-metal complexes under mild conditions. Chem. Lett. 1976, 5, 863–864. [Google Scholar] [CrossRef]
  9. Denise, B.; Sneeden, R. Hydrogenate CO2. Chemtech 1982, 12, 108–112. [Google Scholar]
  10. Denise, B.; Sneeden, R. Hydrocondensation of CO2. J. Organomet. Chem. 1981, 221, 111–116. [Google Scholar] [CrossRef]
  11. Koinuma, H.; Kawakami, F.; Kato, H.; Hirai, H. Hydrosilylation of carbon dioxide catalysed by ruthenium complexes. J. Chem. Soc. Chem. Commun. 1981, 213–214. [Google Scholar] [CrossRef]
  12. Suess-Fink, G.; Reiner, J. Anionische rutheniumcluster als katalysatoren bei der hydrosilylierung von kohlendioxid. J. Organomet. Chem. 1981, 221, C36–C38. [Google Scholar] [CrossRef]
  13. Deglmann, P.; Ember, E.; Hofmann, P.; Pitter, S.; Walter, O. Experimental and theoretical investigations on the catalytic hydrosilylation of carbon dioxide with ruthenium nitrile complexes. Chem.–Eur. J. 2007, 13, 2864–2879. [Google Scholar] [CrossRef] [PubMed]
  14. Hoberg, H.; Schaefer, D. Nickel(0)-induzierte C-C-verknüpfung zwischen kohlendioxid und ethylen sowie mono- oder di-substituierten alkenen. J. Organomet. Chem. 1983, 251, c51–c53. [Google Scholar] [CrossRef]
  15. Schubert, G.; Pápai, I. Acrylate formation via metal-assisted C- C coupling between CO2 and C2H4: Reaction mechanism as revealed from density functional calculations. J. Am. Chem. Soc. 2003, 125, 14847–14858. [Google Scholar] [CrossRef] [PubMed]
  16. Pápai, I.; Schubert, G.; Mayer, I.; Besenyei, G.; Aresta, M. Mechanistic details of nickel (0)-assisted oxidative coupling of CO2 with C2H4. Organometallics 2004, 23, 5252–5259. [Google Scholar] [CrossRef]
  17. Li, J.; Jia, G.; Lin, Z. Theoretical studies on coupling reactions of carbon dioxide with alkynes mediated by nickel(0) complexes. Organometallics 2008, 27, 3892–3900. [Google Scholar] [CrossRef]
  18. Tidwell, T.T. Ketenes; John Wiley & Sons: Hoboken, NJ, USA, 2006. [Google Scholar]
  19. Borrmann, D. Methoden zur Herstellung und Umwandlung von Ketenen; G. Thieme Verlag Stuttgart: Stuttgart, Germany, 1968; Volume 7, pp. 94–95. [Google Scholar]
  20. Tidwell, T.T. Ketene chemistry: The second golden age. Acc. Chem. Res. 1990, 23, 273–279. [Google Scholar] [CrossRef]
  21. Tidwell, T.T. The first century of ketenes (1905–2005): The birth of a versatile family of reactive intermediates. Angew. Chem. Int. Ed. 2005, 44, 5778–5785. [Google Scholar] [CrossRef]
  22. Tidwell, T.T. Ketene chemistry after 100 years: Ready for a new century. Eur. J. Org. Chem. 2006, 2006, 563–576. [Google Scholar] [CrossRef]
  23. Seikaly, H.R.; Tidwell, T.T. Addition reactions of ketenes. Tetrahedron 1986, 42, 2587–2613. [Google Scholar] [CrossRef]
  24. Allen, A.D.; Tidwell, T.T. New Directions in Ketene Chemistry: The Land of Opportunity. Eur. J. Org. Chem. 2012, 2012, 1081–1096. [Google Scholar] [CrossRef]
  25. Staudinger, H.; Kupfer, O. Über Reaktionen des Methylens. III. Diazomethan. Ber. Dtsch. Chem. Ges. 1912, 45, 501–509. [Google Scholar] [CrossRef]
  26. McCarney, C.C.; Ward, R.S. An improved method for the preparation of monoalkylketens. J. Chem. Soc. Perkin Trans. 1 1975, 1600–1603. [Google Scholar] [CrossRef]
  27. Ungvári, N.; Kégl, T.; Ungváry, F. Octacarbonyl dicobalt-catalyzed selective carbonylation of (trimethylsilyl)diazomethane to obtain (trimethylsilyl)ketene. J. Mol. Catal. A Chem. 2004, 219, 7–11. [Google Scholar] [CrossRef]
  28. Ungvári, N.; Fördős, E.; Kégl, T.; Ungváry, F. Reactions of triphenylphosphane-substituted ethoxycarbonylcarbene-bridged dicobalt carbonyl complexes with carbon monoxide or 13CO: An experimental and theoretical study. Inorg. Chim. Acta 2009, 362, 1333–1342. [Google Scholar] [CrossRef]
  29. Ungvári, N.; Fördős, E.; Balogh, J.; Kégl, T.; Párkányi, L.; Ungváry, F. Triphenylphosphane-modified cobalt catalysts for the selective carbonylation of ethyl diazoacetate. Organometallics 2010, 29, 3837–3851. [Google Scholar] [CrossRef]
  30. Barcs, B.; Kollár, L.; Kégl, T. Density Functional Study on the Mechanism of Nickel-Mediated Diazo Carbonylation. Organometallics 2012, 31, 8082–8097. [Google Scholar] [CrossRef]
  31. Kégl, T.R.; Kollár, L.; Kégl, T. DFT Study on the Mechanism of Iron-Catalyzed Diazocarbonylation. Molecules 2020, 25, 5860. [Google Scholar] [CrossRef]
  32. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 16 Revision C.01; Gaussian Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  33. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef]
  34. Weigend, F.; Ahlrichs, R. Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. Phys. Chem. Chem. Phys. 2005, 7, 3297–3305. [Google Scholar] [CrossRef] [PubMed]
  35. Glendening, E.D.; Badenhoop, J.K.; Reed, A.E.; Carpenter, J.E.; Bohmann, J.A.; Morales, C.M.; Landis, C.R.; Weinhold, F. NBO 7.0; Theoretical Chemistry Institute, University of Wisconsin: Madison, MA, USA, 2018. [Google Scholar]
  36. Bader, R.F.W. Atoms in Molecules—A Quantum Theory; Oxford University Press: Oxford, UK, 1990. [Google Scholar]
  37. Keith, T.A. AIMAll, Version 19.10.12; TK Gristmill Software: Overland Park, KS, USA, 2019. Available online: https://aim.tkgristmill.com/(accessed on 17 March 2020).
  38. ADF2023; SCM, Theoretical Chemistry. Vrije Universiteit: Amsterdam, The Netherlands, 2023. Available online: http://www.scm.com (accessed on 1 October 2023).
  39. Mitoraj, M.; Michalak, A. Donor-Acceptor Properties of Ligands from the Natural Orbitals for Chemical Valence. Organometallics 2007, 26, 6576–6580. [Google Scholar] [CrossRef]
  40. Mitoraj, M.; Michalak, A. Natural orbitals for chemical valence as descriptors of chemical bonding in transition metal complexes. J. Mol. Model. 2007, 13, 347–355. [Google Scholar] [CrossRef] [PubMed]
  41. Mitoraj, M.P.; Michalak, A.; Ziegler, T. A Combined Charge and Energy Decomposition Scheme for Bond Analysis. J. Chem. Theory Comput. 2009, 5, 962–975. [Google Scholar] [CrossRef] [PubMed]
  42. Mitoraj, M.P.; Zhu, H.; Michalak, A.; Ziegler, T. On the origin of the trans-influence in square planar d8-complexes: A theoretical study. Int. J. Quantum Chem. 2009, 109, 3379–3386. [Google Scholar] [CrossRef]
  43. Mitoraj, M.P.; Michalak, A. σ-donor and π-acceptor properties of phosphorus ligands: An insight from the natural orbitals for chemical valence. Inorg. Chem. 2010, 49, 578–582. [Google Scholar] [CrossRef] [PubMed]
  44. Cukrowski, I.; de Lange, J.H.; Mitoraj, M. Physical nature of interactions in ZnII complexes with 2, 2’-bipyridyl: Quantum theory of atoms in molecules (QTAIM), interacting quantum atoms (IQA), noncovalent interactions (NCI), and extended transition state coupled with natural orbitals for chemical valence (ETS-NOCV) comparative studies. J. Phys. Chem. A 2014, 118, 623–637. [Google Scholar] [PubMed]
  45. Xu, S.; Li, M.; Pei, G.; Zhao, P.; Zhao, X.; Wu, G.; Kong, C.; Yang, Z.; Ehara, M.; Yang, T. Stabilities, Electronic Structures, and Bonding Properties of 20-Electron Transition Metal Complexes (Cp) 2TMO and their One-Dimensional Sandwich Molecular Wires (Cp = C5H5, C5(CH3)H4, C5(CH3)5; TM = Cr, Mo, W). J. Phys. Chem. A 2021, 125, 721–730. [Google Scholar] [CrossRef] [PubMed]
  46. Poater, J.; Vermeeren, P.; Hamlin, T.A.; Bickelhaupt, F.M.; Solà, M. On the existence of collective interactions reinforcing the metal-ligand bond in organometallic compounds. Nat. Commun. 2023, 14, 3872. [Google Scholar] [CrossRef] [PubMed]
  47. Parr, J.M.; White, A.J.; Crimmin, M.R. Calcium-stabilised transition metal bis (formyl) complexes: Structure and bonding. Chem. Commun. 2023, 59, 9840–9843. [Google Scholar] [CrossRef]
  48. Katlenok, E.A.; Kuznetsov, M.L.; Cherkasov, A.V.; Kryukov, D.M.; Bokach, N.A.; Kukushkin, V.Y. Metal-involved C... d z2-Pt II tetrel bonding as a principal component of the stacking interaction between arenes and the platinum (ii) square-plane. Inorg. Chem. Front. 2023, 10, 3916–3928. [Google Scholar] [CrossRef]
  49. Gironés, X.; Ponec, R.; Roithová, J. Program WinFermi, version 2.0; Prague, Czech Republic, 2002.
  50. Kégl, T.; Ponec, R.; Kollár, L. Theoretical Insights into the Nature of Nickel Carbon Dioxide Interactions in Pt(PH3)2(η2-CO2). J. Phys. Chem. A 2011, 115, 12463–12473. [Google Scholar] [CrossRef] [PubMed]
  51. Kégl, T.R.; Carrilho, R.M.; Kégl, T. Theoretical insights into the electronic structure of nickel (0)-diphosphine-carbon dioxide complexes. J. Organomet. Chem. 2020, 924, 121462. [Google Scholar] [CrossRef]
  52. Pálinkás, N.; Kollár, L.; Kégl, T. Nature of the Metal-Ligand Interactions in Complexes M (PH3) 2 (η2-L)(M = Ni, Pd, Pt; L = CO2, COS, CS2): A Theoretical Study. ChemistrySelect 2017, 2, 5740–5750. [Google Scholar] [CrossRef]
  53. Greenwood, N.N.; Earnshaw, A. Chemistry of the Elements, 2nd ed.; Elsevier: Amsterdam, The Netherlands, 1995. [Google Scholar]
  54. Fradera, X.; Austen, M.A.; Bader, R.F. The Lewis model and beyond. J. Phys. Chem. A 1999, 103, 304–314. [Google Scholar] [CrossRef]
  55. Ponec, R.; Cooper, D.L. Anatomy of bond formation. Bond length dependence of the extent of electron sharing in chemical bonds from the analysis of domain-averaged Fermi holes. Faraday Discuss. 2007, 135, 31–42. [Google Scholar] [CrossRef]
Figure 1. The generally assumed coordination modes of CO2 to transition metals.
Figure 1. The generally assumed coordination modes of CO2 to transition metals.
Molecules 29 00324 g001
Figure 2. Structure and NPA charges of carbon dioxide calculated at the PBEPBE/def2-TZVP level; HOMO and LUMO orbitals; Laplacian distribution in the plane of the molecule, along with the 2 ρ ( r ) = 0 envelope; and the MESP and electron density at the bond critical point, as well as the delocalization indices and bond ellipticity.
Figure 2. Structure and NPA charges of carbon dioxide calculated at the PBEPBE/def2-TZVP level; HOMO and LUMO orbitals; Laplacian distribution in the plane of the molecule, along with the 2 ρ ( r ) = 0 envelope; and the MESP and electron density at the bond critical point, as well as the delocalization indices and bond ellipticity.
Molecules 29 00324 g002
Figure 3. The DAFH eigenvectors of the CO2 O fragment, each corresponding to an electron pair, and their associated eigenvalues.
Figure 3. The DAFH eigenvectors of the CO2 O fragment, each corresponding to an electron pair, and their associated eigenvalues.
Molecules 29 00324 g003
Figure 4. The three-center four-electron π electron system of carbon dioxide (left); the natural localized orbitals (NLMOs); and the dominant donor–acceptor NBO interaction in the case of the NLMO functioning as the second π bond.
Figure 4. The three-center four-electron π electron system of carbon dioxide (left); the natural localized orbitals (NLMOs); and the dominant donor–acceptor NBO interaction in the case of the NLMO functioning as the second π bond.
Molecules 29 00324 g004
Figure 5. The NPA charge distribution, HOMO (left) and LUMO (right) orbitals, and Laplace distribution with the 2 ρ ( r ) = 0 envelope surface; the MESP and electron density at the bond critical point value; as well as the delocalization indices and bond ellipticities of carbonyl sulfide (a), carbon disulfide (b), ketene (c), diazomethane (d), and allene (e).
Figure 5. The NPA charge distribution, HOMO (left) and LUMO (right) orbitals, and Laplace distribution with the 2 ρ ( r ) = 0 envelope surface; the MESP and electron density at the bond critical point value; as well as the delocalization indices and bond ellipticities of carbonyl sulfide (a), carbon disulfide (b), ketene (c), diazomethane (d), and allene (e).
Molecules 29 00324 g005
Figure 6. The Natural Resonance Structures of diazomethane and their associated NRT probabilities.
Figure 6. The Natural Resonance Structures of diazomethane and their associated NRT probabilities.
Molecules 29 00324 g006
Figure 7. Laplacian (top) and kinetic energy density (bottom) maps of the Ni(PH3)2(ketene) and Ni(PH3)2(allene) complexes. Red dots indicate critical points on the kinetic energy density paths.
Figure 7. Laplacian (top) and kinetic energy density (bottom) maps of the Ni(PH3)2(ketene) and Ni(PH3)2(allene) complexes. Red dots indicate critical points on the kinetic energy density paths.
Molecules 29 00324 g007
Figure 8. The PBEPBE/def2-TZVP level’s calculated structures of the Ni(PH3)2(ketene), Ni(PH3)2(allene), and Ni(PH3)2(diazomethane) complexes.
Figure 8. The PBEPBE/def2-TZVP level’s calculated structures of the Ni(PH3)2(ketene), Ni(PH3)2(allene), and Ni(PH3)2(diazomethane) complexes.
Molecules 29 00324 g008
Figure 9. The π -acceptor (top) and π -donor (bottom) NOCV deformation densities of the Ni(PH3)2(CH2NN) 5 (left) and Ni(PH3)2(NNCH2) 4 (right) complexes and the corresponding complementary NOCV orbitals.
Figure 9. The π -acceptor (top) and π -donor (bottom) NOCV deformation densities of the Ni(PH3)2(CH2NN) 5 (left) and Ni(PH3)2(NNCH2) 4 (right) complexes and the corresponding complementary NOCV orbitals.
Molecules 29 00324 g009
Table 1. Wiberg bond indices and QTAIM parameters calculated at the bond critical point for Ni-PH3 ketene, allene, and diazomethane complexes.
Table 1. Wiberg bond indices and QTAIM parameters calculated at the bond critical point for Ni-PH3 ketene, allene, and diazomethane complexes.
Bond ρ (r) ε V G WBIBond ρ (r) ε WBI
Ni(PH3)2(OCCH2)Ni-C10.1190.2741.5180.465C1-O0.3630.0731.496
(1) C1-C20.3280.3891.784
Ni(PH3)2(H2CCO)Ni-C20.1190.2571.4740.442C1-C20.2940.2221.381
(2)Ni-C10.1031.5431.3690.332C2-O0.4230.0081.891
Ni(PH3)2(H2CCCH2)Ni-C20.1180.3441.4740.386C1-C20.3030.2251.481
(3)Ni-C10.1018.5331.3180.353C2-C30.3420.3021.946
Ni(PH3)2(NNCH2)Ni-N20.1120.6631.2260.315N1-N20.4720.0631.730
(4)Ni-N10.1160.3321.3040.525N2-C0.3230.5671.456
Ni(PH3)2(H2CNN)Ni-C0.0970.0431.3760.264C-N10.2700.4521.228
(5) N1-N20.5960.0252.421
Table 2. Ziegler–Rauk ETS-NOCV analysis of Ni(PH3)2(OCCH2) (1), Ni(PH3)2(H2CCO) (2), Ni(PH3)2(H2CCCH2) (3), Ni(PH3)2(NNCH2) (4), and Ni(PH3)2( η 1 -H2CNN) (5) complexes with orbital energy decomposition and Hirshfeld charges of the ligands ( Q H ). The NOCV eigenvalues are given in parentheses. The ratio of dominant back-donation interaction is given in square brackets, and energy values are in kcal/mol.
Table 2. Ziegler–Rauk ETS-NOCV analysis of Ni(PH3)2(OCCH2) (1), Ni(PH3)2(H2CCO) (2), Ni(PH3)2(H2CCCH2) (3), Ni(PH3)2(NNCH2) (4), and Ni(PH3)2( η 1 -H2CNN) (5) complexes with orbital energy decomposition and Hirshfeld charges of the ligands ( Q H ). The NOCV eigenvalues are given in parentheses. The ratio of dominant back-donation interaction is given in square brackets, and energy values are in kcal/mol.
Complex Δ E str * Δ E orb Δ E int Δ E d Δ E bd Q H
Ni(PH3)2(OCCH2)50.3−113.1−62.8−11.1 (0.31)−90.7 [%] (0.99)−0.362
Ni(PH3)2(H2CCO)41.1−105.2−64.2−15.1 (0.40)−67.2 [%] (0.80)−0.320
Ni(PH3)2(H2CCCH2)33.2−96.7−63.5−15.1 (0.37)−67.8 [%] (0.82)−0.283
Ni(PH3)2(NNCH2)59.6−123.0−63.4−12.1 (0.31)−99.3 [77%] (1.05)−0.365
Ni(PH3)2(H2CNN)12.2−39.2−27.1−14.2 (0.37)−17.6 [58%] (0.52)−0.038
* Total steric interaction: Δ E s t r = Δ E e l e c t r o s t a t i c + Δ E P a u l i .
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kégl, T.R.; Kégl, T. Comparative Analysis of Nickel–Phosphine Complexes with Cumulated Double Bond Ligands: Structural Insights and Electronic Interactions via ETS-NOCV and QTAIM Approaches. Molecules 2024, 29, 324. https://doi.org/10.3390/molecules29020324

AMA Style

Kégl TR, Kégl T. Comparative Analysis of Nickel–Phosphine Complexes with Cumulated Double Bond Ligands: Structural Insights and Electronic Interactions via ETS-NOCV and QTAIM Approaches. Molecules. 2024; 29(2):324. https://doi.org/10.3390/molecules29020324

Chicago/Turabian Style

Kégl, Tímea R., and Tamás Kégl. 2024. "Comparative Analysis of Nickel–Phosphine Complexes with Cumulated Double Bond Ligands: Structural Insights and Electronic Interactions via ETS-NOCV and QTAIM Approaches" Molecules 29, no. 2: 324. https://doi.org/10.3390/molecules29020324

APA Style

Kégl, T. R., & Kégl, T. (2024). Comparative Analysis of Nickel–Phosphine Complexes with Cumulated Double Bond Ligands: Structural Insights and Electronic Interactions via ETS-NOCV and QTAIM Approaches. Molecules, 29(2), 324. https://doi.org/10.3390/molecules29020324

Article Metrics

Back to TopTop