Abstract

Density functional theory (DFT) calculations were performed on a new family of mono- and bimetallic complexes, containing 4,8-(paracyclophane)-1,5-dihydro-s-indacene as the bridging ligand between the two metallic centers and different ancillary ligands. The s-indacene was blocked by substitution of the central benzene ring with the paracyclophane to obtain the syn-conformations. This would force the metallic centers to be close together. It is proposed, due to the calculated molecular and electronic properties such as the reactivity indexes, the delocalized nature of the s-indacenyl ligand, and the electron-rich metals, that these systems could be reactive in a catalytic reaction. The results indicate that the systems with Rh and Re holding ancillary ligands such as bipy and CO show the best properties to be active in a chemical reaction. In this sense, by the assessed geometrical and electronic properties, when compared with a previously reported system, these complexes could be candidates for the reduction of CO2 to oxalate.

1. Introduction

s-Indacene and as-indacene belong to a class of conjugated hydrocarbons characterized by having electrons (). Therefore, according to Hückel’s Rule, they should be antiaromatic [1]. Nevertheless, they have been for decades of significant interest for organic chemists as they have focused on them upon their diverse and fascinating reactivity. Due to the synthetic challenges that these kinds of compounds present and the questions that have arisen over their aromaticity and structure, the potential of these systems in different applications, such as semiconducting devices and catalysis, has been slowly explored.

Even when the unsubstituted molecules are unstable, their dianions and several substituted derivatives stabilize both sterically and electronically, resulting in thermally stable structures at room temperature [2].

Moreover, s-indacene may connect two or more organometallic centers, which may have different oxidation states, and allow interaction between them. Binuclear metal complexes bearing an as-indacenyl or s-indacenyl ligand were first reported in the literature in 1988 [3, 4]. A large variety of transition metals homobinuclear complexes have been described including Fe, Co, Ni, Ru, Rh, Re, and Mn. Also, heterobinuclear complexes with combinations of metals such as Rh-Ru, Fe-Ru, Fe-Rh, Ru-Rh, F-Co, and Cr-Ni have been reported [58]. The most striking feature of this family of binuclear complexes is the ability to delocalize charge among metal centers upon oxidation. In this sense, we believe that bimetallic complexes are useful compounds that might induce catalytic reactions.

On the other hand, literature reports on the selective conversion of CO2 to oxalate (C2O4)2- point to the occurrence of a two-metal-centers mechanism. Bi- or polymetallic compounds have been studied showing a wide variation of results [912]. The most important conclusion of these reports is that, for systems where the coordination of CO2 is not reversible, the ability to transform CO2 into other compounds depends on the extent of the charge transfer from the transition metal to CO2. If the metal transfers an electron, the formation of C-C bonds is possible and the rigidity of the bimetallic system seems to play an essential role in the formation of this bond [13, 14]. In this sense, interestingly, one of the most innovative works in the area of CO2 reduction reported by Angamuthu et al. in 2010 consists in the synthesis of a tetranuclear complex of Cu(I) with the ligand N-(2-mercaptopropyl)-N,N-bis(2-pyridylmethyl)amine, which spontaneously and selectively captures CO2 from the air [15]. Cu(I) gets oxidized to Cu(II) and a bridge of oxalate is formed, which links two binuclear complexes, forming a tetramer. Adding a lithium salt to the solution allows the quantitative precipitation of lithium oxalate. Finally, the Cu(I) complex can be recovered by electrochemical reduction. The crystallographic study of the Cu complex showed that the Cu-Cu distance in the complex was 5.43 Å. This distance is similar to the Rh-Rh distance observed in a syn binuclear complex synthesized by our group [(2,6-diethyl-4,8-dimethyl-s-indacenediide) Rh (COD)] [16] of about 5.5 Å.

The large majority of the examples of binuclear complexes with s-indacene reported in the literature show the two metals centers in an anticonfiguration. This is due mostly to the steric bulk imposed by the ancillary ligands. In this work, we propose that it is possible to avoid the anticonfiguration by blocking one of the faces of the s-indacene introducing a group that can effectively obstruct the access of the metals to both faces forcing them to adopt the syn configuration. This could be easily achieved by constructing a cyclophane structure over the plane of the rings as shown in Figure 1.

