Open Access Article
Alejandra
Acuña
ab,
Nery
Villegas-Escobar
c,
Fernando
Gómez Zamorano
ab,
Nicole
Keim-Gyllen
b,
Sonia
Mallet-Ladeira
d,
Eddy
Maerten
a,
David
Madec
*a and
René S.
Rojas
*b
aLaboratoire Hétérochimie Fondamentale et Appliquée, Université de Toulouse, Centre National de la Recherche Scientifique (CNRS), UMR 5069, 118 Route de Narbonne, 31062 Toulouse Cedex 09, France. E-mail: david.madec@univ-tlse3.fr
bDepartamento de Química Inorgánica, Facultad de Química y de Farmacia, Pontificia Universidad Católica de Chile, Vicuña Mackenna 4860, 7820436, Macul, Santiago, Chile. E-mail: rrojasg@uc.cl
cDepartamento de Físico-Química, Facultad de Ciencias Químicas, Universidad de Concepción, Concepción 4070139, Chile
dInstitut de Chimie de Toulouse (UAR 2599), 118 Route de Narbonne, 31062 Toulouse Cedex 09, France
First published on 21st October 2025
The development of efficient catalysts for CO2 transformation into value-added products is crucial for sustainable chemistry. In this work, we report the synthesis and characterization of a new series of tetra- and penta-coordinated aluminum complexes supported by bis-amidine ligands featuring a rigid naphthalene bridge and their application as catalysts in the reaction of carbon dioxide with epoxides to form cyclic carbonates. Interestingly, while tetra- and penta-coordinated aluminum complexes bearing methyl groups are efficient when used with TBAI as a co-catalyst, their iodo analogs proved to be inefficient in mono-component catalysis, requiring the use of additional TBAI. Computational results of the proposed systems reveal the necessity of achieving a delicate balance between the nucleophilic properties of iodide atoms and the steric configuration of the complexes. Experimental and computational studies offer valuable insights into the role of ligand structure in aluminum-based catalysts, highlighting the potential of main-group metals in CO2 fixation.
As an alternative to transition metals, nowadays the chemistry of main-group element complexes, capable of mimicking transition metal behavior, is experiencing significant development. In particular, recent research has demonstrated that aluminum derivatives, the third most abundant element in the Earth's crust with low toxicity and low production costs, behave quite similarly to transition metal complexes and are able to activate strong bonds in small molecules (H2, N2, CO, and CO2).9,10
Since the choice of ligands is essential in each organometallic process, amidinate ligands, which are endlessly versatile in their potential structure and substitution patterns, have attracted attention for the stabilization of aluminum species.11 Indeed, they present high basicity with a pKa value of approximately 12 and therefore exhibit good affinity for protons and Lewis acids.12,13 Amidinate aluminum complexes have shown promising catalytic activity in various processes, such as ring-opening polymerization, alkyne hydroboration and the conversion of CO2 to cyclic carbonates.11k,l,n
In recent years, our research group has reported a series of aluminum amidinate complexes as catalysts for synthesizing cyclic carbonates from CO2 and epoxides.14–16 The first symmetric aluminum bimetallic complex I (a double aluminum amidinate, see Fig. 1) proved to be an excellent catalyst (conversions between 62 and 100% in 24 h) in the presence of tetrabutylammonium iodide (TBAI) as a co-catalyst.14 Bi- and tri-metallic aluminum-amidinate complexes were also reported (Fig. 1, II and III, respectively) as efficient catalysts for the formation of cyclic carbonates from terminal epoxides and CO2 at 50 °C under 1 bar of CO2 pressure without solvent.15 In 2021, a pentacoordinate bisamidinate iodide aluminum complex IV (Fig. 1) was reported, which is one of the first examples of a neutral non-zwitterionic aluminum complex used as a single-component catalyst for the synthesis of carbonates (catalytic activity at 80 °C and 1 bar of CO2, without a co-catalyst).16
![]() | ||
| Fig. 1 Double aluminum amidinate, bi- and tri-aluminum amidinate and pentacoordinate bis-amidinate complexes described by our group. | ||
As part of our ongoing interest in designing new aluminum complexes for catalytic applications, this work reports the synthesis and structural characterization of a new series of tetra- and penta-coordinated aluminum complexes with amidinate ligands featuring a rigid naphthalene bridge (1–7, Fig. 1). The performance and reactivity of these complexes were evaluated as catalysts in the reaction of carbon dioxide with epoxides to form cyclic carbonates. The influence of the substitution pattern on the amidinate ligand was explored to understand how these variations affect the catalytic activity and selectivity of the corresponding aluminum complexes.
![]() | ||
| Scheme 2 Synthesis of penta-coordinated aluminum (1 and 2) and tetra-coordinated bis-aluminum (3 and 4) complexes. | ||
The bis-aluminum complex 3 was synthesized using 1 equivalent of L1H2 and two equivalents of Al(CH3)3 in dry toluene at room temperature for 1 hour. Complex 3 was obtained as orange crystals in 83% yield. The formation of the corresponding bis-aluminum complexes was established by NMR spectroscopy and corroborated by XRD analysis (Fig. 2). The 1H NMR spectrum exhibited two signals in the negative region at −0.02 and −0.33 ppm, each integrating for 6 H, indicating that each aluminum atom is tetracoordinated with one amidinate ligand and two CH3 groups in its structure. In the 13C NMR spectrum, the signals of the methyl groups coordinated to the aluminum atoms appeared at −5.41 and −9.16 ppm.
Finally, the analogous aluminum complex 4 was obtained when L2H2 is mixed with two equivalents of Al(CH3)3. The NMR data are very similar to those of complex 3 with signals at −0.33 and −0.75 ppm that integrate for 6 H in the 1H NMR spectrum.
To gain more insight into these electronic structures, quantum chemical calculations were performed. Each molecular structure was optimized by employing density functional theory (DFT) calculations at the ωB97XD/def2-TZVP level of theory. Geometric indexes, such as τ5 and τ4, were computed at the DFT level. Molecular electrostatic potentials (MEPs) were calculated for all systems and are shown in Fig. 3. At first glance, this analysis underscores notable variations in electronic distribution that can be attributed to the effects of substitution either at the aluminum center or at the ligand.
![]() | ||
| Fig. 3 Molecular electrostatic potential (MEP) for all complexes (isosurface 0.01 au). NPA charges (qi) for selected atoms are displayed. | ||
The DFT-optimized molecular structures of complexes 1 and 2 confirmed a pentacoordinate Al center with τ5 indexes17 of 0.35 and 0.28, suggesting a distorted square pyramidal geometry at the Al centers. The computed equilibrium Al–C distances (bond orders) were 1.967 Å (0.56) and 1.965 Å (0.58), for complexes 1 and 2, respectively. Fig. 3 reveals that the Al atom in complex 1 exhibits a depletion of the electron density shown as a more positive atomic charge, while the C atom exhibits a more pronounced negative charge compared to complex 2. This observation is further supported by the comparative analysis of bond distances, bond orders, and atomic charges. Specifically, the longer bond distance and lower bond order, coupled with the greater charge separation of the Al–C bond in complex 1, unequivocally indicate a more polar Al–C covalent bond in this system compared to complex 2. Besides the small differences observed, these indexes still highlight the effect of the different R groups in both ligands. The τ4 index18 computed from the XRD and DFT structures was 0.95 and 0.94, which supports the use of the DFT structures for structural comparison. However, for a fair comparison, only the DFT-derived indexes are discussed. The τ4 indexes computed from the DFT structures were found to be 0.95 and 0.96 for complexes 3 and 4, respectively, confirming the distorted tetrahedral structure at Al. Interestingly, the DFT structures of the aluminum complexes reported slightly different Al–C bond distances for each Al center in complex 3 of 1.963 Å and 1.957 Å, respectively. However, the respective bond order remains similar for both distances around 0.67. The indistinguishable bond distances of N–C in the NCN unit and almost equal N–Al bond distances proved the formation of a delocalized system.
In a previously reported article, we have established that THE incorporation of a nucleophile (iodide) at the aluminum center could be an interesting strategy to obtain a single-component catalyst (avoiding the use of ammonium iodide as a cocatalyst). Indeed, the reaction mechanism has been shown to require an epoxide–iodide exchange at the aluminum atom, with an activation free energy of ∼20 kcal mol−1 for styrene oxide.16
For the pentacoordinate aluminum complexes, the formation of the respective iodide derivatives was achieved by reacting the corresponding complex (1 or 2) with one equivalent of I2 using a method previously described (Scheme 3).
Compounds 5 and 6 were obtained as brown solids in 95% and 96% yields, respectively. The products were fully characterized by NMR spectroscopy. These reactions were easily monitored by 1H NMR, with the disappearance of the signal corresponding to the methyl group at the aluminum center.
The NMR patterns of both iodo-complexes 5 and 6 again indicate molecular symmetry.
The optimized DFT structures of iodine-containing complexes (complexes 5 and 6) show Al–I bond distances of 2.565 Å (a bond order of 0.82) for complex 5 and a slightly shorter bond length of 2.547 Å (a bond order of 0.85) for complex 6. Notably, the NBO charges at the iodine centers remain consistent, with only minor variations in the electron population at the aluminum center. This significant electron density difference between Al and I suggests a polarized covalent bond, an essential characteristic for their function as single-component catalysts (vide infra). The MEP of complexes 5 and 6 (Fig. 3) indicates charge accumulation on the iodine atom, with delocalization across the ligand backbone. This delocalization is more pronounced in complex 6, extending between the naphthalene and phenyl groups due to the complex's relatively planar conformation.
In the case of the bis-aluminum complexes 3 and 4, attempts to replace the methyl groups with iodide using I2 with the present procedure were unsuccessful in our hands, resulting in mixtures of products that were difficult to identify. As an alternative synthetic route, BI3 was used as a source of iodides to carry out the substitution. The reaction between 3 and 2 equivalents of BI3 was carried out in benzene at room temperature for 1 hour (Scheme 4), resulting in the isolation of compound 7 as yellow crystals (91% yield). As before, the evident disappearance of signals in the negative region of the 1H NMR spectrum was a good indicator of reaction completion. Complex 7 presents an extremely low solubility in common solvents, preventing us from acquiring the corresponding 13C NMR spectrum. However, the formation of the bis-(iodo)complex was established by DRX analysis (Fig. 4), which revealed that the aluminum atoms exhibit a distorted tetrahedral geometry.
For the bimetallic complex 7, the MEP surface reveals a distinctive electronic distribution influenced by the presence of two Al–I bonds. By comparing the atomic charges at the Al and I centers, we gain insights into the polarity of these bonds.19 The results show that the Al–I bonds in complex 7 are approximately 20% less polar than those in complexes 5 and 6. Consequently, this reduced polarity of the two covalent bonds in complex 7 may suggest lower reactivity, limiting its efficacy as a single-component catalyst for cyclic carbonate formation.
Surprisingly, the replacement of methyl groups in complex 4 with iodides to form the halogenated complex led to a thermally unstable species. Indeed, at low temperature (−78 °C), the 1H spectrum indicates that the reaction proceeds similarly to the previously reported complexes 1–3. Above this temperature, decomposition occurred, making it impossible to isolate the iodocomplex.
![]() | ||
Fig. 5 Conversion of cyclic carbonates 9a–h employing catalysts 1–7. a Determined by 1H NMR spectroscopy of the crude reaction mixture. | ||
Generally, complexes 3 and 4 exhibited slightly better efficiency than complexes 1 and 2. For the mono-component catalysts 5, 6, and 7, the results were disappointing. Indeed, no catalytic activity could be detected for catalysts 5 and 7, and extremely low conversions were observed with catalyst 6 for three substrates (8, 3 and 9% conversion for substrates 9a, 9e, 9g, respectively). Specifically, the absence or the poor catalytic activity of compounds 5, 6 and 7 can be attributed to the poor polarity of the Al–I bonds as shown by the MEP surface (Fig. 3) that probably prevents (or limits) the liberation of the iodide atom. The X-ray data obtained in the case of complex 7 corroborate this hypothesis, with Al–I bond lengths rather short (between 2.4847 and 2.4996 Å) when compared to the one observed in a similar active one-component system (2.553 Å) with a similar ligand reported by our research group.16
To further understand the absence of reactivity of 5 and the limited reactivity of 6, we examined the geometries of these complexes by DFT when functioning as single-component catalysts. A critical comparison with the reactive system IV (Fig. 1) suggests that the pivotal step in these new complexes is the initial stage of the mechanism, involving the epoxide–iodide exchange at the aluminum center. The transition-state structure (TSIV+a) of complex IV during this exchange with styrene oxide (a) is depicted in Fig. 6. Concurrently, steric maps illustrating the computed free volume in the frontal regions of systems 5 and 6 are also presented in Fig. 6.
It is important to note that complexes IV and 5 differ only in the R substituent linked to the amidine moiety, with methyl in IV and tbutyl in complex 6. Although this substitution might seem minor, the bulkier tbutyl substituent in complex 5 significantly hampers its reactivity as a single-component catalyst. The steric maps in Fig. 7 further detail the approximate entry channel for the epoxide and the exit channel for the iodide atom, consistent with the transition state structure shown on the left. As illustrated, the tbutyl group on the amidine moiety, along with the methyl groups on the phenyl rings, obstructs both the entry channel for the epoxide and the exit channel for the iodide, thereby inhibiting the reaction mechanism.
![]() | ||
Fig. 7 Conversion of cyclic carbonates 9a–h employing catalysts 5, 6, and 7. a Determined by 1H NMR spectroscopy of the crude reaction mixture. | ||
Conversely, the low reactivity of complex 6 towards the selected epoxide can be attributed to more open entry and exit channels, coupled with a larger free volume. This is achieved by retaining the bulky tbutyl groups while positioning the methyl group on the phenyl rings at the para position. Computational analysis of the systems suggests that achieving an optimal balance between the nucleophilic nature of the iodide atoms and the steric demands of the complexes, which facilitates favorable entry and exit channels for the exchange of epoxide and iodide atoms, is crucial for the future design of single-component catalysts.
To investigate whether the lack of catalytic activity of complexes 5, 6, and 7 as one-component systems is due to the non-dissociation of the iodide ligand, we employed these complexes in the presence of TBAI as a cocatalyst, keeping all the other experimental conditions unchanged.
The conversions are depicted in Fig. 7, demonstrating that the simple addition of TBAI as a cocatalyst leads to an active system for every studied substrate. The efficiency of the mono- and bimetallic catalysts is very similar for every substrate. It is worth noting that the catalytic activity of iodoaluminum complexes is slightly lower than that of their methylated analogs.
These new catalytic systems show interesting reactivities, comparable to those previously described in the literature.11k–n,14–16 The recyclability of the catalyst, due to the type of catalysis (homogeneous catalysis) and the experimental conditions (without solvent, mol%…), is not possible.
21 were synthesized following reported procedures. 1D and 2D NMR spectra were recorded with the following spectrometers for 1H and 13C: Bruker Avance II 300 MHz, Avance III HD 400 MHz, and Avance I and II 500 MHz spectrometers. The chemical shifts were counted positively relative to the low field and expressed in parts per million (ppm). Mass spectrometric analyses were performed using a Maldi micro MX micro-mass spectrometer with a pyrene matrix (ratio of product/matrix
:
1/100). Melting points were measured with a capillary electrothermal Stuart SMP40 apparatus, and samples were prepared in the glovebox before the analysis. Single-crystal X-ray data were collected at low temperature (193(2) K) on a Bruker D8 VENTURE diffractometer equipped with a PHOTON III detector and using MoKα radiation (λ = 0.71037 Å). Phi and omega scans were performed for data collection, and an empirical absorption correction was applied.22 The structures were solved by the intrinsic phasing method (ShelXT)23 and refined by the full-matrix least-squares method on F2.24 All non-hydrogen atoms were refined with anisotropic displacement parameters. All hydrogen atoms were refined isotropically at calculated positions using a riding model, except for those of the water molecules in the structure of compound 3, which were not located in the difference Fourier maps. For 3, some electron density was difficult to model; therefore, the SQUEEZE tool25 of PLATON was used to remove the corresponding electron density from the intensity data.
Deposition numbers 2477633 (3) and 2477634 (7) contain the supplementary crystallographic data for this paper.
Complexes 1–7 were fully optimized at the ωB97XD/def2-TZVP level of theory, with the resulting structures confirmed by frequency calculations at the same level. Population and bond order analyses were conducted using the NBO3.1 program, interfaced with Gaussian16 Rev B.01. Bond orders were specifically calculated using the Mayer definition, recommended as a suitable method for characterizing bonds in metal complexes due to its incorporation of the overlap matrix in the bond order determination.26 Molecular electrostatic potentials were computed for all systems at the same level of theory.
(CH3)3), 29.5 (C(
H3)3), 19.7 (CH3), −4.3 (CH3–Al).
(CH3)3), 31.1 (C(
H3)3), 21.0 (CH3), −4.0 (CH3–Al).
(CH3)3), 29.8 (C(
H3)3), 20.7 (CH3), 20.5 (CH3), −5.4 (CH3–Al), −9.2 (CH3–Al). MS (Maldi-TOF) m/z: 643.21 ([M]+), 629.17 ([M − CH3]+), 615.16 ([M − 2(CH3)]+).
(CH3)3), 31.0 (C(
H3)3), 20.9 (CH3), −7.9 (CH3–Al), −11.0 (CH3–Al).
(CH3)3), 29.2 (C(
H3)3), 20.4 (CH3).
(CH3)3), 30.7 (C(
H3)3), 21.1 (CH3).
CCDC 2477633 and 2477634 contain the supplementary crystallographic data for this paper.27a,b
| This journal is © The Royal Society of Chemistry 2025 |