Open Access Article
Huili Wang†
a,
Qiqi Zhang†b,
Haizhu Yu*b,
Kane Jacobc,
Azzedine Bousseksou*c,
Philippe Hapiot
d,
Jean-Luc Pozzo
*a and
Didier Astruc
*a
aUniversité de Bordeaux, ISM, UMR CNRS N°5255, 33405, Talence Cedex, France. E-mail: didier.astruc@u-bordeaux.fr
bDepartment of Chemistry and Center for Atomic Engineering of Advanced Materials, Anhui Province Key Laboratory of Chemistry for Inorganic/Organic Hybrid Functionalized Materials, Anhui University, Hefei, Anhui 230601, People's Republic of China. E-mail: yuhaizhu@ahu.edu.cn
cUniversité de Toulouse, LCC, UPR CNRS N°8241, 31077, Toulouse Cedex, France. E-mail: azzedine.bousseksou@lcc-toulouse.fr
dUniv. Rennes CNRS, ISCR (Institut des Sciences Chimiques de Rennes), UMR 6226, F-35000, Rennes, France
First published on 13th April 2026
The first Fe(0)-catalyzed alkyne carboxylation using CO2 is reported using Fe(0)–phenanthroline complexes, and DFT calculations suggest a spin crossover from a triplet to a quintet state during the reaction.
First-row, late transition metal complexes with strong field ligands such as cyclopentadienyl and CO mostly fulfill the 18-electron rule, and they are “closed-shell” systems.16,17 On the other hand, inorganic nitrogen and oxygen ligands exert weak fields, so this weak ligand field leaves a relatively lower energy (Δ) between the three lower-energy eg orbitals and the two higher-energy t2g d orbitals in the case of an octahedral or a pseudo-octahedral geometry. The low Δ value in inorganic complexes opens the possibility of different orbital occupations by electrons, leading to the possibilities of two spin states, low spin or high spin.18,19 Different reactivities, including catalytic behaviours, are expected from these two spin states, which has already been exploited as a key feature since the 2000s in bio-inorganic catalytic systems.20–23 More recently, various studies have demonstrated the effect of the spin state and spin crossover in catalysis,24–37 but not yet for carboxylation reactions. Complexes with the general formula Fe(II)(1,10-phenanthroline)Cl2 bearing bulky substituents at the 2,9 positions have been reduced by the Zhu group to substituted Fe(0)(phenanthroline) complexes showing catalytic activities for reactions of alkenes and alkynes.27,37–43 This family of Fe(0) complexes are used in this study as alkyne carboxylation catalysts. The Fe electronic structure therein is 4s24p63d8, i.e. the eight d electrons occupy the five d orbitals. The chloro ligands have now been replaced by solvent ligands in the reduction process, and, according to both Pauli and Hund's rules, the three lower d orbitals are doubly occupied with electrons of opposite spins, and the two higher orbitals are singly occupied with electrons of the same spin, corresponding to a triplet spin ground state (S = 1; 2S + 1 = 3). Zhu et al. have taken their catalytic results into account using a two-spin-state reactivity, the phenanthroline acting as a redox ligand that could possibly accept an electron from Fe(0) to yield a quintet-state Fe(I) complex of phenanthroline radical anion.38 Under these conditions, among the five d orbitals, three become singly occupied, whereas one electron occupies a ligand orbital (S = 2). Such spin delocalization was reported in 2010–2012 with other iron catalysts containing a redox-active ligand; however, the relationship with reactivity effects was not rationalized.38
Herein, we report, using this family of iron–phenanthroline complexes, the first Fe(0)-catalyzed alkyne carboxylation reaction utilizing CO2 (eqn (1)), and DFT calculations suggest spin crossover during the catalyzed carboxylation process.
![]() | (1) |
In order to obtain the Fe(0)(BTIPP) complex (BTIPP = 2,9-bis(2,4,6-triisopropylphenyl)-1,10-phenanthroline ligand), the Zhu group reduced the [Fe(II)(BTIPP)Cl2] complex using EtMgBr in THF and recorded the X-ray crystal structure of the Fe(II) precursor and a diolefin Fe(0) adduct.37,39
In the present work, the same and related Fe(II) precursors have been synthesized by stoichiometric reactions in THF between FeCl2 and the phenanthroline derivatives with substituents at the 2,9 positions (chloro, phenyl, mesityl and triisopropylphenyl), followed by precipitation from addition of excess hexane to a concentrated THF solution of the complexes and recrystallization from ether. The corresponding reduced black Fe(0) complexes have been generated in the same way as Zhu et al.37,39 by reduction of the Fe(II) precursor in DMF using EtMgBr in ether. Their cyclic voltammograms consistently show two successive irreversible cathodic waves at −1.2 and −2.1 V vs. Fc+/Fc on an Au disk electrode corresponding to the successive reductions of the first and second Cl ligands (SI). The resulting reactive Fe(0) catalysts that are very sensitive have not been isolated, but utilized in situ for carboxylation of various alkynes.
After reduction of these Fe(II) complexes to Fe(0) catalysts conducted using EtMgBr in ether/DMF, it is speculated that Fe is in the zero oxidation state with two or three solvent molecules (from the solvent used in the synthesis) in addition to the chelating L2 phenanthroline ligand, which would involve, respectively, 16-electron (2 weak solvent ligands) or 18-electron (3 weak solvent ligands) complexes. According to the DFT calculation results, in the presence of the sterically bulky BTIPP and DMF solvent, the coordination of three DMF molecules is unlikely. The two DMF molecules favorably adopt one η1- and one η2-coordination, and the Fe(0) complex is in a triplet spin state (Fig. 1).
![]() | ||
| Fig. 1 Relative energies (in kcal mol−1) of the BTIPP Fe(0) structures with coordination of two DMF ligands. The superscripts T, Q and S denote triplet, quintet and singlet states, respectively. | ||
Specifically, in the reduction step, the two chloro (X-type) ligands of the 14-electron Fe(II) precursors are replaced by the solvent:
![]() | (2) |
(L-type ligand, eqn (2)), yielding the Fe(N^N)L2 or Fe(N^N)L3 complexes.
During the catalytic process, the solvent ligands are readily displaced by alkyne and CO2 substrates. DFT calculations confirm that this ligand exchange is highly facile (Scheme 1).
![]() | ||
| Scheme 1 Gibbs free energy changes for the replacement of the DMF solvent by two alkynes (top) or an alkyne and CO2 (bottom). | ||
The alkyne carboxylation reactions catalyzed by the Fe(0) complexes were conducted using the four Fe(0) complexes at 25 °C and 80 °C, but it was found that the yields were approximately the same at both temperatures without any improvement at 80 °C, and therefore the yields obtained using different arylacetylene substrates at 25 °C are provided in Table 1. The best alkyne carboxylation yield (95%) was obtained using the bulkiest catalyst Fe(0)BTIPP with phenylacetylene, but, interestingly, all the other alkynes bulkier than phenylacetylene provided much lower yields, clearly demonstrating the predominance of steric effects over electronic substituent effects. With the other less bulky catalysts, the electron withdrawing substituents on the arene were beneficial to the reaction, in contrast to the arylacetylene bearing an electron-donating substituent. The catalysts bearing CF3 substituents led, due to steric bulk, to lower yields than those obtained with fluoroaryl substituents. The parent Fe(0)phenanthroline catalyst without substituents at the 1,9 positions gave very low yields (10% or less) with the various arylacetylenes, confirming the importance of the crucial steric effect in 2,9-disubstituted 1,10-phenanthroline positions observed with the Fe(0)BTIPP catalyst. With 1-butyne and 1-pentyne, the reaction provided modest yields; however, with the latter catalysts, mutual steric effects disfavored the reaction (Table 1). On the basis of the aforementioned experimental observations, DFT calculations were conducted on the modeling phenylacetylene carboxylation mechanism as well as the possible occurrence and effect of spin crossover during catalysis (see the SI for details).
As shown in Fig. 2, the catalytic cycle begins with the triplet bis-alkyne adduct TIn1, whose energy lies below that of its quintet counterpart QIn1 and all other states. A spin crossover from a triplet to a quintet state could then occur via the MECP (minimum energy crossing point), whose energy is 31.0 kcal/mol higher than that of TIn1 to generate its quintet state QIn1. After that, the dissociation of one alkyne from QIn1 to generate the mono-alkyne coordinated quintet intermediate QIn2, and this step is slightly exergonic (ΔG0 = −4.5 kcal mol−1). From QIn2, the subsequent deprotonation and carboxylation proceed smoothly on the quintet potential energy surface, ultimately delivering the quintet carboxylate product QIn4. Ligand exchange between QIn4 and the alkyne substrates could regenerate the active catalyst TIn1, thereby finishing the catalytic cycle. After QIn2, the close energy of the triplet state intermediates and transition states to their quintet counterparts indicates that the deprotonation, carboxylation, and catalyst regeneration along the triplet potential energy surface could be competitive.
Throughout this catalytic cycle (an illustrative diagram is provided in Fig. 3), the rate-determining step is the spin crossover step, i.e. the conversion from the triplet bis-alkyne adduct TIn1 to the quintet state. Compared to this mechanism, all the other mechanistic possibilities, such as the prior coordination of CO2 onto the Fe(0) center (i.e. the formation of T/Q/SIn5) or the involvement of the singlet mono-alkyne coordinated Fe(0) structure SIn2, and the CO2-coordinated structures (T/Q/SIn6 and QIn7– in supporting information) can be excluded due to their relatively high energy (Fig. 2). Specifically, similar to Zhu's observation,37,38 the spin distribution of the key species involved in the most feasible pathway indicate that the substituted phenanthroline acts as a redox ligand that accepts an electron from Fe(0) to provide a formal quintet-state Fe(I) complex of the phenanthroline radical anion (Table S1).
In summary, this study reports the first Fe(0)-phenanthroline catalyzed alkyne carboxylation using CO2 as a carbon source. Steric hindrance at the 2,9-positions of the ligand appears crucial for achieving high yields (up to 95%). DFT calculations reveal a unique catalytic pathway involving spin crossover from a triplet to a quintet state, which serves as the rate-determining step. Using earth-abundant iron, this work provides a sustainable and efficient alternative to precious or/and toxic metal catalysts for synthesizing propionic acid derivatives and a spin crossover phenomenon during the catalytic process.
H. Wang: investigation, data curation, methodology, software, validation, writing – original draft; Q. Zhang: DFT calculations; K. Jacob: investigation, data curation, software, validation; A. Bousseksou: methodology, validation, resources, funding acquisition, review and editing; Philippe Hapiot: methodology, validation, review and editing; H. Yu: methodology, software, data curation, validation, resources, funding acquisition, review and editing; J.-L. Pozzo: conceptualization, methodology, resources, project administration, review and editing; D. Astruc: conceptualization, resources, project administration, review and editing.
Supplementary information (SI) is available. See DOI: https://doi.org/10.1039/d6cc00893c.
Footnote |
| † Both authors (HW and QZ) contributed equally to this work. |
| This journal is © The Royal Society of Chemistry 2026 |