Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Impact of Sprouting Process on the Protein Quality of Yellow and Red Quinoa (Chenopodium quinoa)
Next Article in Special Issue
Aluminium 8-Hydroxyquinolinate N-Oxide as a Precursor to Heterometallic Aluminium–Lanthanide Complexes
Previous Article in Journal
Assessment of Tilapia Skin Collagen for Biomedical Research Applications in Comparison with Mammalian Collagen
Previous Article in Special Issue
{GdIII7} and {GdIII14} Cluster Formation Based on a Rhodamine 6G Ligand with a Magnetocaloric Effect
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Photocatalytic Reduction of CO2 into CO with Cyclometalated Pt(II) Complexes of N^C^N Pincer Dipyridylbenzene Ligands: A DFT Study

by
Antonia Sarantou
and
Athanassios Tsipis
*
Laboratory of Inorganic Chemistry, Department of Chemistry, University of Ioannina, 45110 Ioannina, Greece
*
Author to whom correspondence should be addressed.
Molecules 2024, 29(2), 403; https://doi.org/10.3390/molecules29020403
Submission received: 23 December 2023 / Revised: 9 January 2024 / Accepted: 12 January 2024 / Published: 14 January 2024
(This article belongs to the Special Issue Exclusive Feature Papers in Inorganic Chemistry, 2nd Edition)

Abstract

:
In this work, density functional theory (DFT) calculations were employed to study the photocatalytic reduction of CO2 into CO using a series of Pt(II) square planar complexes with the general formula [Pt(5-R-dpb)Cl] (dpb = 1,3-di(2-pyridyl)benzene anion, R = H, N,N-dimethylaniline,T thiophene, diazaborinine). The CO2-into-CO conversion process is thought to proceed via two main steps, namely the photocatalytic/reduction step and the main catalytic step. The simulated absorption spectra exhibit strong bands in the range 280–460 nm of the UV-Vis region. Reductive quenching of the T1 state of the complexes under study is expected to be favorable since the calculated excited state redox potentials for the reaction with sacrificial electron donors are highly positive. The redox potentials reveal that the reductive quenching of the T1 state, important to the overall process, could be modulated by suitable changes in the N^C^N pincer ligands. The CO2 fixation and activation by the three coordinated Pt(II) catalytically active species are predicted to be favorable, with the Pt–CO2 bond dissociation energies D0 in the range of −36.9–−10.3 kcal/mol. The nature of the Pt–CO2 bond of the Pt(II) square planar intermediates is complex, with covalent, hyperconjugative and H-bonding interactions prevailing over the repulsive electrostatic interactions. The main catalytic cycle is estimated to be a favorable exergonic process.

1. Introduction

The problem of global overheating dates back long before the industrial revolution. The excessive use of fossil fuels and the massive release of greenhouse gases into the atmosphere have turned the term “global warming” into an uncomfortable reality of modern life. CO2, although not a primary greenhouse gas, nevertheless has a significant contribution to global warming. In 2012, global CO2 emissions, which come mainly from the burning of fossil fuels, reached a total of 34.5 billion tons [1].
The steadily increasing concentration of CO2 in the atmosphere is now a real problem. The CO2 emissions from the massive consumption of fossil fuels are largely responsible for global climate change: an increase in the concentration of CO2 in the atmosphere causes global warming, as well as ocean acidification from the uptake of atmospheric CO2 [2].
Man’s growing needs for energy lead to the ominous conclusion that CO2 emissions will not remain stagnant but will increase even more. This results in the scientific community focusing on how the abundance of CO2 will be used to produce materials of commercial interest. Only 1‰ of the total abundance of CO2 is currently used for chemical synthesis, which is mainly due to its chemical inertness (it is highly oxidized in nature) but also to the fact that capturing and storing CO2 is an expensive process [3].
Many studies have focused on how to deal with the problem of the continuous increase in atmospheric CO2 while, at the same time, there are ongoing efforts and studies not only to capture CO2 but also to catalytically convert it into “fine chemicals”, as well as CO, CH3OH and HCOOH [3]. A recent study by Audisio et al. demonstrated that CO, produced via a very fast photocatalytic CO2-into-CO conversion (<10 min), could be valorized to produce radiopharmaceuticals [4].
Among various strategies to achieve this goal, the greatest and most attractive challenge is the efficient conversion of CO2 into useful compounds using sunlight as an energy source [5]. Studies on this issue have shown the photochemical and electrochemical reduction of CO2 into CO, into CH3OH and also into HCOOH, using transition metal electrodes, metal complexes, semiconductors, as well as some organic molecules [6,7,8,9].
Much of the research has focused on how to sequester CO2 from coal plant emissions, which account for 76% of all global emissions [10]. Recent advances in catalyst design combined with the reduced cost of clean energy has made the use of CO2 more attractive as a feedstock for chemical production. Particular interest has gathered in the conversion of CO2 into CO as a step in the petrochemical production process [11].
Mainly, there are two approaches to CO2 reduction, namely either via heterogeneous or homogeneous catalysis. Both approaches can be achieved either electrocatalytically or photocatalytically. Heterogeneous electrocatalytic CO2 reduction is carried out using electrolysis with the aid, for example, of nanostructured materials [12], while heterogeneous photocatalytic CO2 reduction mainly utilizes various types of semiconductors as photocatalysts [13]. On the other hand, homogeneous electrochemical/photochemical CO2 reduction is mainly based on transition metal complexes [14]. So far, coordination complexes of numerous metals have been studied as catalysts, such as the first row Mn, Fe, Co, Ni and Cu, also the second row Ru, Rh, Pd, as well as the third row W, Re, Os, Ir [14].
Since 1970, a great number of works have appeared in the literature concerning electrochemical/photochemical CO2 reduction using aza macrocyclic and polypyridine complexes [15]. Along these lines, complexes of the earth-abundant Fe, Co and Ni transition metals with porphyrin-like ligands have attracted much interest. Thus, for example, a porphyrin-like aminopyridyl macrocyclic Co complex was found to be highly selective (98%) for CO production [16,17]. Fe tetraphenyl porphyrins have been demonstrated to be very active for the conversion of CO2 into CO [15]. Water-soluble Fe-porphyrins showed very high activity (close to 100%) for CO production from CO2 during electrolysis under low overpotential [18]. Phthalocyanine complexes of another earth-abundant, first-row Group VIII metal, namely Ni, were shown to very efficient for CO production, with nearly 100%, just like similar Co complexes [19]. In addition, other aza macrocyclic and polypyridine complexes have shown great potential for CO2 conversion. Lioret-Fillol et al., in a combined experimental/theoretical work, studied extensively the mechanism of CO2 reduction using a very promising Co aminopyridine catalyst [20]. Iron and cobalt quaterpyridine complexes proved to be highly efficient and selective for photocatalytic CO2 reduction [21]. Complexes of Ni and Co with the aza-crown ether cyclam ligand have been extensively studied for CO2 conversion [15]. In a pioneering study, Tinnemans et al. employed Ni as well as Co cyclams to photocatalytically reduce CO2 [22]. Since then, numerous studies on cyclam complexes have appeared. A simple Ni cyclam complex exhibited high catalytic activity for CO2 reduction, and this was attributed to the complex stabilization due to the presence of a cyclam ring, as well as to an acidic N-H proton [23,24,25]. Finally, studies have demonstrated that complexes of Mn with electron-rich ligands such as bipyridine are effective in CO2 reduction [26].
Up to now, Ru and Re complexes have received the major focus and are the most widely studied. Probably, the interest in these complexes stems from the inspiring seminal works of Lehn et al. [27,28,29], who studied electro/photocatalytic CO2 reduction using Re(I) and Ru(II) complexes. Since then, many other Re(I) and Ru(II) complexes have been studied for electro/photocatalytic CO2 reduction [2,30,31].
The reduction of CO2 into either CO or HCOOH proceeds via the following chemical equations:
CO2 + 2e + 2H+ → CO + H2O
CO2 + 2e + 2H+ → HCOOH
The reduction potentials of the above 2e/2H+ reductions are −0.52 and −0.61 V vs. SHE [31] and the e- is provided either by the electrodes in the case of electrocatalysis or by a “sacrificial donor”, e.g., TEOA in the case of photocatalysis. At least for Re(I) and Ru(II) complexes, the proposed mechanism involves, in almost all cases, the formation of a coordinatively unsaturated, 17e, very reactive intermediate [2]. Thus, for example, in the reduction of CO2 using cis-[Ru(bpy)2(CO)X]n+, a mixture of CO/HCOO is obtained [32]. Two mechanisms are proposed to operate in conjunction, having in common the very reactive, five-coordinated cis-[Ru(bpy)2(CO)]0 intermediate. The product selectivity depends upon the reactivity of this intermediate towards CO2 or H+ [2]. The so-called η1–CO2 complex mechanism operates upon CO2 capture through the cis-[Ru(bpy)2(CO)]0 intermediate, forming the cis-[Ru(bpy)2(CO)(η1-CO2)] adduct, which, upon protonations, yields CO. On the other hand, in the so-called hydride mechanism, operating simultaneously with the η1–CO2 complex mechanism, the cis-[Ru(bpy)2(CO)]0 intermediate captures H+, forming a hydrido intermediate, which upon reduction and CO2 insertion into the Ru–H bond, yields HCOO. However, it should be noticed that the unsaturated five-coordinated [Ru(bpy)2(CO)]0 intermediate has not yet been detected [2]. Accordingly, there have been proposed two other mechanisms for the reduction of CO2 [29].
Ishitani et al. [33], using various spectroscopic techniques, demonstrated the existence of the fac-[ReI(bpy)(CO)3(OCH2CH2NR2)] formed upon the coordination of the TEOA anion to the metal center of the 17e five-coordinated fac-[ReI(bpy)(CO)3]+ intermediate. They also demonstrated the existence of the fac-[ReI(bpy)(CO)3(O(CO)OCH2CH2NR2)] produced upon the insertion of CO2 into the Re–O bond of the fac-[ReI(bpy)(CO)3(OCH2CH2NR2)] intermediate. A similar CO2 insertion product has been observed also for the [Ru(dmb)2(CO)2]2+ complex in DFM-TEOA solution [34]. The CO2-into-CO reduction upon the coordination of the TEOA anion and subsequent CO2 insertion has also been scrutinized for Re(I) complexes with the aid of DFT calculations [35].
Finally, based upon both experimental [36,37] as well as theoretical [38] studies, there is also another proposed mechanism for CO2 reduction using transition metal complexes, thought to operate via the dimerization of the five-coordinated, very reactive intermediate. Hence, Muckerman et al. [38], employing DFT electronic structure calculations, showed that CO2 is inserted between two five-coordinated [Re(dmb)(CO)3]· (dmb = 4,4′-dimethyl-bpy) radical intermediates, bridging their metal centers. Through a second CO2 addition to this dimeric species, and via a transition state, CO is produced.
Although, so far, complexes of a large number of transition metals have been studied for electro/photocatalytic CO2 reduction, Pt complexes have attracted very little attention. Ceballos et al. [39] demonstrated electrocatalytic CO2 reduction into formate using the [Pt(dmpe)2](PF6)2 complex with high Faradaic efficiency and low overpotential. The vast majority of studies on CO2 reduction involving Pt metals have been devoted to heterogeneous catalysis. For example, Zhang et al. [40] dealt with the photocatalytic reduction of CO2 into CH4 under the influence of ultraviolet radiation, with water as solvent on a TiO2 catalyst with Pt. Quite numerous other studies have also appeared in the literature [41,42,43,44,45,46,47] concerning heterogeneous CO2 transformation. Finally, an enormous amount of work has been devoted to the ubiquitous photocatalytic water splitting for H2 production using Pt nanoparticles [48,49].
To the best of our knowledge, neither experimental nor theoretical studies have been undertaken so far on photocatalytic CO2-into-CO reduction employing Pt(II) square planar complexes as the catalysts. Accordingly, we instigated scrutinizing, by means of electronic structure calculations, the photocatalytic reduction of CO2 into CO using square planar Pt(II) cyclometalated complexes bearing N^C^N-coordinating substituted dipyridylbenzene ligands of the general formula [Pt(5-R-dpb)Cl] (dpb = 1,3-di(2-pyridyl)benzene anion, R = H, N,N-dimethylaniline, thiophene, bodipy). Our targets are (1) to study the photophysical properties of these complexes in order to investigate whether they could be used as both photosensitizers as well as catalysts, (2) to explore the mechanism of the catalytic reaction and (3) to probe how the nature of the N^C^N-substituted bipyridyl ligand affects both the photophysical properties and mechanism of action of the Pt(II) complexes under study.

