Tuning the degree of CO2 activation by carbon doping Cun− (n = 3–10) clusters: an IR spectroscopic study

Copper clusters on carbide surfaces have shown a high catalytic activity towards methanol formation. To understand the interaction between CO2 and the catalytically active sites during this process and the role that carbon atoms could play in this, they are modeled by copper clusters, with carbon atoms incorporated. The formed clusters CunCm− (n = 3–10, m = 1–2) are reacted with CO2 and investigated by IR multiple-photon dissociation (IR-MPD) spectroscopy to probe the degree of CO2 activation. IR spectra for the reaction products [CunC·CO2]−, (n = 6–10), and [CunC2·CO2]−, (n = 3–8) are compared to reference spectra recorded for products formed when reacting the same cluster sizes with CO, and with density functional theory (DFT) calculated spectra. The results reveal a size- and carbon load-dependent activation and dissociation of CO2. The complexes [CunC·CO2]− with n = 6 and 10 show predominantly molecular activation of CO2, while those with n = 7–9 show only dissociative adsorption. The addition of the second carbon to the cluster leads to the exclusive molecular activation of the CO2 on all measured cluster sizes, except for Cu5C2− where CO2 dissociates. Combining these findings with DFT calculations leads us to speculate that at lower carbon-to-metal ratios (CMRs), the C can act as an oxygen anchor facilitating the OC 
<svg xmlns="http://www.w3.org/2000/svg" version="1.0" width="13.200000pt" height="16.000000pt" viewBox="0 0 13.200000 16.000000" preserveAspectRatio="xMidYMid meet"><metadata>
Created by potrace 1.16, written by Peter Selinger 2001-2019
</metadata><g transform="translate(1.000000,15.000000) scale(0.017500,-0.017500)" fill="currentColor" stroke="none"><path d="M0 440 l0 -40 320 0 320 0 0 40 0 40 -320 0 -320 0 0 -40z M0 280 l0 -40 320 0 320 0 0 40 0 40 -320 0 -320 0 0 -40z"/></g></svg>
 O bond rupture, whereas at higher CMRs the carbon atoms increasingly attract negative charge, reducing the Cu cluster’s ability to donate electron density to CO2, and consequently its ability to activate CO2.