Therefore, in order to have the right template structure for coupling reactions in those compounds and to favor the reaction, the right Metal-Metal (M-M) distance, electronic configuration, and lability of the ligand coordinated to the metal center (ancillary ligands) are needed. Specifically, binuclear metal frameworks could play an important role as templates for these types of reactions, such as selective reductive coupling of two molecules of CO2 to oxalate.

Inspired by the work of Bouwman et al. and the versatile indacene derivatives, we decided to propose templates with synthetical feasibility, consisting in a family of binuclear metal complexes.

The aim of this work is to evaluate the effect of different transition metals (Ir, Rh, and Re), different ancillary ligands (bipy and CO), and -R groups (CH3 and CH2CH3) (Figure 1) over the electronic properties, assessing the Frontier Molecular Orbitals (FMOs) location, the geometrical parameters highlighting the M-M distance, the global reactivity indexes, the atomic metal charges, and the Fukui functions. This preliminary study could lead experimentalist to have a hint about which metallic s-indacenyl complexes, ancillary ligands, and peripheral substituents might show suitable properties for a variety of chemical applications.

2. Computational Details

Quantum chemical calculations were performed in order to further investigate the electronic structure of a family of compounds in the ground state. Density functional theory (DFT) as implemented in Gaussian 09 [17] software was used to perform the calculations of the studied systems. The Becke three-parameter hybrid functional combined with Lee-Yang-Parr correlation B3LYP [18, 19] was employed. The relativistic effects of the heavy metals were incorporated using the pseudopotential LANL2DZ [2022] and the Gaussian basis sets 6-31G(d) were used for the nonmetallic atoms [23]. With the aim of getting the minima points over the potential energy surface of the proposed compounds, the molecular structures were fully optimized without symmetry constriction and the frequency calculations were performed at the same level of theory to confirm that all the optimized systems are stationary minima points. Implicit solvation effects were incorporated using the polarized continuum model (PCM, ) for tetrahydrofuran [24]. The local reactivity indexes and the Fukui functions computed at the ground state geometry are also assessed and analyzed.

3. Results and Discussion

The optimized ground state minima of the series of compounds displayed in Figure 1, ligand (s-Ic) in the neutral and dianion states, and the mononuclear and binuclear coordination systems were obtained at the B3LYP/6-31G(d) + PCM () level of theory. Selected geometrical parameters of the corresponding minima are listed in Table 1.

The ligand s-Ic shows a planar structure and one plane face is blocked by substitution of the central benzene ring with a saturated alkyl chain formed by 10 -CH2- units. This substitution is a proposal for experimentalist to synthesize the compounds adopting exclusively the syn conformation after metallic coordination.

Geometrical parameters such as C-C and C-H of the s-Ic coordinating different metals (not shown in Table 1) do not change significantly when all the studied compounds are compared as some values show differences in the third decimal digit or in less proportion in the second decimal digit. This indicates that the different metals and ancillary ligands do not lead to significant changes in the structural geometry of the s-Ic.

Regarding the M-M distances, as observed in Table 1, the binuclear transition metals frameworks show the appropriate length between the two metal centers to act as a template for the selective reduction of two molecules of CO2 to oxalate. Specifically, the results show that the systems with Rh and Re holding ancillary ligands such as bipy and CO, containing -CH3 as -R substituent, present the most accurate M-M distance around 5.43 Å (Table 1). This distance corresponds to the Cu-Cu distance of the complex reported by Bouwman et al. that catalyzed the transformation of CO2 to oxalate.

The ancillary ligands contribute steric effects towards the distances of the bonds between the metal and the s-Ic, while as bigger the ancillary ligand is, the shorter the distances among the metal and the carbon atoms are (Table 1). In this sense, it is also worth mentioning that while the -R substituents in the s-Ic are changed, some differences over the M-M distances were observed; however, no clear trend is found.

These facts lead us to conclude that the M-M distance can be tuned by adjusting different ancillary ligands and substituents over the s-Ic ligand.