2. Results and Discussion

The molecular structure of the square planar Pt(II) complexes employed for the CO2-into-CO photocatalytic reduction in this study is depicted in Scheme 1.
The CO2-into-CO photocatalytic reduction is thought to proceed via two steps: (1) the photoexcitation/reductive quenching step and (2) the actual catalytic cycle. The two steps are depicted in Figure 1 along with the molecular structures of all the species participating in the overall process. Let us start by examining first the initial step, which henceforth we shall call the photoexcitation/reductive quenching step, and next we will deal with the second step, which we shall call the catalytic cycle step.

2.1. Photoexcitation/Reductive Quenching Step

In the first step, the Pt(II) catalyst, upon irradiation, is excited to its first singlet excited state, S1. For a third-row transition metal complex like those under study, it is expected for the transition to occur from the S1 state to the first triplet excited state T1 via Intersystem Crossing (ISC) due to Spin–Orbit Coupling (SOC). Then, if the lifetime τ of the T1 state of the complex is large enough, reductive quenching occurs, where the complex, upon accepting an e from a donor D, e.g., TEOA or TEA, yields the One-Electron-Reduced (OER) complex (Figure 1). Notice that the experimentally determined τ of the T1 state of 1 has been found to be 7.2 μs in a CH2Cl2 solvent at 298 °C [50], which is large enough to permit reductive quenching. The optimized geometries of 14 in their S0 and T1 states, as well the geometries of the respective OER species, along with the selected structural parameters calculated at the PBE0-GD3BJ/Def2-TZVP level, in the DCM solvent are given in Figure 2, as well as in Figures S1–S3 of the Supplementary Materials. The calculated structural parameters of 1 in its S0 ground state are in excellent agreement with the respective structural parameters derived from the X-ray structural analysis of this complex [51]. Thus, for example, the difference between the calculated and X-ray derived bond lengths around the coordination sphere is in the range of only 0.001–0.016 Å, while the calculated bond angles differ by only 0.3–0.5° from the respective X-ray experimental values. Upon excitation to the T1 excited state, no significant differences in the structural parameters could be observed (Figure 2b and Figures S1–S3). Accordingly, the changes in the Pt–Cl and Pt–C coordination bonds are in the ranges 0.03–0.001 Å and 0.009–0.019 Å, respectively, while the changes in the Pt–N coordination bond lie within the range 0.0–0.015 Å. On the other hand, the bond angles around the coordination sphere change slightly by about 1° in all cases, 14. Finally, the same holds also for the organic framework of the N^C^N pincer ligand, where the structural changes upon S0 → T1 excitation are negligible. Next, the one-electron reduction of the T1 state causes more significant structural changes. Thus, there are more obvious structural differences between the optimized geometries of the T1 state and the OER species (Figure 2c and Figures S1–S3). The most striking structural change is observed in the Pt–Cl coordination bond, which is elongated by 0.007–0.032 Å in the OER species as compared to the T1 state species. The same holds also if we compare the OER species with the S0 state species. Also, the Pt–C coordination bond is elongated by 0.035 Å in contrast to the two Pt–N coordination bonds, which are shortened.
Finally, the bond angles around the coordination sphere do not change significantly upon one-electron reduction.

2.1.1. Absorption Spectra

Since the excitation of the precursor complexes 14 is a prerequisite in photocatalytic CO2 conversion, we set out to study their absorption spectra in DCM by means of TDDFT calculations. Accordingly, we employed the two most popular functionals (vide infra), for the simulation of the absorption spectra in DCM of the complexes under study, namely PBE0 and CAM-B3LYP, using the Def2-TZVP basis set in both cases. The simulated absorption spectra of 14 calculated using the TDDFT/PBE0/Def2-TZVP computational protocol are depicted in Figure 3, while those obtained using the TDDFT/CAM-B3LYP/Def2-TZVP computational protocol are given in Figure S4 of the Supplementary Materials.
The simulated absorption spectrum of 1 calculated at the PBE0/Def2-TZVP level in DCM is in agreement with the experimental absorption spectrum derived by Williams et al. for 1 in DCM [50]. Thus, the simulated spectrum of 1 exhibits a band at 380 nm, reproducing excellently the experimental maximum found at 380 nm in the experimental spectrum as well. On the other hand, the simulated absorption spectrum of 1 at the CAM-B3LYP/Def2-TZVP level (Figure S4) shows a maximum at 344 nm, underestimating the experimental maximum by about 36 nm. The simulated absorption spectra using both computational protocols employed in this study exhibit also higher energy maxima below 300 nm, underestimating the respective maxima found in the experimental spectra at 330 nm and below. It should be noticed that overall, the simulated absorption spectra of 14 at the CAM-B3LYP/Def2-TZVP level are blue-shifted by 30–40 nm compared to those calculated at the PBE0/Def2-TZVP level in DCM. Let us further analyze the spectra calculated at the PBE0/Def2-TZVP level due to their better agreement with the experiment. Perusal of Figure 3 reveals that the simulated absorption spectra at the PBE0/Def2-TZVP level in DCM of 1, 2, 3 and 4 exhibit low-energy bands at 380, 460, 420 and 428 nm, respectively. In addition, the simulated absorption spectra of all complexes exhibit high-energy bands in the region of 250–300 nm.
In terms of Natural Transition Orbitals (NTOs), the band peaking at 380 nm arises mainly from an electronic transition at 381 nm, which is due to the electron excitation between the NTO pair depicted in Figure 4. Based upon the shape of these NTOs, the band shows a mixed assignment comprising of Metal-to-Ligand Charge Transfer (MLCT), Intra-Ligand (IL) and Ligand-to-Ligand Charge Transfer (LL’CT) characters (MLCT/IL/L L′CT). The same assignment holds also true for the low-energy bands in the region 420.
The calculated absorption spectra of 2 and 3 exhibit a band around 460 nm. The former exhibits a band around 460 nm while the latter exhibits a band around 420 nm, arising from electronic transitions at 452 and 414 nm, respectively. Based upon the shapes of the NTO pairs (Figures S5 and S6 of Supplementary Materials) participating in these electronic transitions, these low-energy bands are assigned as MLCT/IL/LL′CT. The exception is complex 4, where the band around 428 nm arises mainly from an electronic transition at 420 nm, assigned as MLCT/IL.
Next, the simulated absorption spectra exhibit high-energy bands in the UV region (Figure 3). The spectrum of complex 1 shows two such high-energy bands around 250 and 280 nm. The latter arises mainly from an electronic transition at 288 nm while the former is mainly due to an electronic transition at 249 nm. The shape of the NTOs involved in these transitions indicates that the high-energy bands in the UV are of MLCT/IL/LL′CT character. On the other hand, the absorption spectra of complexes 24 exhibit just one distinct high-energy band in the UV region. This high-energy band appears around 295 nm for 2 and 3 and around 284 nm for 4. For complexes 2 and 3, these bands are also of the highest intensity in their absorption spectra. For 2, the most intense high-energy band arises mainly from an electronic transition at 298 nm, involving almost solely an excitation between a pair of NTOs (Figure S5). Based upon the shape of these NTOs, the band at 285 nm for 2 is assigned as MLCT/IL. Conversely, for 3, the high-energy band arises mainly from an electronic transition at 279 nm, which is, however, due to excitations involving two NTO pairs (Figure S6). The band around 295 nm for 3 is assigned as MLCT/IL/LL′CT in accordance with the shape of the NTOs involved.