Introduction
The continuous growth of carbon emissions has led, and continues to lead, to an accumulation of CO 2 in the atmosphere, and the concomitant global warming. To mitigate this, an urgent solution is needed to reduce atmospheric CO 2 concentrations. At the same time, CO 2 is the largest carbon source on Earth, and it could be used as a feedstock for many valuable chemicals, such as alcohols and other higher hydrocarbons. The bottleneck lies in the high kinetic and thermal stability of the CO 2 molecule, expressed in the C]O bond energy of more than 7 eV. 1 One promising way for CO 2 utilization is its hydrogenation to simple alcohols like methanol. Industrially this process takes place over a Cu/ZnO/Al 2 O 3 catalyst at temperatures of 200-300 C and pressures of 50-100 bar. 2 The energy required to drive this reaction leads to additional CO 2 emissions and elevated costs, while the one-pass selectivity towards methanol formation is rather low, with larger concentrations of CO and H 2 O formed. Thus, for direct CO 2 conversion, a more selective catalyst is required. Many different catalyst materials have been tested to reduce reaction temperatures and increase the methanol yield.
One issue that seems not to affect relatively well-functioning CO 2 hydrogenation catalysts is poisoning due to coke formation, a process that hampers the reverse reaction during methanol decomposition. 3 Rather, metal carbide catalyst materials have shown enhanced activity in CO 2 activation. This enhancement has been rationalized by increased catalyst stability upon adding carbon to the metal, and by the modication of the metal's electronic and structural properties, changing its catalytic activity. 4,5 Transition metal carbides (TMCs) have even demonstrated chemical activity similar, and sometimes even better, than that of the platinum group metals in the transformation of hydrocarbons and oxygencontaining molecules. [6][7][8] TMCs are also considered to be promising catalysts for CO 2 conversion by H 2 . 8 However, experimentally, it was shown that, for example, Mo 2 C is quite aggressive and leads to CO 2 dissociation, forming CO and H 2 O, via the reverse water gas shi (rWGS) reaction. 9,10 Based on density functional theory (DFT) calculations, it is suggested that Mo 8 C 12 nanoparticles are better catalysts, considering their high stability and moderate chemical activity in comparison to that of bulk Mo (too low) and Mo 2 C (too high). 11 Similar conclusions were drawn for Ti carbides. 11 Some transition metal catalysts with a carbon-to-metal ratio (CMR) below 1 easily break the C-O bond leading to CO formation. 6,7 An increase in carbon content may enhance catalyst stability, it also leads to a decrease of the chemical activity, for instance observed for TiC and MoC, due to the ligand effect expressed in a downshi of the metal d-band center affecting the charge transfer from the metal to CO 2 and, thus, also adsorption energy. 7,[11][12][13] Accordingly, the activity of a metal carbide catalyst towards CO 2 dissociation can be tuned by varying the CMR of the catalyst. 14 TMC catalytic activity can also be altered by the deposition of other metals on the TMC surface. 5,7,10,14 Cu-promoted carbides exhibit enhanced selectivity towards methanol formation in comparison to both bare copper and bare TMC surfaces, like TaC, SiC, TiC, Mo 2 C, and MoC. 5,7,14 It is also pointed out that the Cu-TMCs interface plays a crucial role. The hydrogenation of the reaction intermediates, such as HCO and H 2 CO, is energetically more favorable on the Cu sites, making small Cu clusters more suitable catalysts for methanol formation. 5,14 Cu 4 / TiC has an even higher CO 2 adsorption energy, and therefore methanol formation rate, than the commercial Cu/ZnO catalyst. 5,15,16 However, most Cu/TMC catalyzed reactions proceed via CO formation and its subsequent hydrogenation, rather than via direct CO 2 hydrogenation and a formate intermediate, as is the widely accepted reaction pathway over Cu catalysts. 2 The high selectivity towards CO results in an elevated CO yield at the expense of methanol production. Therefore, it is important to understand what the direct inuence of carbon is on the deposited clusters and how it inuences CO 2 activation.
The ne details of a reaction proceeding on the active sites can be studied at the molecular level using isolated clusters. Mass-selective gas-phase spectroscopy allows studying the interaction of CO 2 with clusters of well-dened elemental composition. For example, photoelectron and IR spectroscopy have been used to study the interaction between CO 2 and metal anions. [17][18][19][20][21][22][23][24][25][26][27] This interaction can be characterized by four binding motifs: bidentate (h 2 ), metal formate, oxalate, and dissociative adsorption. The variety in CO 2 adsorption motifs is related to the difference in the degree of electron transfer from the metal center into the antibonding 2p u LUMO of CO 2 . 28 For Cu À and other coinage metal ions in particular, it was found that CO 2 binds in a formate-like fashion, where the metal ion replaces the hydrogen. 17,19 The reaction of CO 2 with isolated anionic metal clusters was also studied with ow tube reactor mass spectrometry, 29,30 and IR spectroscopy. 31,32 The nature of CO 2 adsorption on these clusters appears to be sizedependent, for example on Pt n À clusters, where Pt 4 À shows CO 2 activation, whereas Pt 5 À , Pt 6 À and Pt 7 À show dissociation. 31 Another example is Co n À , where the Co À ion binds two CO 2 molecules in a bidentate conguration, 27 and clusters with n > 7 dissociate CO 2 . 32 In this work, we investigate how the carbon doping of anionic Cu clusters will inuence the degree of CO 2 activation, where the anionic charge state was chosen to facilitate electron transfer to the CO 2 molecule. Anionic copper clusters containing 5-10 Cu atoms were doped with one or two C atoms and reacted with CO 2 . The formed [Cu n C m $CO 2 ] À species, with n ¼ 3-10 and m ¼ 1-2, were studied by mass-selective IR spectroscopy. The obtained spectra are compared to reference IR spectra of the clusters reacted with CO and with DFT calculated spectra.