In terms of the electronic structure, as consequence of the formation of the s-Ic dianion species with different R groups, the HOMO-LUMO gaps get smaller in comparison with the gap of the neutral species (Figures S1-S2 in Supplementary Material available online at https://doi.org/10.1155/2017/9101720). For the isolated ligand, it is also observed that the FMOs structures present contribution of the s-Ic ligand, except for the LUMO+1 of the dianionic compounds which is localized over the alkyl chain binding across the six-edge ring of s-Ic (Figures S1-S2). The change of the R groups does not show transcendental effect over the FMOs.

In cases of the metallic complexes, the localization of the FMOs does not change considerably when the comparison is made between the monometallic and the bimetallic species (see complexes with Rh as example in Figure 2), as they are formed by contribution of the orbitals extended over the s-Ic and over the transition metals.

The energy levels diagrams display that the FMOs of all the mononuclear and binuclear complexes are more stable than the free substituted s-Ic ligand. It is worth noticing that all studied systems with ligands CO show more stable occupied orbitals than the system with bipy.

It was also observed that the different studied peripheral R substituents of the s-Ic ring (-H, -CH3, and -CH2CH3) show no contribution to the FMOs.

In general, all the complexes that have a bimetallic structure exhibit FMOs localized mainly over the ancillary ligands and the metal centers, although a slight variation of the FMO isosurfaces is observed specifically in the LUMO. This LUMO, in cases of voluminous ancillary ligands, is composed mainly for the ancillary ligands orbitals and the contribution of the metals decreases (Figures S3–S20 in the Supplementary Material).

In these systems, the mononuclear and the binuclear, the metals would have a dual role. The metal would be receiving electrons, therefore stabilizing the reduced fragments for the formation of the C-C bond between two CO2 molecules. This is supported by the fact that the HOMO is localized over the fragment containing the metal and the ancillary ligands (Figure 2 and Figures S3–S20 in Supplementary Material). The localization of the FMOs over these regions of the molecules suggests their high reactivity towards any catalytic activity.

To assess the atomic charges in the metals of the studied family of compounds, a Mülliken analysis was carried out and the values are reported in Table 2. The formation of the complexes involves a charge transfer from the s-Ic and ancillary ligands to the metals as the atomic charges get less positive with respect to the oxidation state of the metals (+1 for all metals in all cases). Significant charge transfer is observed for compounds containing Rh and CO as ancillary ligands, where it is known that the back-bonding effect of the CO ligand is present. In this sense, it is possible to state that in these systems the electron pair of each CO ligand donates more charge than it receives, which is evidenced in the slight negative charge over the metals. This behavior is followed by Re compounds containing CO which show atomic charges around 0.16 a.u. These values are significantly lower than the values obtained with Ir as metal center.

On the other hand, in cases of monometallic complexes, the charge transfer to the metals is bigger for the mononuclear anionic systems than for the neutral mononuclear ones, except the anionic mononuclear system containing Ir and bipy as ancillary ligands.

Taking into account these results, it is observed that Rh and Re bimetallic compounds might act as templates to carry out catalytic applications such as CO2 reduction reactions. These metals when forming part of these types of systems get rich in electron charge, which is in concordance with the fact that the occupied FMOs are localized over the metals. These facts lead us to consider these systems as better reducing agents than the compounds with Ir.

In order to rationalize the chemical reactivity, a set of global reactivity descriptors were computed in the framework of the conceptual density functional theory. The estimated chemical potential (), chemical hardness (), and electrophilicity () have been analyzed [2529]. This set of values is reported in Table 3. Moreover, the HOMO and LUMO energies and the HOMO-LUMO energy gaps are reported. For this analysis, we restricted the discussion to the complexes that are considered the best templates for catalytic reactions as stated above, that is, systems with Rh and Re, CO and bipy as ancillary ligands, and -CH3 as peripheral substituent.

Parr defined µ in (1) as the infinitesimal change of energy of the system with respect to the electron number N at a constant external potential of the nuclei (), that is, the potential created by the nuclei. is related with its electronegativity as and is associated with the feasibility of a system to exchange electron density with the environment at the ground state [2529].

is calculated using the second-order derivative of the chemical potential (see (2)) and can be interpreted as the resistance of a molecule to change the electron density in presence of other species.In numerical applications, these reactivity indexes are calculated following approximations using Koopmans’ theory and finite differences leading to (3). Here, I and A correspond to the ionization energy and electron affinity, while and correspond to the orbital energies of the HOMO and LUMO, respectively [2528]:The ω index as defined in (4) measures the tendency of a molecule to receive electronic charge from a donor species. Therefore, ω is considered as a sort of “electrophilicity power” [2729]:When bimetallic complexes are compared (see Table 3), the values of of systems with Rh and bipy (–2.76 eV) and complex with Re and bipy and CO (–2.99 eV) are higher than values presented by systems with Rh, Re, and CO as ancillary ligands (–3.47 eV and –4.01 eV, resp.). This suggests bigger reactivity of the first two named complexes towards electron donation. Hence, these complexes tend to react as source of electrons that enhance their possible behavior as catalytic templates for reduction reactions, such as the reduction of 2CO2 to .

It is worth mentioning that of the monometallic neutral species adopts values similar to of the bimetallic complexes, which indicates that these systems can also act as strong electron donors. Finally, the anionic mononuclear species display the highest contemplated in the range of to  eV.

The decreases in the order mononuclear neutral > bimetallic > anionic mononuclear are reasonable with respect to the impact of the global reactivity/polarization indexes. The obtained trend indicates that the last two named systems are more reactive. In case of , as expected, the anionic mononuclear systems show the lowest values and increase for the neutral and bimetallic compounds.

In order to clarify the dependency of the reactivity with the molecular structure, we also analyzed the condensed Fukui functions [25, 30, 31]. In Figure 3, we included the wavefunctions for nucleophilic () and electrophilic () attacks. Considering that molecules may react with sites that show high or low electron density to carry out a reaction, the plots depict that the metallic centers are able to follow a reaction pathway to carry out a reduction, as they have high density of electrons. These results are consistent with the known reactivity of ferrocene (and many other organometallic complexes) with any reducing species, where the reaction involves the formation of ferrocenium. Therefore, the reactivity of the herein studied molecules will depend on the interactions occurring in the metal of the complex with other species. Although these results are not a definitive evidence on how the mechanism of the chemical reactions might be, they are useful to understand the chemical behavior of the family of complexes.

4. Conclusions

A theoretical study of a new family of organometallic complexes was carried out. As it is observed, the results indicate that the bimetallic systems with Rh and Re holding ancillary ligands such as bipy and CO show the best properties to be active in a chemical reaction. This is supported by the fact that the most electron-rich portion of the molecule corresponds to the metallic centers and the s-Ic ligand. This is shown by the localization of the FMOs, the Mülliken charges, and the Fukui function plots. Furthermore, these complexes show the highest reactivity indexes. Specifically, the results display that the bimetallic complexes containing Rh and bipy and Re with bipy and CO would tend to be strong electron donors in a chemical reaction. This behavior is also observed for their monometallic neutral species.

Finally, in addition to all previously discussed results, it has been stated that an appropriate M-M distance around 5.43 Å is an important parameter for a bimetallic compound to be able to perform the reductive coupling of two molecules of CO2 to oxalate. All the studied complexes show an M-M distance around this value, which would suggest that these complexes might show activity in this type of reactions.

This study will continue employing more ligands and metals, including the formation of oxalate, and then experiments will be carried out.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors acknowledge the financial support from the following projects: Fondecyt (Grants 1141138 and 1161416), CONICYT, Operating Expenses (Grant no. 21110092), and Milenio (RC120001 and REDES150042).

Supplementary Materials

The structures of the neutral ligand are reported in Figure S1. The structures of the dianionic ligand are reported in Figure S2. The structures of the neutral monometallic, monoanionic monometallic and bimetallic system with Ir and bipyridine as ancillary ligand are reported in Figures S3–S5. The structures of the neutral monometallic, monoanionic monometallic and bimetallic system with Rh and bipyridine as ancillary ligand are reported in Figures S6–S8. The structures of the neutral monometallic, monoanionic monometallic and bimetallic system with Re and bipyridine and CO as ancillary ligands are reported in Figures S9–S11. The structures of the neutral monometallic, monoanionic monometallic and bimetallic system with Ir and CO as ancillary ligand are reported in Figures S12–S14. The structures of the neutral monometallic, monoanionic monometallic and bimetallic system with Rh and CO as ancillary ligand are reported in Figures S15–S17. The structures of the neutral monometallic, monoanionic monometallic and bimetallic system with Re and CO as ancillary ligand are reported in Figures S18–S20.

  1. Supplementary Material