2.1.2. Excited State Electrochemistry

One of the most important steps in photocatalytic CO2 conversion using transition metal catalysts is the reductive quenching of their T1 state upon receiving an e from a sacrificial donor, e.g., TEOA. Therefore, we instigated studying the T1 excited state reduction potentials using DFT calculations. The ground-state redox potentials could be calculated by employing the Born–Haber cycle [52] depicted schematically in Scheme 2.
The standard absolute ground-state reduction potential E0red is calculated using the following equation:
E0red = −ΔG°(soln. redox)/ZF,
where F is the Faraday constant (23.061 kcal per volt gram equivalent) and Z is the unity for the one-electron redox processes. ΔG°(soln., redox) is obtained from the following equation:
ΔG°(soln., redox) = ΔG°(gas, redox) + ΔG°(solv.,[Pt(5-R-dpb)Cl]•−) − ΔG°(solv.,[Pt(5-R-dpb)Cl], S0)
In Table 1 are tabulated the values of the reduction potentials calculated for complexes 14 calculated at the PBE0/def2-TZVP level. In Table 1 are given also the reduction potentials of the sacrificial electron donors D TEOA and TEA at the same level of theory.
Inspection of Table 1 reveals that the oxidizing ability of the Pt(II) complexes in their T1 state follows the series 2 > 1 > 3 > 4, i.e., 2 is anticipated to accept an e in the T1 state and be the most easily reduced amongst all the complexes under study. In contrast, the oxidizing ability of the Pt(II) complexes in their S0 ground state follows exactly the same order as that observed for the T1 state, following the order 2 > 1 > 3 > 4. We could draw the same conclusion based also upon the values of the overall redox potential Ε0 calculated for the redox reaction of 14 in the T1 state with TEOA or TEA, which is favorable since it is positive in all cases. The most favorable redox reaction will be with 2, followed by 1, 3 and 4.

2.2. Catalytic Cycle Step

The initial photoexcitation/reductive quenching step ends with the formation of the OER species. The latter, upon losing a chloride ligand, yields a three-coordinated intermediate ImA, which starts the actual catalytic cycle (Figure 1). ImA could be considered the “true” catalyst while the initial Pt(II) complex could be considered a “precatalyst”. Based upon earlier studies [2], the catalytic cycle step could proceed upon the one-electron reduction of ImA with the aid of an electron donor D like, for example, TEOA, yielding ImB. Next, there is one of the most important steps, i.e., the CO2 addition to the metal center of ImB. After two subsequent protonations (e.g., with TEOA, which can act as proton source as well [2]), we arrive, through the transition state TS, to the carbonyl intermediate ImF. A second electron reduction of the latter yields ImG, from which we obtain CO and the initial “true” catalyst ImA.

2.2.1. The η1–CO2 Complex

Since the CO2 capture/activation step, upon coordination to the metal center of the catalyst, is considered to be of paramount importance in electro/photocatalytic CO2 conversion [53,54], we set out to study in depth the structural, bonding and electronic properties of ImC.

Structural Properties

The optimized geometry of 1_ImC calculated at the PBE0-D3/Def2-TZVP level, in DCM solvent, is depicted in Figure 5a, while those of 2_ImC, 3_ImC and 4_ImC are given in Figures S5–S7 of the Supplementary Material.
Perusal of Figure 5a and Figures S5–S7 reveals that, in all complexes, CO2 is coordinated to the Pt metal center via the C atom in an η1-bonding mode. The complex formed via the η1–CO2 bonding mode is considered to be of pivotal importance to the majority of the electron transfer reactions for photo/electrocatalytic CO2 conversion [2,53]. The Pt–CO2 bond length is estimated to be 2.132–2.135 Å, indicative of a bonding interaction. Upon coordination to the Pt metal center, the CO2 molecule is activated and changes from linear to bended, with a <O–C–O bond angle equal to 120.8°, while the C–O bonds are equal to 1.270 Å, elongated by about 0.110 Å as compared to the “free” CO2 molecule (1.160 Å).

Bonding and Electronic Properties

The covalent nature of the Pt–CO2 bond is reflected in the shape of HOMO calculated for 1_ImC and depicted schematically in Figure 5b (the HOMOs of 2_ImC, 3_ImC and 4_ImC). The HOMO is a bonding MO mainly constructed by the in-phase combination of the Pt dz2 AO with the π* MO of CO2, in line with previous studies [51]. The bond dissociation energy D0 of the Pt–CO2 bond was found to be −36.8, −36.9, −35.1 and −10.3 kcal/mol for 1_ImC, 2_ImC, 3_ImC and 4_ImC, respectively. Obviously, substitution into position 5 of the N^C^N pincer ligand has no significant impact on D0(Pt–CO2). The exception is complex 4, where the introduction of the diazaborinine group reduces D0(Pt–CO2) by a third as compared to the rest of the complexes. Nevertheless, based upon the magnitude of D0(Pt-CO2), the interaction between ImB and CO2 is expected to be a relatively favorable interaction.
To further analyze the Pt–CO2 bond in ImC, we employed the atoms in molecules (AIM) method. In Figure 5c are depicted the bond-critical points (BCPs) found using the AIM method for 1_ImC (similar BCPs are detected also for the intermediate ImC formed by complexes 24). According to Bader’s [55,56] theory, the presence of a BCP between two atoms indicates bond formation. Inspection of Figure 5c reveals the existence of a BCP between the Pt metal center and C atom of the CO2 ligand, indicating a Pt–CO2 bonding interaction. It has been proposed [57,58], based upon the values of certain properties of BCPs, that the bonding interactions could be classified into three categories: (a) pure closed-shell interactions (e.g., ionic bonds, hydrogen bonds and van der Waals interactions) characterized by |VBCP|/GBCP < 1 (∇2ρBCP > 0 and HBCP > 0); (b) pure open-shell (covalent) interactions characterized by |VBCP|/GBCP > 2 (∇2ρBCP < 0 and HBCP < 0) and (c) intermediate bonds with 1 < |VBCP|/GBCP < 2 (i.e., ∇2ρBCP > 0 and HBCP < 0) where VBCP is the potential energy density at the BCP, GBCP is the kinetic energy density at the BCP, ∇2ρBCP is the Laplacian of electron density ρ at the BCP and finally HBCP is the energy density at the BCP. Accordingly, the calculated values of these properties for the BCP found between Pt and CO2 in intermediates 1_ImC4_ImC are given in Table 2.
Inspection of Table 1 reveals that the AIM parameters calculated for the Pt–CO2 BCP do not vary significantly between the ImC intermediates of 14. Based upon the values of AIM parameters given in Table 1, we can assume that Pt–CO2 interaction falls into the third category, i.e., intermediate nature arising from an interplay of covalent, electrostatic, charge transfer and probably weak dispersion force components. However, the latter should be excluded for the Pt–CO2 bond, in terms of the Reduced Density Gradient (RDG) defined using the following equation [59]:
RDG(r) = |∇ρ(r)|/2(3π2)1/3ρ(r)4/3
where ρ(r) is the electron density ρ at point r. The 3D isosurface map of the RDG, depicted in Figure 5d, reveals the absence of any weak interactions between the Pt metal center and CO2 (notice that the green isosurface represents weak interactions regions). Similar RDG maps are observed for complexes 24 as well. In addition, inspection of Figure 5c reveals the existence of two BCPs located between the O atoms of CO2 and the H atoms of the N^C^N pincer ligand. Also, the RDG function (Figure 5d) reveals non-covalent interactions between those atoms. Therefore, these H bond interactions further contribute to the overall interaction of CO2 with the Im_B intermediate.
Another useful method for analyzing the chemical bond is the NBO method, which is a partitioning scheme for the electronic charge distribution in a molecule. Thus, we employed NBO calculations to analyze the Pt–CO2 bond from this point of view. The results of the NBO analysis are given in Table 3.
The natural charges obtained using NBO analysis are found to be 0.235–0.260 for Pt and 0.494–0.497 for the C atom of the CO2 ligand. This results in a repulsive electrostatic interaction between the Pt and CO2 in line with previous findings [53,60]. The NBO analysis reveals also the existence of a bonding BD σ NBO, which is depicted in Figure 6a for 1_ImC (similar BD NBOs are observed also for the rest of the respective intermediates).
The linear combinations of the BD[σ(Pt–CO2)] NBOs found for 1_ImC4_ImC are given in Table 2. Thus, for example, the σ(Pt–C) bonding NBO of 1_ImC has an occupation number equal to 1.857 |e|, arising from the interaction of the sp2.23d1.08 hybrid orbital of Pt (23.15% s, 51.66% p and 25.03% d character) with the sp1.88 hybrid orbital of the C donor atom of CO2 (34.72% s and 65.25% p character), and is described as σ(Pt–C) = 0.558hPt + 0.830hC. In Figure 6b are also depicted the shapes of the NBOs participating in the donor– acceptor hyperconjugative interactions related to the Pt–CO2 bond in 1_ImC. According to NBO analysis [61], there is a stabilization of a system due to the charge transfer interactions between specific donor-acceptor NBOs. The stabilization energy, ΔE(2), arising from these hyperconjugative (stereoelectronic) interactions is given using the following equation:
ΔE(2) = qiFij/(εiεj)
where qi is the occupancy of the donor NBO, Fij is the off-diagonal Fock matrix elements and εi and εj are the energy of the donor and acceptor NBOs, respectively. The two hyperconjugative interactions related to the Pt–CO2 bond given in Figure 6b involve σ(Pt–C) and σ*(Pt–C) NBOs and stabilize the system overall by 86.1 kcal/mol. The overall stabilization due to hyperconjugative interactions for 2_ImC, 3_ImC and 4_ImC is 84.5, 85.5 and 80.5 kcal/mol, respectively.
Another method for analyzing a chemical bond is Charge Decomposition Analysis (CDA) [62]. We employed the latter to study the Pt–CO2 interaction for 1_ImC as a representative case. In contrast to the Co–CO2 interaction in Co complexes, where charge transfer has been proposed [63,64], the net charge transfer from the Pt fragment to the CO2 is calculated to be marginal, amounting to only 0.048 |e| (the donation d and backdonation b terms were found to be −0.023 and −0.071 |e|, respectively, and are somewhat balanced). Finally, the charge polarization term r has a negative value equal to −0.350, indicating electronic charge removal from the overlapping into the non-overlapping regions upon the formation of 1_ImC. An orbital interaction diagram for the formation of 1_ImC from the fragments [Pt(dpb)] (Fragm. 1) and CO2 (Fragm. 2) is depicted in Scheme 3.
Perusal of Scheme 3 reveals that the bonding HOMO of the complex 1_ImC, reflecting the covalent Pt-CO2 interaction, is composed of 34% of LUFO+1 of Fragm. 1 and 55% of LUFO of Fragm. 2. In other words, the HOMO of 1_ImC is constructed by the in-phase interaction of the, mainly of the Pt dz2 character, LUFO+1 of the [Pt(dpb)] fragment, with a π* MO of the CO2 fragment.