Methods
The experiment is shown schematically in Fig. 1. Carbon-doped anionic copper clusters Cu n C m À are produced in a Smalley-type laser ablation source. 33 An isotopically enriched Cu-65 (99.9%) target is ablated by the focused second harmonic of a Nd:YAG laser (532 nm, 14 mJ) in the presence of a He carrier gas that is injected by a pulsed valve (General Valve series 9, stagnation pressure 6 bar). Carbon-doped copper clusters were formed due to the presence of trace carbon in the carrier gas or residual pump oil vapour. The formed clusters are reacted with CO 2 or CO by injecting the pure gases through a second pulsed valve with a 1 bar stagnation pressure into an extension of the clustering channel. The Cu n C m À clusters react with CO 2 or CO at room temperature in a reaction channel, forming [Cu n C m $CO 2 ] À or [Cu n C m $CO] À complexes with n ¼ 3-10, and m ¼ 1-2. The bracket notation implies that no prior knowledge of the adsorption form of CO or CO 2 is assumed. The gas mixture of helium and clusters is expanded into vacuum, forming a molecular beam that is collimated by a 2 mm diameter skimmer and an 8 Â 2 mm slit aperture. The ions are then irradiated by IR laser light, and subsequently pulse-extracted by high voltage plates into an orthogonal time-of-ight mass spectrometer (TOF-MS), where they are mass separated and detected by a multichannel plate (MCP) detector. For the spectroscopic experiments, the molecular beam overlaps between the extraction plates at a 35 angle with the IR laser beam of the free-electron laser for intracavity experiments (FELICE). 34 IR light is produced in the form of macropulses, consisting of an approximately 10 ms pulse train of transform-limited ps duration micropulses, and is characterized by a spectral bandwidth of $0.7% of the central frequency. The frequency range covered in this work is 320-2120 cm À1 , with typical macropulse energies of 0.06-0.8 J, and uences of 0.7-6.1 J cm À2 . The experiment is operated at twice the IR laser frequency, and mass spectra with and without irradiation are recorded to correct for the uctuations during cluster production. IRMPD spectra are obtained by calculating the depletion yield Y(n) at frequency n, dened as: where I on (n) is the ion intensity of the irradiated complex, I off the ion intensity without irradiation, and P the macropulse energy. Spectra displaying the depletion I on (n)/I off can be found in Fig. S1 † of the ESI. For the structural assignment of selected spectra (for Cu 6 C À and Cu 7 C À ), DFT calculations of different isomers were done using the Q-Chem 5.3 program package. 35 Stable structures were identied using the TPSSh/def2-TZVP + D3 level of theory, 36 as described in the ESI. † A detailed description of the search procedure is given in the ESI. † This combination of functional, basis set and dispersion correction was chosen to be able to directly compare the results with our previous work on the adsorption of CO 2 on cationic Cu clusters, 37 where adsorption was

Paper
Faraday Discussions limited to physisorption. The accuracy of this method was carefully evaluated compared to a CCSD(T)/def2-QZVPPD benchmark. 37 In other studies it has been found that the precise mode of physisorption found in these calculations can markedly be inuenced by the choice of dispersion correction. 38 However, since the interaction for the species under study here is signicantly stronger, as seen by the activated and dissociated products discussed below, we believe that these inuences are not decisive here. This is conrmed by further DFT based computations with D3 or D4 and without dispersion correction (Table S1 in the ESI †).
Only the lowest spin multiplicities were considered in the calculations, so doublet for Cu 6 C À , and singlet for Cu 7 C À . Harmonic vibrational frequencies of these structures are convoluted with 20 cm À1 Gaussian line shape function and compared to the experimental data for nal assignment. No frequency scaling has been applied for the comparison.

Results and discussion
To understand how CO 2 binds to anionic carbon-doped copper clusters, IRMPD spectra of products formed upon reacting Cu n C m À (n ¼ 3-10) with different carbon loading m ¼ 1-2 are recorded. A mass spectrum of all species formed is shown in Fig. 2. From this gure, one can see that for each size n, a distribution of Cu n C m À clusters is formed, that for larger clusters (n $ 9) is essentially limited to m ¼ 0-3. Pure Cu n À clusters are lower in intensity than the carbon doped clusters, which we attribute to a more facile formation of a CuC À seed as a nucleation core than Cu 2 À . In the inset, a close-up of the mass spectrum in the region close to Cu 6 À and Cu 7 À is shown, with individual mass peaks corresponding to Cu n Cm (CO 2 ) p À labeled by (n and m) or (n, m and p). It can be seen that the pure Cu Fig. 2 Mass spectrum of the products formed when reacting anionic Cu n C m À clusters with CO 2 . Mass peaks in the inset are labeled by (n and m) or (n, m, and p) for Cu n C m (CO 2 ) p À .

Faraday Discussions Paper
clusters are produced in this mass spectral region too, albeit at lower intensity than the Cu n C À , Cu n C 2 À , and Cu n C 3 À signals. Unfortunately, CO 2 adsorption on the pure Cu n À clusters is negligible, precluding spectroscopic characterization.