2.2.2. Energetic Reaction Profiles

The free energy reaction profiles calculated for 14, corresponding to the respective catalytic cycle step (Figure 1), are given in Figure 7. The catalytic cycle starts with ImA, which is obtained from OER upon the loss of the chloride ligand (vide supra). Experimental studies have already proposed the loss of the halide ligand from the Re and Ru complexes used in catalytic CO2 conversion [2]. Our study indicates that this should also hold for the Pt complexes under investigation. Therefore, the latter, upon reduction in their T1 excited state, accept an electron in the SOMOβ depicted in Scheme 4. Putting an electron in this antibonding MO would weaken the Pt–Cl bond, thus facilitating its dissociation. This is further corroborated by the fact that the single electron in the OER species resides on the antibonding SOMO (Scheme 4), destabilizing the Pt–Cl bond and favoring its dissociation.
Next, inspection of the energetic profiles reveals that the CO2-into-CO conversion by the catalytically active species ImA is a strongly exergonic reaction. Upon reduction, ImA yields the reactive intermediate ImB, which subsequently captures a CO2 molecule to yield the η1–CO2 complex ImC. Successive protonation of the latter results in the formation of the immensely stabilized intermediate ImE. Next, ImE is converted into the carbonyl intermediate ImF through a concerted mechanism and via a transition state, TS, surmounting an activation barrier of around 38 kcal/mol. Finally, upon a second electron reduction, we obtain CO while regenerating the catalytically active species ImA.
It should be noticed that the dissociation of the Pt–CO upon the reduction of the ImF intermediate should be favorable since it accepts an electron in its antibonding LUMO (Scheme 4). Also, the single electron in the ImG intermediate resides on the antibonding SOMOβ (Scheme 4), destabilizing the Pt–CO bond. These findings are in line with experimental studies proposing CO dissociation upon the reduction of various complexes acting as CO2 reduction catalysts [2].
The optimized geometries, along with selected structural parameters, of all the intermediates and TSs involved in the catalytic cycle are given in Figures S7–S10. Perusal of the latter reveals that the substitution of one H atom of 1 by an R substituent (Scheme 1) has no significant impact on the structural parameters of the [Pt(5-R-dpb)Cl] complexes.

3. Computational Methods

Full geometry optimization has been performed for all the species under study, without symmetry constraints, using the 1997 hybrid functional of Perdew, Burke and Ernzerhof [65,66,67,68,69,70], as implemented in the program Gaussian16W [71]. This functional uses 25% exchange and 75% weighting correlation and is denoted as PBE0. Dispersion interactions were accounted for by using the D3 version of Grimme dispersion with Becke−Johnson damping [72]. The Def2-TZVP basis set for all atoms was used for the geometry optimizations. The computational protocol will be hereafter denoted as PBE0-GD3BJ/Def2-TZVP. All stationary points have been identified as minima (number of imaginary frequencies, NImag = 0). Natural bond orbital (NBO) population analysis was carried out, employing the methodology by Weinhold [63] as implemented in the G16W software, version C.01. The atoms in molecules (AIM) method of Bader [57], the reduced gradient density (RDG) method [73] and the CDA method [62] were used as implemented in the Multiwfn software, version 3.8 [74]. The Gibbs free energy was calculated to be 298.15 K and 1 atm pressure. Solvent effects were calculated using the Polarizable Continuum Model (PCM) using the integral equation formalism variant (IEF-PCM), which is the default method of G16W (self-consistent reaction field (SCRF)) [75], while Dichloromethane (DCM) was used as the solvent. Time-dependent density functional theory (TD-DFT) calculations [76,77,78] were performed on the ground-state S0 equilibrium geometries in the DCM solvent employing two computational protocols, namely PBE0/Def2-TZVP and CAM-B3LYP/Def2-TZVP, taking into account the first 30 excited states. Both PBE0 and CAM-B3LYP are among the most popular functionals for the calculation of absorption spectra [79]. Also, the def2-TZVP basis set has been successfully employed for the calculation of the absorption spectra of the Pt(II) square planar complexes [80]. Finally, Natural Transition Orbitals were used, offering a greatly simplified description of the electronic transitions. [81].

4. Conclusions

The photocatalytic reduction of CO2 into CO by [Pt(5-R-dpb)Cl] (dpb = 1,3-di(2-pyridyl)benzene anion, R = H, N,N-dimethylaniline, thiophene, bodipy) Pt(II) square planar complexes has been studied using DFT electronic structure calculations. The overall process is thought to proceed via two main steps: (i) photoexcitation/reduction of the initial complexes and (ii) the main catalytic cycle, producing CO. The TDDFT-simulated absorption spectra showed that these complexes exhibit strong absorptions in the range of 280–380 of the UV region, as well as weak absorptions in the range of 410–460 nm. The exception is complex 4, with a diazaborinine-substituted N^C^N pincer ligand, where this phenomenon is reversed, with the strongest absorption appearing around 420 nm and the weakest around 280 nm. Thus, it is anticipated that changing the pincer ligand could provide complexes that absorb in the visible region, making them ideal for use as photocatalysts for CO2-into-CO conversion. Next, the excited state reduction potentials, E0*red, dictating the reductive quenching of the Pt(II) complexes in their T1 state have been calculated using DFT. Based upon the E0*red values, we conclude that the oxidizing ability of the Pt(II) complexes in their T1 state follows the series 2 > 1 > 3 > 4, meaning that 2 is expected to be more easily reduced amongst all the complexes under study. In general, based upon the calculated redox potentials E0*redox in the range of 6–8 V, all the complexes are expected to be easily reduced with a sacrificial electron donor such as TEOA or TEA. During reduction, the complexes 14 in the T1 state accept an e in their SOMO, being antibonding in the region of the Pt–Cl bond. In the OER product species, resulting from the reduction of the T1 state complexes, the single e resides also in an antibonding SOMO, destabilizing and facilitating the Pt–Cl bond dissociation process.
The second step of photocatalytic CO2-into-CO conversion, namely the catalytic cycle, starts with a three-coordinated Pt(II) ImA intermediate, obtained from OER upon Pt–Cl dissociation. Upon one-e reduction, ImA captures CO2, and after successive protonations and reductions, yields CO. The latter is obtained from ImF, which accepts an e in its LUMO, being antibonding in the region of the Pt–CO bond, while this e resides in a similar antibonding SOMO of the reduced ImG intermediate. Consequently, the Pt–CO bond is expected to be broken in the final stage of the catalytic cycle.
The Pt–CO2 bond of the ImC intermediate, crucial to the overall process, has been scrutinized. It is revealed that this bond is relatively strong with D0 in the range of −36.9–−10.3 kcal/mol, exhibiting a complex nature comprising mainly covalent and hyperconjugative interactions, which compensate for the repulsive electrostatic interactions. The Pt–CO2 bond involves also H bonding with the N^C^N pincer ligand. Overall, the catalytic cycle is estimated to be a strongly exergonic process. Taking into account that 1 exhibits a very high T1 excited state lifetime τ, a whole new series of Pt(II) complexes could be synthesized bearing suitable pincer ligands in order to absorb in the visible range, making them ideal for CO2-into-CO conversion using sunlight.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules29020403/s1, Figure S1: Optimized geometries of 2 in S0 ground state 2_S0, in T1 excited state 2_T1 and of One-Electron-Reduced form 2_OER, calculated at the PBE0-GD3BJ/Def2-TZVP levels in DCM solvent. Figure S2: Optimized geometries of 3 in S0 ground state 3_S0, in T1 excited state 3_T1 and of One-Electron-Reduced form 3_OER, calculated at the PBE0-GD3BJ/Def2-TZVP levels in DCM solvent. Figure S3: Optimized geometries of 4 in S0 ground state 4_S0, in T1 excited state 4_T1 and of One-Electron-Reduced form 4_OER, calculated at the PBE0-GD3BJ/Def2-TZVP levels in DCM solvent. Figure S4: Simulated absorption spectra of 14 in DCM at the TDDFT/CAM-B3LYP/Def2-TZVP levels. Figure S5: NTO pairs for the most significant electronic transitions of the simulated absorption spectra of 2 calculated at TDDFT/PBE0/Def2-TZVP levels (hole h+ on the left, electron e on the right). Figure S6: NTO pairs for the most significant electronic transitions of the simulated absorption spectra of 3 calculated at TDDFT/PBE0/Def2-TZVP levels (hole h+ on the left, electron e on the right). Figure S7: NTO pairs for the most significant electronic transitions of the simulated absorption spectra of 4 calculated at TDDFT/PBE0/Def2-TZVP levels (hole h+ on the left, electron e on the right). Figure S8: Optimized geometries of all species involved in the catalytic cycle of CO2-into-CO conversion by 1 calculated at the PBE0-GD3BJ/Def2-TZVP levels in DCM solvent. Figure S9: Optimized geometries of all species involved in the catalytic cycle of CO2-into-CO conversion by 2 calculated at the PBE0-GD3BJ/Def2-TZVP levels in DCM solvent. Figure S10: Optimized geometries of all species involved in the catalytic cycle of CO2-into-CO conversion by 3 calculated at the PBE0-GD3BJ/Def2-TZVP levels in DCM solvent. Figure S11: Optimized geometries of all species involved in the catalytic cycle of CO2-into-CO conversion by 4 calculated at the PBE0-GD3BJ/Def2-TZVP levels in DCM solvent. Table S1: Cartesian coordinates and energetic data of the optimized geometries of all species calculated at the PBE0-GD3BJ/Def2-TZVP levels in DCM solvent.

Author Contributions