CO 2 activation by Cu n C À
The IRMPD spectra of low-carbon loading products [Cu n C$CO 2 ] À (n ¼ 6-10) are shown in Fig. 3 (panels a-e). The complexes with n ¼ 7-9 show only one main band around 2000 cm À1 , which is fairly close to the frequency of the free C-O stretch vibration at 2143 cm À1 . 39 We can thus speculate that these cluster sizes adsorb CO 2 dissociatively, leading to the formation of CO. In contrast, the spectra of complexes with n ¼ 6 and 10 are dominated by two bands around 720 and 1630 cm À1 , with a weaker band visible around 1100 cm À1 . Earlier, experimental IR studies of the CO 2 À radical embedded in alkali halide matrices have revealed a characteristic band at 1671 cm À1 . 40 The presence of a band in this region in the IRMPD spectra suggests that Cu n C À clusters with n ¼ 6 and 10 induce a distortion of CO 2 , activating it. Bands at similar frequencies were observed when the Pt 4 À ion was reacted with CO 2 , which was also concluded to lead to activation. 31 From here, we focus on the species [Cu n C$CO 2 ] À with n ¼ 6 and 7 since their spectra are representative of two modes of adsorption, activated and dissociated. The suspicion that CO 2 adsorbs dissociatively on CCu 7 À can be veried by reference spectra of carbon monoxide (CO) adsorbed to Cu 7 C À , which is compared to the IRMPD spectrum of [Cu 7 C$CO 2 ] À in Fig. 4b (second panel from the top). For this, we reacted the clusters with CO and recorded IRMPD spectra in the 650-2100 cm À1 range for the [Cu n C$CO] À (n ¼ 4-10) species formed. The resulting spectra for all measured cluster sizes are dominated by a strong band around 2020 cm À1 , characteristic for the C-O stretch. This band shows an excellent agreement with the bands observed when Cu 7 C À is reacted with CO 2 (Fig. 4b, top panel), strengthening the hypothesis that CO 2 adsorbs dissociatively on Cu n C À (n ¼ 7-9). In contrast, the mismatch between spectra for [Cu 6 C$CO 2 ] À Fig. 4 IRMPD spectra of: (a) Cu 6 C À , and (b) Cu 7 C À , reacted with CO 2 (red) and CO (blue), and calculated spectra (black) of the lowest energy structure. All calculated spectra are convoluted with a 20 cm À1 FWHM Gaussian line shape function, and are accompanied by the geometrical structure (Cu, C and O atoms represented by orange, black and red spheres) and relative energy. The green curves are vertical zooms with the multiplication factors indicated. The CO 2 binding energy is given for the lowest energy structure, relative energies for other isomers. and [Cu 6 C$CO] À (Fig. 4a) suggests that CO 2 adsorbs predominantly molecularly on Cu 6 C À and, by extension, also on Cu 10 C À . DFT calculations were performed to further interpret these observations. For this, we rst searched for bare Cu n C m À clusters. In our global optimization routine, we have identied the lowest energy structures shown in Fig. 5. The lowest energy structures found are not unlike those predicted for the pure clusters. Cu 6 C À is similar to the predicted Cu 6 À octahedron, 41 with the C substituting one of the Cu atoms in the octahedron, and the sixth Cu atom bound to one of the Cu-Cu edges. Cu 7 C À looks like the boat structure proposed for Cu 6 À , 42 with a C bound on a hollow site, and the last Cu atom bound both to the Cu-Cu edge and to the C atom. Interestingly, the structures for Cu 5 C 2 À and Cu 6 C 2 À are planar and are similar to the structure predicted for pure Cu 5 À , a 2D trapezoid. 41,42 In both clusters, the C atoms are bound together in a C 2 unit that forms one corner of a planar hexagon for Cu 6 C 2 À , with one Cu atom missing for Cu 5 C 2 À .
Using the structures found for the bare clusters, a search was done for structures of the reaction product of Cu 6 C À and Cu 7 C À . An extensive search has identied over sixty stable structures for each species. Many of these structures have a similar binding motif of C and O atoms to the cluster leading to only minor differences in IR spectra, which are indistinguishable in the currently applied experimental method. Therefore, the structures were grouped, rst using their spectral properties only. Visual inspection of all structures per group allowed to conclude that they all shared similar binding motifs. For the CO 2 reaction products with Cu 6 C À , sixteen unique spectral patterns were identied, for Cu 7 C À only nine. The spectra of the lowest energy representative from each group are compared to the experimental data in Fig. S2. † A selection of the most promising candidates is shown in Fig. 4 (black traces).
In the search, structures are found both with dissociated and molecular CO 2 , where structures with CO 2 dissociated are generally lowest in energy. The spectra of structures with dissociated CO 2 are dominated by high-frequency bands originating from the C-O stretch vibration(s), with frequencies ranging from 1800- Fig. 5 Lowest energy structures of the single C-doped clusters Cu 6 C À and Cu 7 C À , and double C-doped clusters Cu 5 C 2 À and Cu 6 C 2 À .