Conceptualization, A.T.; validation, A.T.; formal analysis, A.T. and A.S.; investigation, A.T. and A.S.; data curation, A.T. and A.S.; writing—original draft preparation, A.T. and A.S.; writing—review and editing, A.T.; visualization, A.T. and A.S.; supervision, A.T.; project administration, A.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available in article and Supplementary Materials.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Lahijani, P.; Zainal, Z.A.; Mohammadi, M.; Mohamed, A.R. Conversion of the greenhouse gas CO2 to the fuel gas CO via the Boudouard reaction: A review. Renew. Sustain. Energy Rev. 2015, 41, 615–632. [Google Scholar] [CrossRef]
  2. Kuramochi, Y.; Ishitani, O.; Ishida, H. Reaction mechanisms of catalytic photochemical CO2 reduction using Re(I) and Ru(II) complexes. Coord. Chem. Rev. 2018, 373, 333–356. [Google Scholar] [CrossRef]
  3. Cokoja, M.; Bruckmeier, C.; Rieger, B.; Herrmann, W.A.; Kühn, F.E. Transformation of carbon dioxide with homogeneous transition-metal catalysts: A molecular solution to a global challenge? Angew. Chemie Int. Ed. 2011, 50, 8510–8537. [Google Scholar] [CrossRef] [PubMed]
  4. Monticelli, S.; Talbot, A.; Gotico, P.; Caillé, F.; Loreau, O.; Del Vecchio, A.; Malandain, A.; Sallustrau, A.; Leibl, W.; Aukauloo, A.; et al. Unlocking full and fast conversion in photocatalytic carbon dioxide reduction for applications in radio-carbonylation. Nat. Comm. 2023, 14, 4451–4461. [Google Scholar] [CrossRef] [PubMed]
  5. Morris, A.J.; Meyer, G.J.; Fujita, E. Molecular approaches to the photocatalytic reduction of carbon dioxide for solar fuels. Acc. Chem. Res. 2009, 42, 1983–1994. [Google Scholar] [CrossRef] [PubMed]
  6. Benson, E.E.; Kubiak, C.P.; Sathrum, A.J.; Smieja, J.M. Electrocatalytic and homogeneous approaches to conversion of CO2 to liquid fuels. Chem. Soc. Rev. 2009, 38, 89–99. [Google Scholar] [CrossRef] [PubMed]
  7. Takeda, H.; Koike, K.; Morimoto, T.; Inumaru, H.; Ishitani, O. Photochemistry and photocatalysis of rhenium(I) diimine complexes. Adv. Inorg. Chem. 2011, 63, 137–186. [Google Scholar]
  8. Costentin, C.; Robert, M.; Savéant, J.M. Catalysis of the electrochemical reduction of carbon dioxide. Chem. Soc. Rev. 2013, 42, 2423–2436. [Google Scholar] [CrossRef]
  9. Das, S.; Wan Daud, W.M.A. Photocatalytic CO2 transformation into fuel: A review on advances in photocatalyst and photoreactor. Renew. Sustain. Energy Rev. 2014, 39, 765–805. [Google Scholar] [CrossRef]
  10. Vandezande, J.E.; Schaefer, H.F. CO2 Reduction Pathways on MnBr(N-C)(CO)3 Electrocatalysts. Organometallics 2018, 37, 337–342. [Google Scholar] [CrossRef]
  11. Paquin, F.; Rivnay, J.; Salleo, A.; Stingelin, N.; Silva, C. Multi-phase semicrystalline microstructures drive exciton dissociation in neat plastic semiconductors. J. Mater. Chem. C 2015, 3, 10715–10722. [Google Scholar] [CrossRef]
  12. Yang, C.; Wang, Y.; Qian, L.; Al-Enizi, A.M.; Zhang, L.; Zheng, G. Heterogeneous Electrocatalysts for CO2 Reduction. Appl. Energy Mater. 2021, 4, 1034–1044. [Google Scholar] [CrossRef]
  13. Li, K.; Peng, B.; Peng, T. Recent Advances in Heterogeneous Photocatalytic CO2 Conversion to Solar Fuels. ACS Catal. 2016, 6, 7485–7527. [Google Scholar] [CrossRef]
  14. Ishida, H. Electrochemical/Photochemical CO2 Reduction Catalyzed by Transition Metal Complexes. In Carbon Dioxide Chemistry, Capture and Oil Recovery; Iyad, K., Shaya, J., Srour, H., Eds.; InTechOpen Ltd.: London, UK, 2018; pp. 17–40. [Google Scholar]
  15. Boutin, E.; Merakeb, L.; Ma, B.; Boudy, B.; Wang, M.; Bonin, J.; Anxolabèhère-Mallart, E.; Robert, M. Molecular catalysis of CO2 reduction: Recent advances and perspectives in electrochemical and light-driven processes with selected Fe, Ni and Co aza macrocyclic and polypyridine complexes. Chem. Soc. Rev. 2020, 49, 5772–5809. [Google Scholar] [CrossRef] [PubMed]
  16. Chapovetsky, A.; Do, T.H.; Haiges, R.; Takase, M.K.; Marinescu, S.C. Proton-Assisted Reduction of CO2 by Cobalt Aminopyridine Macrocycles. J. Am. Chem. Soc. 2016, 138, 5765–5768. [Google Scholar] [CrossRef] [PubMed]
  17. Chapovetsky, A.; Welborn, M.; Luna, J.M.; Haiges, R.; Miller, T.F.; Marinescu, S.C. Pendant Hydrogen-Bond Donors in Cobalt Catalysts Independently Enhance CO2 Reduction. ACS Cent. Sci. 2018, 4, 397–404. [Google Scholar]
  18. Costentin, C.; Robert, M.; Saveant, J.M.; Tatin, A. Efficient and selective molecular catalyst for the CO2-to-CO electrochemical conversion in water. Proc. Natl. Acad. Sci. USA 2015, 112, 6882–6886. [Google Scholar] [CrossRef] [PubMed]
  19. Furuya, N.; Matsui, K. Electroreduction of carbon dioxide on gas-diffusion electrodes modified by metal phthalocyanines. J. Electroanal. Chem. 1989, 271, 181–190. [Google Scholar] [CrossRef]
  20. Fernández, S.; Franco, F.; Casadevall, C.; Martin-Diaconescu, V.; Luis, J.M.; Lloret-Fillol, J. A Unified Electro- and Photocatalytic CO2 to CO Reduction Mechanism with Aminopyridine Cobalt Complexes. J. Am. Chem. Soc. 2020, 142, 120–133. [Google Scholar] [CrossRef]
  21. Guo, Z.; Cheng, S.; Cometto, C.; Anxolabéhère-Mallart, E.; Ng, S.-M.; Ko, C.C.; Liu, G.; Chen, L.; Robert, M.; Lau, T.-C. Highly Efficient and Selective Photocatalytic CO2 Reduction by Iron and Cobalt Quaterpyridine Complexes. J. Am. Chem. Soc. 2016, 138, 9413–9416. [Google Scholar] [CrossRef]
  22. Tinnemans, A.H.A.; Koster, T.P.M.; Thewissen, D.H.M.W.; Mackor, A. Tetraaza-macrocyclic cobalt(II) and nickel(II) complexes as electron-transfer agents in the photo(electro)chemical and electrochemical reduction of carbon dioxide. Recl. Trav. Chim. Pays-Bas 1984, 103, 288–295. [Google Scholar] [CrossRef]
  23. Beley, M.; Collin, J.-P.; Ruppert, R.; Sauvage, J.-P. Nickel (II)-cyclam: An extremely selective electrocatalyst for reduction of CO2 in water. J. Chem. Soc. Chem. Commun. 1984, 1315–1316. [Google Scholar] [CrossRef]
  24. Beley, M.; Collin, J.-P.; Ruppert, R.; Sauvage, J.-P. Electrocatalytic reduction of carbon dioxide by nickel cyclam2+ in water: Study of the factors affecting the efficiency and the selectivity of the process. J. Am. Chem. Soc. 1986, 108, 7461–7467. [Google Scholar] [CrossRef]
  25. Collin, J.-P.; Jouaiti, A.; Sauvage, J.-P. Electrocatalytic properties of (tetraazacyclotetradecane) nickel (2+) and Ni2(biscyclam)4+ with respect to carbon dioxide and water reduction. Inorg. Chem. 1988, 27, 1986–1990. [Google Scholar] [CrossRef]
  26. Nandal, N.; Jain, S.L. A review on progress and perspective of molecular catalysis in photoelectrochemical reduction of CO2. Coord. Chem. Rev. 2022, 451, 214271. [Google Scholar] [CrossRef]
  27. Lehn, J.-M.; Ziessel, R. Photochemical generation of carbon monoxide and hydrogen by reduction of carbon dioxide and water under visible light irradiation. Proc. Natl. Acad. Sci. USA 1982, 79, 701–704. [Google Scholar] [CrossRef]
  28. Hawecker, J.; Lehn, J.-M.; Ziessel, R. Efficient photochemical reduction of CO2 to CO by visible light irradiation of systems containing Re(bipy)(CO)3X or Ru(bipy)32+–Co2+ combinations as homogeneous catalysts. Chem. Commun. 1983, 536–538. [Google Scholar] [CrossRef]
  29. Hawecker, J.; Lehn, J.-M.; Ziessel, R. Photochemical and Electrochemical Reduction of Carbon Dioxide to Carbon Monoxide Mediated by (2,2′-Bipyridine)tricarbonylchlororhenium(I) and Related Complexes as Homogeneous Catalysts. Helv. Chim. Acta 1986, 69, 1990–2012. [Google Scholar] [CrossRef]
  30. Elgrishi, N.; Chambers, M.B.; Wang, X.; Fontecave, M. Molecular polypyridine-based metal complexes as catalysts for the reduction of CO2. Chem. Soc. Rev. 2017, 46, 761–796. [Google Scholar] [CrossRef]
  31. Yamazaki, Y.; Takeda, H.; Ishitani, O. Photocatalytic reduction of CO2 using metal complexes. J. Photochem. Photobiol. C Photochem. Rev. 2015, 25, 106–137. [Google Scholar] [CrossRef]
  32. Voyame, P.; Toghill, K.E.; Méndez, M.A.; Girault, H.H. Photoreduction of CO2 Using [Ru(bpy)2(CO)L]n+ Catalysts in Biphasic Solution/Supercritical CO2 Systems. Inorg. Chem. 2013, 52, 10949–10957. [Google Scholar] [CrossRef]
  33. Morimoto, T.; Nakajima, T.; Sawa, S.; Nakanishi, R.; Imori, D.; Ishitani, O. CO2 Capture by a Rhenium(I) Complex with the Aid of Triethanolamine. J. Am. Chem. Soc. 2013, 135, 16825–16828. [Google Scholar] [CrossRef] [PubMed]
  34. Tamaki, Y.; Morimoto, T.; Koike, K.; Ishitani, O. Photocatalytic CO2 reduction with high turnover frequency and selectivity of formic acid formation using Ru(II) multinuclear complexes. Proc. Natl. Acad. Sci. USA 2012, 109, 15673–15678. [Google Scholar] [CrossRef] [PubMed]
  35. Tsipis, A.C.; Sarantou, A.A. DFT insights into the photocatalytic reduction of CO2 to CO by Re(I) complexes: The crucial role of the triethanolamine “magic” sacrificial electron donor. Dalton Trans. 2021, 50, 14797–14809. [Google Scholar] [CrossRef] [PubMed]
  36. Hayashi, Y.; Kita, S.; Brunschwig, B.S.; Fujita, E. Involvement of a Binuclear Species with the Re-C(O)O-Re Moiety in CO2Reduction Catalyzed by Tricarbonyl Rhenium(I)Complexes with Diimine Ligands: Strikingly Slow Formation of the Re-Re and Re-C(O)O-Re Species from Re(dmb)(CO)3S (dmb = 4,4’-Dimethyl-2,2’-bipyridine, S = Solvent). J. Am. Chem. Soc. 2003, 125, 11976–11987. [Google Scholar] [PubMed]
  37. Sullivan, B.P.; Bolinger, C.M.; Conrad, D.; Vining, W.J.; Meyer, T.J. One- and Two-electron Pathways in the Electrocatalytic Reduction of CO2 by fac-Re(bpy)(CO)3Cl (bpy = 2,2’-bipyridine). J. Chem. Soc. Chem. Commun. 1985, 20, 1414–1416. [Google Scholar] [CrossRef]
  38. Agarwal, J.; Fujita, E.; Schaefer, H.F., III; Muckerman, J.T. Mechanisms for CO Production from CO2 Using Reduced Rhenium Tricarbonyl Catalysts. J. Am. Chem. Soc. 2012, 134, 5180–5186. [Google Scholar] [CrossRef]
  39. Ceballos, B.M.; Yang, J.Y. Directing the reactivity of metal hydrides for selective CO2 reduction. Proc. Natl. Acad. Sci. USA 2018, 115, 12686–12691. [Google Scholar] [CrossRef]
  40. Zhang, Q.H.; Han, W.D.; Hong, Y.J.; Yu, J.G. Photocatalytic reduction of CO2 with H2O on Pt-loaded TiO2 catalyst. Catal. Today 2009, 148, 335–340. [Google Scholar] [CrossRef]
  41. Wang, W.; An, W.; Ramalingam, B.; Mukherjee, S.; Niedzwiedzki, D.M.; Gangopadhyay, S.; Biswas, P. Size and Structure Matter: Enhanced CO2 Photoreduction Efficiency by Size-Resolved Ultrafine Pt Nanoparticles on TiO2 Single Crystals. J. Am. Chem. Soc. 2012, 134, 11276–11281. [Google Scholar] [CrossRef]
  42. Katsumata, K.I.; Sakai, K.; Ikeda, K.; Carja, G.; Matsushita, N.; Okada, K. Preparation and photocatalytic reduction of CO2 on noble metal (Pt, Pd, Au) loaded Zn-Cr layered double hydroxides. Mater. Lett. 2013, 107, 138–140. [Google Scholar] [CrossRef]
  43. Xie, S.; Wang, Y.; Zhang, Q.; Fan, W.; Denga, W.; Wang, Y. Photocatalytic reduction of CO2 with H2O: Significant enhancement of the activity of Pt–TiO2 in CH4 formation by addition of MgO. Chem. Commun. 2013, 49, 2451–2453. [Google Scholar] [CrossRef] [PubMed]
  44. Xiong, Z.; Lei, Z.; Kuang, C.C.; Chen, X.; Gong, B.; Zhao, Y.; Zhang, J.; Zheng, C.; Wu, J.C.S. Selective photocatalytic reduction of CO2 into CH4 over Pt-Cu2O TiO2 nanocrystals: The interaction between Pt and Cu2O cocatalysts. Appl. Catal. B Environ. 2017, 202, 695–703. [Google Scholar] [CrossRef]
  45. Kočí, K.; Dang Van, H.; Edelmannová, M.; Reli, M.; Wu, J.C.S. Photocatalytic reduction of CO2 using Pt/C3N4 photocatalyts. Appl. Surf. Sci. 2020, 503, 144426. [Google Scholar] [CrossRef]
  46. Tasbihi, M.; Kočí, K.; Edelmannová, M.; Troppová, I.; Reli, M.; Schomäcker, R. Pt/TiO2 photocatalysts deposited on commercial support for photocatalytic reduction of CO2. J. Photochem. Photobiol. A Chem. 2018, 366, 72–80. [Google Scholar] [CrossRef]
  47. Xu, J.; Liu, X.; Zhou, Z.; Deng, L.; Liu, L.; Xu, M. Platinum Nanoparticles with Low Content and High Dispersion over Exfoliated Layered Double Hydroxide for Photocatalytic CO2 Reduction. Energy Fuels 2021, 35, 10820–10831. [Google Scholar] [CrossRef]
  48. Li, X.; Bi, W.; Zhang, Q.; Tao, S.; Chu, W.; Zhang, Q.; Luo, Y.; Wu, C.; Xie, Y. Single-atom Pt as co-catalyst for enhanced photocatalytic H2 evolution. Adv. Mater. 2016, 28, 2427–2431. [Google Scholar] [CrossRef]
  49. Cao, Y.; Wang, D.; Lin, Y.; Liu, W.; Cao, L.; Liu, X.; Zhang, W.; Mou, X.; Fang, S.; Shen, X.; et al. Single Pt atom with highly vacant d-orbital for accelerating photocatalytic H2 evolution. ACS Appl. Energy Mater. 2018, 1, 6082–6088. [Google Scholar] [CrossRef]
  50. Williams, J.A.G.; Beeby, A.; Davies, S.; Weinstein, J.A.; Wilson, C. An Alternative Route to Highly Luminescent Platinum(II) Complexes: Cyclometalation with N^C^N-Coordinating Dipyridylbenzene Ligands. Inorg. Chem. 2003, 42, 8609–8611. [Google Scholar] [CrossRef]
  51. Cárdenas, D.J.; Echavarren, A.M.; Ramírez de Arellano, M.C. Divergent Behavior of Palladium(II) and Platinum(II) in the Metalation of 1,3-Di(2-pyridyl)benzene. Organometallics 1999, 18, 3337–3341. [Google Scholar] [CrossRef]
  52. Demissie, T.B.; Ruud, K.; Hansen, J.H. DFT as a Powerful Predictive Tool in Photoredox Catalysis: Redox Potentials and Mechanistic Analysis. Organometallics 2015, 34, 4218–4228. [Google Scholar] [CrossRef]
  53. Grice, K.A. Carbon dioxide reduction with homogenous early transition metal complexes: Opportunities and challenges for developing CO2 catalysis. Coord. Chem. Rev. 2017, 336, 78–95. [Google Scholar] [CrossRef]
  54. Kinzel, N.W.; Werl, C.; Leitner, W. Transition Metal Complexes as Catalysts for the Electroconversion of CO2: An Organometallic Perspective. Angew. Chem. Int. Ed. 2021, 60, 11628–11686. [Google Scholar] [CrossRef] [PubMed]
  55. Bader, R.F.W. Atoms in Molecules—A Quantum Theory; Oxford University Press: Oxford, UK, 1990. [Google Scholar]
  56. Bader, R.F.W. A Bond Path:  A Universal Indicator of Bonded Interactions. J. Phys. Chem. A 1998, 102, 7314–7323. [Google Scholar] [CrossRef]
  57. Macchi, P.; Sironi, A. Chemical bonding in transition metal carbonyl clusters: Complementary analysis of theoretical and experimental electron densities. Coord. Chem. Rev. 2003, 383, 238–239. [Google Scholar] [CrossRef]
  58. Espinosa, E.; Alkorta, I.; Elguero, J.; Molins, E. From weak to strong interactions: A comprehensive analysis of the topological and energetic properties of the electron density distribution involving X–H⋯F–Y systems. J. Chem. Phys. 2002, 117, 5529–5542. [Google Scholar] [CrossRef]
  59. Johnson, E.R.; Shahar Keinan, S.; Mori-Sánchez, P.; Contreras-García, J.; Cohen, A.J.; Yang, W. Revealing Noncovalent Interactions. J. Am. Chem. Soc. 2010, 132, 6498–6506. [Google Scholar] [CrossRef]
  60. Yin, X.; Moss, J.R. Recent developments in the activation of carbon dioxide by metal complexes. Coord. Chem. Rev. 1999, 181, 27–59. [Google Scholar] [CrossRef]
  61. Reed, A.E.; Curtiss, L.A.; Weinhold, F. Intermolecular Interactions from a Natural Bond Orbital, Donor-Acceptor Viewpoint. Chem. Rev. 1988, 88, 899–926. [Google Scholar] [CrossRef]
  62. Dapprich, S.; Frenking, G. Investigation of Donor-Acceptor Interactions: A Charge Decomposition Analysis Using Fragment Molecular Orbitals. J. Phys. Chem. 1995, 99, 9352–9362. [Google Scholar] [CrossRef]
  63. Sheng, H.; Frei, H. Direct Observation by Rapid-Scan FT-IR Spectroscopy of Two-Electron-Reduced Intermediate of Tetraaza Catalyst [CoIIN4H(MeCN)]2+ Converting CO2 to CO. J. Am. Chem. Soc. 2016, 138, 9959–9967. [Google Scholar] [CrossRef] [PubMed]
  64. Call, A.; Cibian, M.; Yamamoto, K.; Nakazono, T.; Yamauchi, K.; Sakai, K. Highly Efficient and Selective Photocatalytic CO2 Reduction to CO in Water by a Cobalt Porphyrin Molecular Catalyst. ACS Catal. 2019, 9, 4867–4874. [Google Scholar] [CrossRef]
  65. Vetere, V.; Adamo, C.; Maldivi, P. Performance of the ‘parameter free’ PBE0 functional for the modeling of molecular properties of heavy metals. Chem. Phys. Lett. 2000, 325, 99–105. [Google Scholar] [CrossRef]
  66. Adamo, C.; Barone, V. Inexpensive and accurate predictions of optical excitations in transition-metal complexes: The TDDFT/PBE0 route. Theor. Chem. Acc. 2000, 105, 169–172. [Google Scholar] [CrossRef]
  67. Adamo, C.; Barone, V. Toward reliable density functional methods without adjustable parameters: The PBE0 model. J. Chem. Phys. 1999, 110, 6158–6170. [Google Scholar] [CrossRef]
  68. Ernzerhof, M.; Scuseria, G.E. Assessment of the Perdew–Burke–Ernzerhof exchange-correlation functional. J. Chem. Phys. 1999, 110, 5029–5036. [Google Scholar] [CrossRef]
  69. Adamo, C.; Scuseria, G.E.; Barone, V. Accurate excitation energies from time-dependent density functional theory: Assessing the PBE0 model. J. Chem. Phys. 1999, 111, 2889–2899. [Google Scholar] [CrossRef]
  70. Adamo, C.; Barone, V. Toward reliable adiabatic connection models free from adjustable parameters. Chem. Phys. Lett. 1997, 274, 242–250. [Google Scholar] [CrossRef]
  71. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 16W, Revision C.01; Gaussian, Inc.: Wallingford, CT, USA, 2016.
  72. Grimme, S.; Ehrlich, S.; Goerigk, L. Effect of the damping function in dispersion corrected density functional theory. J. Comput. Chem. 2011, 32, 1456–1465. [Google Scholar] [CrossRef]
  73. Wu, P.; Chaudret, R.; Hu, X.; Yang, W. Noncovalent Interaction Analysis in Fluctuating Environments. J. Chem. Theory Comput. 2013, 9, 2226–2234. [Google Scholar] [CrossRef]
  74. Lu, T.; Chen, F. Multiwfn: A multifunctional wavefunction analyser. J. Comput. Chem. 2012, 33, 580–592. [Google Scholar] [CrossRef] [PubMed]
  75. Tomasi, J.; Mennucci, B.; Cammi, R. Quantum Mechanical Continuum Solvation Models. Chem. Rev. 2005, 105, 2999–3093. [Google Scholar] [CrossRef] [PubMed]
  76. Van Gisbergen, S.J.A.; Kootstra, F.; Schipper, P.R.T.; Gritsenko, O.V.; Snijders, J.G.; Baerends, E.J. Density-functional-theory response-property calculations with accurate exchange-correlation potentials. Phys. Rev. A At. Mol. Opt. Phys. 1998, 57, 2556–2571. [Google Scholar] [CrossRef]
  77. Jamorski, C.; Casida, M.E.; Salahub, D.R. Dynamic polarizabilities and excitation spectra from a molecular implementation of time-dependent density-functional response theory: N2 as a case study. J. Chem. Phys. 1996, 104, 5134–5147. [Google Scholar] [CrossRef]
  78. Bauernschmitt, R.; Ahlrichs, R. Treatment of electronic excitations within the adiabatic approximation of time dependent density functional theory. Chem. Phys. Lett. 1996, 256, 454–464. [Google Scholar] [CrossRef]
  79. Laurent, A.D.; Jacquemine, D. TD-DFT Benchmarks: A review. Int. J. Quant. Chem. 2013, 113, 2019–2039. [Google Scholar] [CrossRef]
  80. Latouche, C.; Skouteris, D.; Palazzetti, F.; Barone, V. TD-DFT Benchmark on Inorganic Pt(II) and Ir(III) Complexes. J. Chem. Theory Comput. 2015, 11, 3281–3289. [Google Scholar] [CrossRef]
  81. Martin, R.L. Natural Transition Orbitals. J. Chem. Phys. 2003, 118, 4775–4777. [Google Scholar] [CrossRef]
Scheme 1. The molecular structures of Pt(II) complexes under study.
Scheme 1. The molecular structures of Pt(II) complexes under study.
Molecules 29 00403 sch001
Figure 1. Proposed calculated mechanism for the CO2-into-CO photocatalytic reduction by the [Pt(5-R-dpb)Cl], 15 Pt(II) complexes.
Figure 1. Proposed calculated mechanism for the CO2-into-CO photocatalytic reduction by the [Pt(5-R-dpb)Cl], 15 Pt(II) complexes.
Molecules 29 00403 g001
Figure 2. Optimized geometries of 1 in (a) S0 ground state, (b) in T1 excited state and (c) of One-Electron-Reduced form, 1_OER, calculated at the PBE0-GD3BJ/Def2-TZVP level, in DCM solvent.
Figure 2. Optimized geometries of 1 in (a) S0 ground state, (b) in T1 excited state and (c) of One-Electron-Reduced form, 1_OER, calculated at the PBE0-GD3BJ/Def2-TZVP level, in DCM solvent.
Molecules 29 00403 g002
Figure 3. Simulated absorption spectra of 14 in DCM at the TDDFT/PBE0/Def2-TZVP level.
Figure 3. Simulated absorption spectra of 14 in DCM at the TDDFT/PBE0/Def2-TZVP level.
Molecules 29 00403 g003
Figure 4. NTO pairs for the most significant electronic transitions of the simulated absorption spectra of 1 (hole h+ on the left, electron e on the right).
Figure 4. NTO pairs for the most significant electronic transitions of the simulated absorption spectra of 1 (hole h+ on the left, electron e on the right).
Molecules 29 00403 g004
Scheme 2. Born–Haber cycle for calculation of ground-state reduction potentials E0red (above), Latimer diagram for the calculation of T1 excited state reduction potentials E0*red (middle) and Born–Haber cycle for calculation of sacrificial electron donor D (=TEOA or TEA) reduction potentials (below).
Scheme 2. Born–Haber cycle for calculation of ground-state reduction potentials E0red (above), Latimer diagram for the calculation of T1 excited state reduction potentials E0*red (middle) and Born–Haber cycle for calculation of sacrificial electron donor D (=TEOA or TEA) reduction potentials (below).
Molecules 29 00403 sch002
Figure 5. (a) Optimized geometry of 1_ImC, (b) 3D surface plot of HOMO, (c) BCPs (orange spheres) and (d) 3D surface of RDG function.
Figure 5. (a) Optimized geometry of 1_ImC, (b) 3D surface plot of HOMO, (c) BCPs (orange spheres) and (d) 3D surface of RDG function.
Molecules 29 00403 g005
Figure 6. 3D isosurface plots of (a) σ(Pt–C) BD NBO and (b) of donor and acceptor NBOs participating in the Pt–C hyperconjugative interactions.
Figure 6. 3D isosurface plots of (a) σ(Pt–C) BD NBO and (b) of donor and acceptor NBOs participating in the Pt–C hyperconjugative interactions.
Molecules 29 00403 g006
Scheme 3. Orbital interaction diagram for the formation of 1_ImC (Complex) from fragments [Pt(dpb)] (Fragm. 1) and CO2 (Fragm. 2).
Scheme 3. Orbital interaction diagram for the formation of 1_ImC (Complex) from fragments [Pt(dpb)] (Fragm. 1) and CO2 (Fragm. 2).
Molecules 29 00403 sch003
Scheme 4. (a) SOMOβ of 1_T1 species (left) and SOMO of 1_OER (right) and (b) LUMO of 1_ImF (left) and SOMOβ of 1_G (right).
Scheme 4. (a) SOMOβ of 1_T1 species (left) and SOMO of 1_OER (right) and (b) LUMO of 1_ImF (left) and SOMOβ of 1_G (right).
Molecules 29 00403 sch004
Figure 7. Free energy, ΔG (in kcal/mol), reaction profiles of 14 in DCM solvent calculated at the PBE0-GD3BJ/Def2-TZVP level (numbers in blue for 1, in green for 2, in orange for 3 and in red for 4).
Figure 7. Free energy, ΔG (in kcal/mol), reaction profiles of 14 in DCM solvent calculated at the PBE0-GD3BJ/Def2-TZVP level (numbers in blue for 1, in green for 2, in orange for 3 and in red for 4).
Molecules 29 00403 g007
Table 1. Reduction potentials (in Volts) of 14 in S0 state, E0red, in T1 state, E0*red, along with E0red of TEOA and TEA and the respective values vs. SHE calculated at the PBE0/Def2-TZVP level.
Table 1. Reduction potentials (in Volts) of 14 in S0 state, E0red, in T1 state, E0*red, along with E0red of TEOA and TEA and the respective values vs. SHE calculated at the PBE0/Def2-TZVP level.
ComplexE0redE0*redE0red vs. SHEE0*red vs. SHE 1E0*redox 2
vs. TEOAvs. TEA
1−2.481.800.324.607.437.19
2−2.461.82−0.363.927.457.21
3−2.541.74−0.343.947.377.13
4−3.480.80−1.882.406.436.19
TEOA−5.63 −1.35
TEA−5.39 −1.11
1 SHE = Standard Hydrogen Electrode. 2 Taken as E0red(complex) + E0ox(D) where E0ox(D) = −E0red(D).
Table 2. Topological and energetic properties of ρ(r) calculated at the (3, −1) bond-critical point (BCP) located at the Pt-CO2 bond of intermediates 1_ImC4_ImC.
Table 2. Topological and energetic properties of ρ(r) calculated at the (3, −1) bond-critical point (BCP) located at the Pt-CO2 bond of intermediates 1_ImC4_ImC.
SpeciesρBCP 1ρBCP 2GBCP 3VBCP 3|VBCP|/GBCPHBCP 3GBCP/ρBCP 4
1_ImC0.1240.1370.089−0.1451.629−0.0560.718
2_ImC0.1250.1380.090−0.1471.633−0.0570.720
3_ImC0.1250.1350.090−0.1471.633−0.0570.720
4_ImC0.1250.1300.088−0.1451.648−0.0570.704
1 in eÅ−3. 2 in eÅ−5. 3 in kJ mol−1 (atomic unit volume)−1. 4 in kJ mol−1 electron−1. The bond degree parameter HBCP/ρBCP represents either the covalence (HBCP < 0) or the softening (HBCP > 0) degree of the interaction.
Table 3. Natural charges Q, bonding NBOs BD and hyperconjugative interaction stabilization energy ΔE(2) relevant to the Pt–CO2 in 1_ImC4_ImC.
Table 3. Natural charges Q, bonding NBOs BD and hyperconjugative interaction stabilization energy ΔE(2) relevant to the Pt–CO2 in 1_ImC4_ImC.
SpeciesQPtQC 1BD[σ(Pt-CO2)]ΔE(2)
σ(Pt-CN^C^N) → σ*(Pt-CO2)σ(Pt-CO2) → σ*(Pt-CN^C^N)
1_ImC0.2350.4950.558hPt + 0.830hC58.327.8
2_ImC0.2370.4940.540hPt + 0.842hC58.126.4
3_ImC0.2440.4950.538hPt + 0.843hC58.327.2
4_ImC0.2600.4970.535hPt + 0.845hC56.923.6
1 Referring to the CO2 carbon atom.
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

Sarantou, A.; Tsipis, A. Photocatalytic Reduction of CO2 into CO with Cyclometalated Pt(II) Complexes of N^C^N Pincer Dipyridylbenzene Ligands: A DFT Study. Molecules 2024, 29, 403. https://doi.org/10.3390/molecules29020403

AMA Style

Sarantou A, Tsipis A. Photocatalytic Reduction of CO2 into CO with Cyclometalated Pt(II) Complexes of N^C^N Pincer Dipyridylbenzene Ligands: A DFT Study. Molecules. 2024; 29(2):403. https://doi.org/10.3390/molecules29020403

Chicago/Turabian Style

Sarantou, Antonia, and Athanassios Tsipis. 2024. "Photocatalytic Reduction of CO2 into CO with Cyclometalated Pt(II) Complexes of N^C^N Pincer Dipyridylbenzene Ligands: A DFT Study" Molecules 29, no. 2: 403. https://doi.org/10.3390/molecules29020403

Article Metrics

Back to TopTop