Paper
Faraday Discussions 2100 cm À1 , depending on the binding site of the CO. Spectra of clusters with molecularly bound CO 2 have multiple bands typically at the lower frequencies associated with the CO 2 bending mode, and the C-O bound and C-O free stretch. For both n ¼ 6 and 7, the largest number of structures found is with CO 2 dissociatively adsorbed, where CO is bound in an on-top conguration (m 1 ) to a Cu atom, and the O eliminated from CO 2 bound to another Cu; an example is structure 7E. However, this binding motif is not the energetically most favorable lying at least 1.78 (2.19) eV higher in energy for n ¼ 6 (7) than the lowest energy structure. We attribute this simply to a computational sampling effect, where the probability to nd a Cu atom in a Cu 6 C À cluster is six times higher than nding a C atom.
For n ¼ 6, the spectra of structures from this most populated group are very similar to that of structure 6A, the lowest energy structure, whereas this has a different structural motif. In structure 6A, the cluster structure contains a Cu 5 square pyramid, with the sixth Cu bridge-bound (m 2 ) to one of the edges, and the C bound to a hollow (m 3 ) Cu 3 site. Of the dissociated CO 2 , the CO is m 1 -bound to the cluster's C atom, forming a linear C-C-O moiety, and the O to a hollow Cu site. The spectrum has its main band, the CO stretch vibration, close to 2100 cm À1 , or blue-shied from the observed band frequency by about 60 cm À1 . The same motif of CO binding to the lone C is found for the lowest energy group for n ¼ 7 (e.g. structure 7A). 7A is a Cu 6 octahedron capped by the seventh Cu and the C, both in m 3 conguration. The CO 2 is dissociated with the single O capping a third octahedron plane, and the CO bound to the C dopant.
For both n ¼ 6 and 7 the second group (e.g. 6B and 7B) is formed by structures where CO 2 dissociation leads to an O atom migrating to the C dopant atom, forming two CO molecules that are m 1 -bound to a Cu. For 6B this is in the form of a Cu 6 boat structure with the CO molecules attached to the bow and the stern, while for 7B it is a Cu 7 pentagonal bipyramid. The frequencies of both CO stretch vibrations are close to each other resulting in one merged vibrational band, which is found at frequencies of about 2000 cm À1 , slightly lower than that for the C-C-O group in 6A and 7A. Group 3 (e.g. 6C/7C) is characterized by two CO molecules, with one m 1 -and the other m 2 -bound. In 7C, the same capped octahedron is recognized, with the bridge-bound CO close to the on-top CO. In this case, two distinct C-O stretch bands are seen, one around 2000 cm À1 for the m 1 bound CO, and one around 1900 cm À1 for the m 2 bound CO. For n ¼ 6, group 4 has CO molecules bound in m 2 and m 3 congurations, with again clearly two distinct bands, now at 1800 and 1900 cm À1 , consistent with a weakening of the C-O bond when it coordinates to multiple Cu atoms. An interesting motif found is the linking of two CO molecules resulting in an O-C-C-O chain. However, for both cluster sizes, this last group is lowly populated, relatively high in energy, and presents no match for the experimental spectrum.
Structures with molecularly bound CO 2 are typically much higher in energy (>1.4 eV) for both n ¼ 6 and 7. The structures found are always di-s-bound via the C and the O. Formation of a C-C bond is favorable resulting in the preferential binding of the CO 2 carbon atom to the cluster's carbon atom, with one of the O atoms binding to a Cu atom. The lowest energy structures displaying this binding motif are 6E and 7D.
Structures are assigned based on the comparison between calculated and experimental spectra. The IRMPD spectrum of [Cu 6 C$CO 2 ] À shows two major bands at 718 and 1630 cm À1 with weaker bands at 1094, 1426, 1802, and 2016 cm À1 . This number of bands cannot be explained by a single structure. The major band at 718 cm À1 agrees best with the CO 2 bending mode predicted at 705 cm À1 for structure 6F, where CO 2 is intact, but bent as a result of charge transfer. The CO 2 molecule is bound with a C and the O to two neighboring Cu atoms in a bidentate bridging conguration. Two more intense bands are predicted at 1141 and 1620 cm À1 offering a good match for the experimental bands at 1094 and 1630 cm À1 . These correspond to the C-O stretching vibrations involving the bound and free O atoms, respectively. Although the band frequencies of 6F match the experiment very well, the relative intensities are less convincing. The experimental band at 718 cm À1 has the highest intensity, even though it does not dominate any of the calculated spectra. We do not have a good explanation for this observation; we can speculate that it is related to the excitation mechanism in IRMPD spectroscopy, which may cause discrepancies with calculated intensities.
The search for low energy structures for Cu 6 C À yielded a Cu 5 square pyramid structure shown in Fig. 5, with the C atom fourfold coordinated to the base, and an additional Cu atom m 2 -bound on one of the edges. In structure 6F, this structure is retained, making this a likely entrance complex. Despite the higher energy of structure 6F (+1.86 eV relative to the lowest energy structure) it is the easiest formed. The formation of other more stable isomers necessarily proceeds via either C-O bond cleavage or structural rearrangement of the cluster. Both processes require additional energy to overcome a barrier associated with the transition state.
The depletion spectrum of [Cu 6 C$CO 2 ] À (Fig. S1 †) does not rule out a mixture of different isomers present in the molecular beam, since none of the bands reach 100% depletion. The enhanced intensity of the 1630 cm À1 experimental band could then potentially be explained by the presence of an isomer like 6E with predicted bands at 1163 and 1629 cm À1 . Structure 6E appears also relatively easily formed as an entrance complex, where CO 2 also binds in a bidentate bridging conguration, now with the C to the cluster's carbon atom. However, the bare cluster's bridging Cu is now bound on a hollow site.
The dominant experimental bands can thus be explained by molecularly adsorbed CO 2 , but minor bands above 1750 cm À1 more likely originate from the carbonyl C-O stretch, resulting from CO 2 dissociation. The band at 1802 cm À1 might be due to structure 6C's 1825 cm À1 predicted mode but could also be assigned to the doublet from structure 6D at 1764 and 1825 cm À1 , which could have merged into one band due to band broadening. The 2016 cm À1 band could be the second band of structure 6C at 1990 cm À1 or 6A and 6B's bands at 2069 and 2018 cm À1 , respectively. The offset for 6A makes 6B, 6C and 6D the more likely candidates to explain the high-frequency bands.
All in all, it can be concluded that CO 2 binds mostly molecularly to Cu 6 C À , but in an activated form via the C atom either to C or to Cu atoms. Dissociation may occur, but it is certainly not the dominant motif, which is underlined by the signicantly higher predicted IR intensities of the carbonyl C-O stretch modes.
This dominance of molecular binding is in sharp contrast to what happens aer adsorption of CO 2 on Cu 7 C À . The IRMPD spectrum of [Cu 7 C$CO 2 ] À has one dominant band at 2033 cm À1 and a low-intensity band at 598 cm À1 . From a comparison with low-energy isomers, it is quite clear that CO 2 is dissociated since no bands associated with activated complexes (7D and 7F) are present in the experimental spectrum. Structure 7C has two CO binding motifs, and its spectrum has separate bands at 1915 (m 2 , Cu-Cu) and 2017 cm À1 (m 1 , Cu). Because this does not match the experimental spectrum, we must look for structures with only m 1 binding. Spectra calculated for such structures (7A, 7B, and 7E) are dominated by one band, depending on the adsorption site of CO, at 2078 (m 1 to C), 2017 (twice m 1 to Cu) and 2038 cm À1 (m 1 to Cu), respectively. Each of them could explain the main experimental band at 2033 cm À1 . The low-intensity band at 598 cm À1 could be due to the weak bands around 558 and 600 cm À1 for structures 7A and 7E, respectively. These bands are more than twenty times lower in intensity than the C-O stretch. We therefore assign the spectrum to either 7A or 7E, corresponding to structures where the CO formed from dissociating CO 2 is m 1 -bound either to the cluster's C or Cu atom, and O bound separately to the cluster. Both structures appear not easily formed, because the lowest energy structure for Cu 7 C À (Fig. 5) requires substantial re-arrangement. CO 2 activation by Cu n C 2 À If CO 2 is reacted with Cu clusters containing a second carbon, the picture changes. The spectra for the [Cu n C 2 $CO 2 ] À (n ¼ 3-8) species are presented in Fig. 3 (right panels). Due to a different production efficiency for clusters with two C atoms, a slightly different range of n is presented. In contrast to the Cu n C À species, where only [Cu 6 C$CO 2 ] À and [Cu 10 C$CO 2 ] À show molecular adsorption, all spectra except that for [Cu 5 C 2 $CO 2 ] À show bands around 740, 1150, and 1650 cm À1 indicating activated, but molecular adsorption of CO 2 . For all these species, bands above 1800 cm À1 , attributed to dissociative CO 2 adsorption, have disappeared, suggesting that only molecularly bound CO 2 complexes are le. Only the spectrum of [Cu 5 C 2 $CO 2 ] À shows an intense band at 2034 cm À1 , suggesting dissociative adsorption. From this, it can be concluded that, overall, the addition of the second C atom reduces the cluster reactivity.
Again, we selected two cluster sizes for DFT calculations, [Cu 5 C 2 $CO 2 ] À and [Cu 6 C 2 $CO 2 ] À . From the glancing overview in Fig. 3 they are representative for dissociative and molecular binding of CO 2 . The comparison of the spectra for these structures with spectra of the clusters reacted with CO conrms this conclusion, because the 2034 cm À1 band for [Cu 5 C 2 $CO 2 ] À perfectly matches that observed for [Cu 5 C 2 $CO] À . In contrast, the spectrum for [Cu 6 C 2 $CO] À does not show any counterpart for the bands dominating the spectrum of [Cu 6 C 2 $CO 2 ] À .
The calculated structures were again categorized following their spectral and adsorption motifs, as illustrated in Fig. 6. The full list of the lowest energy structures of the identied groups can be found in Fig. S3. † Interestingly, for higher carbon loading, the most common motif found is still the dissociation of CO 2 to CO and O (5C and 6.2C), both bound to Cu atoms. The lowest energy structures 5A and 6.2A are formed by the adsorption of CO to one of the cluster's C atoms and O to the other, forming a C-C-O and a CO moiety. These structures are comparable to the CO motif of Cu 6 C À and Cu 7 C À , suggesting it is overall more favorable for CO 2 to dissociate and bind to the C atoms of the cluster, rather than to Cu.
The experimental spectrum of [Cu 5 C 2 $CO 2 ] À is dominated by the strong band at 2035 cm À1 and has two more bands at 1667 and 1810 cm À1 . The maximum depletion observed ($45%) allows room for assignment to multiple isomers. All structures shown could potentially be present in the molecular beam since they all have bands that might overlap with the experimentally detected. The experimental band at 2035 cm À1 could be explained either by structure 5A with two close by bands predicted at 2077 and 2027 cm À1 , or by structures 5C, 5D, and 5E with bands at 2063, 2078, and 2012 cm À1 , respectively. The second band, at 1810 cm À1 , could be explained by the main band of structure 5F at 1865 cm À1 , or by the minor bands predicted for structures 5C and 5E. The only serious contender for the weakest band observed at 1667 cm À1 is structure 5B with the band at 1698 cm À1 . All except for 5B, have CO 2 adsorbed dissociatively on the cluster. 5B has molecularly adsorbed CO 2 bound with the one of the C atoms attached to one of the cluster's C atoms, and the O to a Cu.
If we compare the structures shown in Fig. 6a to the lowest energy structure for Cu 5 C 2 À , predicted to have a 2D, wheel-like structure, it appears that all structures apart from 5A could be formed without all too large structural re-arrangements of the original cluster: all retain a planar structure with the two C atoms of the bare clusters close together. Only the formation of 5A upon CO 2 adsorption requires a major structural re-arrangement. In contrast to the spectrum of [Cu 5 C 2 $CO 2 ] À , the spectrum of [Cu 6 C 2 $CO 2 ] À does not exhibit any bands above 1780 cm À1 , which is indicative that CO 2 is adsorbed molecularly. Therefore, structures 6.2A, 6.2B and 6.2C can be ruled out since they all have high-intensity bands at 1800 cm À1 or above. Three of the four Fig. 6 IRMPD spectra of (a) Cu 5 C 2 À and (b) Cu 6 C 2 À reacted with CO 2 (red) and CO (blue), and calculated spectra (black) of the lowest energy structure. For details, see caption Fig. 4.

Paper
Faraday Discussions experimental bands (742, 1173, and 1648 cm À1 ) could quite well be explained by isomer 6.2D with bands at 748, 1179, and 1685 cm À1 . The band at 1648 cm À1 could also be explained by the 1648 cm À1 band of structure 6.2E or by structure 6.2F's 1634 cm À1 band, but these structures lack intense bands at lower frequencies. Only the experimental band at 1484 cm À1 , a shoulder of the 1648 cm À1 band, cannot be explained by isomer 6.2D. Interestingly, it appears quite similar to the weaker band at 1426 cm À1 in the [Cu 6 C$CO 2 ] À spectrum, which may have a similar origin. Another parallel drawn from the spectrum of [Cu 6 C$CO 2 ] À is that the relative intensities of the middle band are not reecting the calculated ratios.

Reactive potential energy surface
It is now of interest to understand what factors could determine dissociation or lack thereof. For this, we attempted to calculate reaction pathways for two representative clusters. The complexity associated with reconstructing the reactive pathway of CO 2 adsorption and subsequent activation over carbon-doped Cu clusters is enormous. First, as experimental evidence for the structures of the bare clusters is unknown, the starting point is ill-dened. Nevertheless, using the minima found in Fig. 5, we calculated a reaction pathway for CO 2 dissociation. It is clear that the structural phase-space for CO 2 initial adsorption is large, and the pathways to dissociation are plenty. We therefore limit ourselves to potential pathways that could illuminate whether there could exist a difference in barrier for even or odd numbers of Cu atoms.
In the top panel of Fig. 7, the adsorption and subsequent dissociation of CO 2 over Cu 6 C À is shown. Adsorption leads to the formation of structure 6F, already thought to be an entrance complex, with an adsorption energy of 0.73 eV. The CO 2 is adsorbed via the C to one of the Cu atoms forming the base of the octahedron, and the O atom to the capping extra Cu atom. From here, CO is abstracted towards the base Cu atom, aer which it is bound to the next Cu-Cu bridge, overall gaining little energy with respect to the entrance complex. The transition state associated with abstracting the CO is relatively high at 1.18 eV above the reactants. Such a barrier is insurmountable, certainly under the roomtemperature near-thermal conditions of the reaction channel used here.
The bottom panel of Fig. 7 shows the same reaction path, now calculated for CO 2 adsorbing onto Cu 7 C À . The entrance complex found is much higher in energy than anything shown in Fig. 4b, at about 2.5 eV higher than structure 7A, the lowest energy structure. The cluster is an octahedral Cu 5 C, with the sixth Cu (denoted Cu(6)) bridge-bound to a Cu-Cu edge, and the seventh to the newly formed Cu-Cu(6) edge. CO 2 is adsorbed via its C atom to an octahedron Cu atom, and via an O atom to the other Cu-Cu(6) edge. CO abstraction then leads to an energetically still not very favorable structure with the CO bound on top of a Cu atom, and the O atom on a hollow Cu-Cu-Cu site. The initial adsorption energy of CO 2 onto Cu 7 C À is with 0.93 eV not much different from that of Cu 6 C À . In contrast, the dissociation reaction over Cu 7 C À is much more facile than over Cu 6 C À , and with an energy barrier 0.55 eV lower than the energy of the reactants, this dissociation reaction is well possible.
These two reaction paths are by no means a complete description of the reactions taking place. However, they do allow to draw initial conclusions. (1) The observation of dominant molecular adsorption in the experimental spectrum of [Cu 6 C$CO 2 ] À is consistent with the high barrier calculated from structure 6F in Fig. 7; (2) nding a dissociation barrier below the energy of the reactants for Cu 7 C À demonstrates that dissociation should occur, also consistent with the observed spectrum. It may not be this pathway, as lower barrier pathways may exist, but if an exothermic reaction with a barrier below the reactants exists, the system should nd it.
So what is now the inuence of the carbon atom on the dissociation propensity? We recall that for Cu n C À the dominant adsorption form is dissociative and for Cu n C 2 À molecular. If extra carbon atoms have a decisive inuence, one would compare the outcome of the CO 2 adsorption reaction for clusters with the same number of atoms, where Cu atoms are exchanged for C atoms. For this comparison, we evaluate the spectra in Fig. 3, comparing Cu 6 C À with Cu 5 C 2 À , Cu 7 C À with Cu 6 C 2 À , and so on. If anything, one sees that this comparison gives opposite outcomes: Cu 6 C À molecular, Cu 5 C 2 À dissociative; Cu 7 C À dissociative, Cu 6 C 2 À molecular; Cu 8 C À dissociative, Cu 7 C 2 À inconclusive; Cu 9 C À dissociative, Cu 8 C 2 À molecular. Then, if we compare clusters with the same number of Cu atoms: Cu 6 C À and Cu 6 C 2 À molecular, Cu 7 C À dissociative, Cu 7 C 2 À inconclusive; Fig. 7 Reactive potential energy surface for CO 2 dissociation over the lowest energy structures found for Cu 6 C À and Cu 7 C À , calculated at the TPSSh/def2-TZVPD level of theory. Energies (in eV) are given with respect to the bare clusters and free CO 2 .
Cu 8 C À dissociative, Cu 8 C 2 À molecular. The only trend, if one may call it this: the addition of the second carbon to the cluster inhibits the dissociative adsorption for all cluster sizes except for n ¼ 5. We can only speculate that the addition of an extra C atom, which has the higher electron affinity of C (2.55) compared to Cu (1.90), 43 leads to a decrease of the charge located on the Cu atoms thereby reducing the electron donation from the cluster into the CO 2 antibonding orbital.
Of course, such a hypothesis requires a more detailed investigation, including spectroscopic characterization of CO 2 activation by the pure Cu n À clusters that was elusive in this work.

Conclusion
IRMPD spectroscopy of the CO 2 adsorbed on the anionic copper carbides reveals a size-dependent bonding nature. Single-carbon doped Cu clusters Cu n C À with n ¼ 6 and 10 show predominantly activated adsorption, but dissociative adsorption of CO 2 for n ¼ 7-9. These results were veried with control experiments of CO adsorbed on the same cluster sizes and with DFT calculations. The exact position of the C and O atoms on the clusters could not be determined precisely, since the bands for specic structures are very similar in frequency. Therefore, only the main binding motifs have been analyzed. Doubly-carbon doped Cu clusters Cu n C 2 À with n ¼ 3-8 show predominantly molecular adsorption, excepting Cu 5 C 2 À which shows dissociative adsorption of CO 2 . DFT calculations of the dissociation reaction of CO 2 over Cu 6 C À and Cu 7 C À are consistent with the observed IR spectra, showing only activation for the former and dissociation for the latter. To understand why the addition of a second C atom apparently throttles CO 2 dissociation requires further extensive computational work, evaluating reactive pathways over other cluster compositions. We speculate that the higher electron affinity of C relative to that of Cu will reduce the capacity of Cu atoms in higher C/Cu ratio clusters to donate electron density to CO 2 , reducing the activation capabilities.

Conflicts of interest
There are no conicts to declare.