Open Access Article
This Open Access Article is licensed under a
Creative Commons Attribution 3.0 Unported Licence

Cluster size dependent coordination of formate to free manganese oxide clusters

Sandra M. Lang *a, Thorsten M. Bernhardt a, Joost M. Bakker b, Robert N. Barnett c and Uzi Landman c
aInstitute of Surface Chemistry and Catalysis, University of Ulm, Albert-Einstein-Allee 47, 89069 Ulm, Germany. E-mail: sandra.lang@uni-ulm.de
bRadboud University, Institute for Molecules and Materials, FELIX Laboratory, Toernooiveld 7, 6525 ED Nijmegen, The Netherlands
cSchool of Physics, Georgia Institute of Technology, Atlanta, GA 30332-0430, USA

Received 22nd August 2023 , Accepted 14th November 2023

First published on 14th November 2023


Abstract

The interaction of free manganese oxide clusters, MnxOy+ (x = 1–9, y = 0–12), with formic acid was studied via infrared multiple-photon dissociation (IR-MPD) spectroscopy together with calculations using density functional theory (DFT). Clusters containing only one Mn atom, such as MnO2+ and MnO4+, bind formic acid as an intact molecule in both the cis- and trans-configuration. In contrast, all clusters containing two or more manganese atoms deprotonate the acid's hydroxyl group. The coordination of the resulting formate group is strongly cluster-size-dependent according to supporting DFT calculations for selected model systems. For Mn2O2+ the co-existence of two isomers with the formate bound in a bidentate bridging and chelating configurations, respectively, is found, whereas for Mn2O4+ the bidentate chelating configuration is preferred. In contrast, the bidentate bridging structure is energetically considerably more favorable for Mn4O4+. This binding motif stabilizes the 2D ring structure of the core of the Mn4O4+ cluster with respect to the 3D cubic geometry of the Mn4O4+ cluster core.


1 Introduction

Understanding the principles governing the oxidation of water and the rational design of novel water oxidation catalyst for sustainable (e.g., solar) energy storage are major challenges in the field of alternative energy production. Such initiatives, typically called “artificial photosynthesis”,1–5 are often inspired by the natural oxygen evolving complex (OEC) of photosystem II. The OEC itself continues to be a subject of intense research, and has been found to comprise an inorganic CaMn4O5 cluster which is stabilized and linked to the protein backbone by several amino acid residues.6 Some of these amino acids, such as glutamic acid or aspartic acid, bind to the cluster by bridging two Mn (or Ca–Mn) atoms with the deprotonated carboxyl group (bridging μ-oxo ligands).

The natural inorganic cluster of the OEC has inspired significant efforts aiming at the synthesis of biomimetic bi- and tetra-nuclear manganese oxide as well as calcium–manganese oxide complexes.7–17 For forming and stabilizing such clusters a large number of different, sometimes complex, O- and N-donor ligands, often forming five- or six-membered rings with the metal atoms, have been reported.8,15 In some cases, small carboxylic acid residues, such as acetate or pivalate ligands, have been used (see e.g. ref. 16,18,19).

To probe fundamental aspects of the water oxidation reaction, an increasing number of research groups study small metal oxide clusters isolated in the gas phase to model water activation and oxidation catalysts.20–26 In our project, we investigate gas-phase clusters consisting of elements also contained in the natural OEC, i.e., Mn, Ca, and O. In the first steps of this approach, we have prepared bare, non-ligated manganese27–30 and calcium–manganese oxide clusters31–33 of different size and composition, and studied their water activation and oxidation properties via reactivity measurements and via infrared spectroscopy. We now aim to increase the model system's complexity by attaching small carboxylic acids mimicking a ligand environment that could potentially stabilize the cluster structure and/or serve as proton/electron transfer media. Carboxylic acids represent the simplest mimics of the carboxyl group of amino acids and are thus suitable initial molecules for such studies.

In this contribution, we present IR-MPD spectra of a series of manganese oxide clusters containing up to nine Mn atoms ligated with one formic acid (HCOOH) molecule. We show that all clusters containing more than one Mn atom deprotonate formic acid. However, the binding of the resulting formate group is strongly cluster-size-dependent and can differ from the bridging μ-oxo configuration typically found for Mn2O2- and Mn4O4-based molecular catalysts. The consequences of cluster fluxionality32–34 and differences from the ligation with acetic acid35,36 are discussed.

2 Methods

2.1 Experimental methods

Cationic manganese oxide clusters were produced by pulsed laser ablation of a rotating manganese target using the second harmonic of a Nd:YAG laser. The ablation took place in a 3 mm diameter and 60 mm long growth channel in the presence of a short pulse of helium carrier gas seeded with 0.25% oxygen. To form cluster-HCOOH complexes, a mixture of about 0.2% HCOOH in helium was introduced via a second pulsed valve 50 mm downstream in a flow tube reactor (3 mm diameter, 45 mm long).

The reaction mixture was then expanded into vacuum forming a molecular beam before entering the intracavity region where it was irradiated by the IR laser beam of the Free Electron Laser for Intra Cavity Experiments (FELICE, 710–1800 cm−1; 10 μs pulse duration; spectral width set to approximately 0.4% FWHM of the central frequency) crossing it at an angle of 35°. A few μs after the interaction with FELICE, all clusters were extracted into a reflectron time-of-flight mass spectrometer by a set of pulsed high voltage plates and detected with a microchannel plate detector.37,38

To correct for long-term source fluctuations, the experiment was operated at twice the FELICE repetition rate, allowing for the recording of reference mass spectra in between successive FELICE pulses. Whenever FELICE was in resonance with an IR active vibrational mode of a given cluster, multiple IR photons were absorbed sequentially, leading to heating of the complex and finally to its fragmentation. The IR-MPD spectra shown in this contribution represent the depletion yield Y([small nu, Greek, tilde]) at wavenumber [small nu, Greek, tilde], calculated as image file: d3cp04035f-t1.tif, where I([small nu, Greek, tilde]) and I0 are the mass peak intensities with and without laser light, respectively, and P([small nu, Greek, tilde]) the macropulse energy. To reduce the IR fluence with which the complexes are irradiated, the whole instrument can be translated up to 300 mm from the focus position leading to a 30-fold reduction in intensity but increased overlap between the laser and molecular beam and thus an increased signal to noise ratio. All spectra presented in this work were recorded 290 mm from the focus. To further reduce the IR fluence and increase the spectral resolution, the overlap between laser and molecular beams can be purposely misaligned such that the molecular beam only observed the lower intensity part of the laser beam.

2.2 Theoretical methods

The atomic arrangements of the model complexes were explored by employing the Born–Oppenheimer spin density-functional theory molecular dynamics (BO-SDFT-MD)39 method with norm-conserving pseudopotentials43 and the generalized gradient approximation (GGA) for electronic exchange and correlations, using the Perdew, Burke, Ernzerhof (PBE)44 functional. We have used a plane-wave basis with a kinetic energy cutoff Ec = 62 Ry (843.6 eV), which shows convergent results. This corresponds to a real-space grid spacing of 0.4a0 (with a0 the Bohr radius); the real-space grid spacing for the density (and potential) was 0.133a0 corresponding to Ec = 555 Ry.

In the construction of the Mn pseudopotentials the valence electrons, 3d5 and 4s2, were characterized by core radii rc(s) = 2.35a0 and rc(d) = 2.35a0, with the s orbital treated as local. The Mn pseudopotential had a non-linear core correction with ∼10% of the [Ar] core charge. For the oxygen atom pseudopotential the valence 2s2 and 2p4 electrons were treated with rc(s) = rc(p) = 1.45a0, with the p orbital treated as local. The BO-SDF-MD method is particularly suitable for investigations of charged systems since it does not employ a supercell (i.e., no periodic replication of the ionic system is used). In all calculations, the dependence on spin multiplicity has been checked and the results that we report correspond to the spin multiplicities with the lowest energies. In the following we denote the difference between the number of spin-up and spin-down electrons by μ = NN.

The vibrational eigenmode frequencies and atomic displacement eigenvectors were determined in the harmonic approximation by direct diagonalization of the dynamical matrix, constructed by finite differencing of the forces at different points near the equilibrium geometry.45 From the DFT calculations for the displacements we also evaluate for each of the atoms the effective-charge tensor (the derivative of the dipole of the system with respect to the displacement of the atom), which together with the eigenmode vectors yields the oscillator strengths of the vibrational modes and the IR intensities.46 To correct for the general slight mismatch between the experimental and calculated vibrational frequencies, we applied an empirical scaling factor of 1.054 to all calculated vibrational modes. This factor is based on a comparison of the experimental and calculated bands I and II (isomers 2,2-A, 2,2-B, 2,4-A, 4,4-A) shown in Fig. 3–5.

3 Results and discussion

3.1 Deprotonation of formic acid on MnxOy+ (x ≥ 2)

In a previous study we have shown that small carboxylic acids such as formic, acetic, or propionic acid, bind to manganese cations via the oxygen atom of the carboxyl group resulting in a monodentate structure with the acid in both the trans- (cf. left structure of Fig. 1) and the cis-configuration (indicated by the arrow).47 In contrast, acetic acid was shown to dissociate (deprotonate) on some manganese oxide clusters forming a hydroxyl and an acetate group, the latter typically in a bidentate bridging or chelating structure (cf. middle and right structure in Fig. 1).36 In case of molecular binding, the C[double bond, length as m-dash]O stretch mode (νC[double bond, length as m-dash]O) of the carboxyl group is most characteristic (typically between 1690 cm−1 and 1560 cm−1) and even allows to distinguish between the trans- and cis-configuration.36,47 For deprotonated acids, νC[double bond, length as m-dash]O is no longer observed; instead, the symmetric and asymmetric stretch modes of the OCO group (νsOCO and νasOCO) become the diagnostic modes (νsOCO typically between 1480 cm−1 and 1380 cm−1; νasOCO typically between 1680 cm−1 and 1500 cm−1).40–42
image file: d3cp04035f-f1.tif
Fig. 1 Binding modes of small carboxylic acids to manganese oxide clusters and most characteristic vibrational modes.36,40–42 Purple and red balls indicate Mn and O atoms of the cluster. The left mode indicates the trans-configuration; the arrow indicates the position of the H atom in the cis-configuration.

Fig. 2 shows the IR-MPD spectra of series of MnxOy(HCOOH)+ (x = 1–9) complexes in the 1200–1800 cm−1 spectral region which is most diagnostic for the different adsorption configurations. From this figure it is evident that the spectra of the mono-manganese complexes Mn(HCOOH)+, MnO2(HCOOH)+ and MnO4(HCOOH)+ are similar, showing a broad, and sometimes structured, band around 1700 cm−1, whereas the spectra of all clusters with x ≥ 2 are dominated by two main bands centered around 1550 cm−1 and 1370 cm−1, respectively. It should be noted that the first band is slightly red-shifted to around 1510 cm−1 for Mn2O2(HCOOH)+ and Mn2O3(HCOOH)+.


image file: d3cp04035f-f2.tif
Fig. 2 IR-MPD spectra of MnxOy(HCOOH)+ complexes (labeled x, y to the right and left of the spectra, respectively). The blue and red dotted line are drawn to guide the eye. The black tick marks on the y-axis indicate the zero depletion yield for each spectrum.

The double peak of Mn(HCOOH)+ around 1700 cm−1 can be assigned to νC[double bond, length as m-dash]O of intact molecularly bound HCOOH in the trans and cis-configurations, respectively.47 The co-existence of these two isomers appears to be maintained upon co-adsorption of one and two oxygen atoms (the two peaks are not resolved for MnO2(HCOOH)+, but the peak width indicates the overlap of two modes). This strong band clearly disappears for all other cluster sizes indicating the deprotonation of the carboxyl group. To gain more insight into the binding motif (bidentate bridging vs. chelating) we will in the following analyze the spectra of selected cluster sizes in more detail.

3.2 Mn2O2(HCOOH)+: bidentate chelating and bridging structures

The top panel of Fig. 3 displays the IR-MPD spectrum of Mn2O2(HCOOH)+ in the whole investigated spectral region (710–1800 cm−1) together with the simulated vibrational spectra for two isomeric structures. These structures contain deprotonated formic acid with the formate group in a bidentate chelating (isomer 2,2-A; ΔE = −2.43 eV with respect to the reactants, μ = 9) and bidentate bridging (isomer 2,2-B, +0.07 eV, μ = 9) configuration, respectively. The IR-MPD spectrum shows two bands centered at 1517 cm−1 (labeled I) and 1360 cm−1 (labeled II), both potentially broadened on the red edge, and a broad structured band between 700 and 900 cm−1 (labeled III).
image file: d3cp04035f-f3.tif
Fig. 3 IR-MPD spectrum (top row) of Mn2O2(HCOOH)+ together with calculated vibrational spectra of two isomeric structures containing a formate group. For isomer 2,2-A the spectrum for the (ground-state) high spin (μ = 9, black curve), and for a higher-energy (by 0.09 eV compared to the low-spin ground-state) low spin isomer (μ = 1, red dashed curve) are shown in the second panel from the top. For isomer 2,2-B the spectrum of the (ground-state) high spin state (μ = 9) is shown in the third panel. The additional spectrum shown in the bottom panel represents the weighted sum of the spectra of isomer A and B (summed with a ratio of 1[thin space (1/6-em)]:[thin space (1/6-em)]2). Mn, O, C, and H atoms are depicted as purple, red, black, and white spheres.

The calculated vibrational spectrum of isomer 2,2-A exhibits a strong mode at 1499 cm−1 (νasOCO), two less intense modes at 1381 cm−1 (νsOCO) and 1296 cm−1 (CH wagging), a low intensity mode at 1051 cm−1 (CH out-of-plane bending) as well as two modes at 868 cm−1 (OH wagging) and 816 cm−1 (OCO bending). νasOCO seems a plausible candidate for band I, and the modes at 868/816 cm−1 might explain band III, although the former, despite the scaling factor of 1.054, still lies 18 cm−1 below the experimental band maximum. However, the low intensity modes at 1381/1296 cm−1 cannot satisfactorily account for band II. From our calculations we observe a rather small energy difference (0.09 eV) between the high spin state (ground state) and the lower spin (μ = 1) isomer. Furthermore, the structures of these two spin states are found to be very similar (imperceptible difference on the scale of the figure (see second panel of Fig. 3)). This leads to highly similar vibrational spectra for the two isomers (black and red curves) showing that the vibrational frequencies are rather insensitive to the spin state for these system.

For isomer 2,2-B (+0.07 eV for μ = 9 and +0.25 eV for μ = 1) the νasOCO stretch is predicted at 1529 cm−1 and νasOCO and the CH bending are shifting closer together. Furthermore, the OH wagging and the OCO bending vibrations red-shift to the edge of the investigated spectral region providing a better match with the maximum of the experimental band III. Thus, bands I, II, and III are better described by isomer 2,2-B, but that cannot completely account for the overall width of band III. So, neither isomer 2,2-A nor 2,2-B give a perfect match with the experimental spectrum. The situation improves when considering the co-existence of both isomers. The bottom panel of Fig. 3 displays the weighted sum of the spectra of isomer 2,2-A and 2,2-B (summed with a ratio of 1[thin space (1/6-em)]:[thin space (1/6-em)]2). It should be noted that this ratio is rather randomly chosen since the relative peak intensities do not seem to allow for a quantitative determination of the contribution of both isomers. However, the summed spectrum seems to explain all experimentally observed features, even the potential broadening of bands II and III. The only mode that is not observed in the IR-MPD spectrum is the low intensity mode predicted at 1051 cm−1 for isomer 2,2-A and 1057 cm−1 for isomer 2,2-B. The lack of this mode can be attributed to its very low intensity. Thus, we conclude that the IR-MPD spectrum of Mn2O2(HCOOH)+ is best described by the co-existence of isomers 2,2-A and 2,2-B, which is plausible given the small energy difference between them.

3.3 Mn2O4(HCOOH)+: bidentate chelating structure

As a second model system we have theoretically studied the complex Mn2O4(HCOOH)+ in more detail. The IR-MPD spectrum (Fig. 4; compared to Fig. 2 this spectrum is obtained at reduced IR macropulse fluence) shows five bands centered around 1538 cm−1 (labeled I), 1369 cm−1 (II), 1261 cm−1 (III), 1164 cm−1 (IV), and 810 cm−1 (V). The calculated vibrational spectrum of isomer 2,4-A (ΔE = −3.27 eV with respect to the reactants for μ = 7; +0.15 eV for μ = 1), with the formate in a bidentate chelating configuration, exhibits four intense modes at 1492 cm−1 (νasOCO), 1383 cm−1 (νsOCO), 1183 cm−1 (O–O stretch), and 802 cm−1 (a combined OH wagging and OCO bending motion) as well as a low intensity mode at 1287 cm−1 (CH wagging). These modes can account for all the experimentally observed bands, even for the low intensity band III. However, similar to Mn2O2(HCOOH)+, band I (νasOCO) for Mn2O4(HCOOH)+ is predicted too low in frequency by about 46 cm−1 despite the scaling factor of 1.054. Also, it could be argued that only the mode predicted at 802 cm−1 is not sufficient to explain the broad band V.
image file: d3cp04035f-f4.tif
Fig. 4 IR-MPD spectrum (top) of Mn2O4(HCOOH)+ together with calculated vibrational spectra of two isomeric structures containing a formate group. Mn, O, C, and H atoms are depicted as purple, red, black, and white spheres. The IR-MPD spectrum has been obtained with reduced IR macropulse fluence to achieve better resolution. The spectrum obtained at higher fluence is shown in Fig. 2.

A second isomer 2,4-B (μ = 7) with the formate in a bidentate bridging position is 0.36 eV higher in energy and is clearly not providing a better match with the IR-MPD spectrum. In particular, the blue-shift of the O–O stretch mode to 1229 cm−1 caused by the rotation of the O2 molecule out of the plane of the Mn2O2 cluster core breaks the agreement with band IV. Furthermore, νsOCO is considerably red-shifted to 1329 cm−1, no longer matching band II very well. Finally, bands III (weak) and V (very intense) are not reproduced by this isomer. Thus, we conclude, that the IR-MPD spectrum of Mn2O4(HCOOH)+ best matches that of isomer 2,4-A with the formate in a bidentate chelating configuration.

3.4 Stabilization of the 2D Mn4O4+ cluster core

As a last model system, we theoretically studied the tetra-manganese complex Mn4O4(HCOOH)+. We have previously predicted that the bare Mn4O4+ cluster has a 2D ring-like structure which can gradually change to a 3D cubic geometry upon adsorption of multiple water molecules29,30 due to the high fluxionality of this cluster.27,32 The cubic structure of the bare cluster was predicted to be 0.25 eV higher in energy. Here we aim to gain insight into the structural changes of the cluster upon dissociative adsorption of formic acid.

Fig. 5 shows the IR-MPD spectrum of Mn4O4(HCOOH)+ together with the calculated vibrational spectra of several isomers containing a ring-like and a cubic-like structure. The spectrum is similar to that for Mn2O2(HCOOH)+ (Fig. 3) and shows three intense bands centered around 1554 cm−1 (I), 1369 cm−1 (II), and 800 cm−1 (IV). In addition, a low intensity band is observed around 1046 cm−1 (III), close to where a band was theoretically predicted for Mn2O2(HCOOH)+ but not experimentally observed.


image file: d3cp04035f-f5.tif
Fig. 5 IR-MPD spectrum (top row) of Mn4O4(HCOOH)+ together with calculated vibrational spectra of several 2D and 3D isomeric structures. Mn, O, C, and H atoms are depicted as purple, red, black, and white spheres.

The lowest energy isomer containing deprotonated formic acid (isomer 4,4-A, μ = 19) found in our theoretical structure search is based on a ring-like Mn4O4+ cluster core. The resulting formate group bridges two Mn atoms and the hydrogen has migrated to the bridging oxygen atom opposite the formate. The bridge-binding of the formate induces the relocation of the oxygen atom initially bridging the two Mn atoms to the middle of the Mn4 square. In contrast, transfer of the hydrogen to the oxygen atom originally bridging the two Mn atoms where the formate has adsorbed (isomer 4,4-B, μ = 19) prevents the oxygen atom from being pushed inside the Mn4 ring. Instead, the formate group rotates out of the Mn4O4 plane allowing the hydroxyl to bridge the same Mn atoms, but opposite the formate. When the formate binds in a bidentate chelating configuration (isomer 4,4-C, μ = 19) the ring-like structure is also retained although it is considerably distorted. Isomers 4,4-B and 4,4-C are considerably higher in energy than isomer 4,4-A.

The vibrational spectrum of isomer 4,4-A shows two strong bands at 1588 cm−1 (νasOCO) and 1400/1386 cm−1 (νsOCO and CH bending) which could account for bands I and II. If so, theory overestimates the frequency of the OCO vibrations, which might indicate that the chosen scaling factor of 1.054 is too large for the Mn4O4 complex. In addition, a mode at 1040 cm−1 is predicted, which is in agreement with band III. The broader band IV is not so well matched by modes at 907 cm−1 and 751 cm−1. The vibrational spectrum of isomer 4,4-B is similar to that of isomer 4,4-A and could equally well explain bands I and III (with the frequency of band I even higher). However, the relative intensity of band II to band I seems to give this spectrum a poorer match with the experiment. On the other hand, the low frequency part may be a better match although the intensity of modes predicted below 800 cm−1 seems at odds with the rather low intensity band IV. Although this isomer is 1.19 eV higher in energy than isomer 4,4-A, both are thermodynamically accessible at −3.33 eV and −2.14 eV with respect to the reactants, respectively.

The same holds for isomer 4,4-C at 2.0 eV above isomer 4,4-A. Its spectrum contains a νasOCO that is red-shifted to 1530 cm−1, now slightly lower than band I. This discrepancy is similar to the case of Mn2O2(HCOOH)+ discussed above. Again, the two low intensity bands around 1400 cm−1 do not appear a perfect match for band II. Interestingly, the more intense modes predicted below 850 cm−1 do match the overall width of band IV, but we hesitate about the predicted strengths.

Finally, we briefly discuss the binding of formic acid to the cubic isomer of Mn4O4+. In this case, we also find that bridge-binding of the formate is energetically more favorable than a chelating geometry. Most remarkably, the lowest energy cubic isomer that we found (isomer 4,4-D, μ = 19) is 0.86 eV higher in energy than isomer 4,4-A, although the cubic structure of the bare Mn4O4+ is only 0.25 eV less stable than the ring structure. This shows that adsorption of one formic acid molecule in a bridging configuration considerably stabilizes the 2D ring structure. In contrast, binding of the formate in a bidentate chelating geometry is energetically more favorable for a 3D cubic cluster core than a 2D ring-shaped cluster (compare isomer 4,4-C with isomers 4,4-E,F,G, all μ = 19). Nevertheless, all isomers formed with a cubic cluster core are considerably less stable than isomer 4,4-A.

For all isomers 4,4-D,E,F,G, containing a cubic cluster core, νasOCO is, like for 4,4-C, slightly lower than for 4,4-A and 4,4-B. None of the calculated spectra provide an excellent match with the IR-MPD spectrum. For example, the low intensity modes around 1300–1400 cm−1 of isomers 4,4-D and 4,4-G cannot account for the intensity of band II, whereas the two more intense modes of isomers 4,4-E and 4,4-F in the same spectral range should probably have led to a splitting or at least a considerable broadening of mode II. Although they all provide a plausible match for band III, none has a convincing match for band IV, but so neither have any of the spectra calculated for 4, 4-A and 4,4-B, or 4,4-C.

Consequently, we base our assignment predominantly on the comparison for bands I-III and find the overlap with isomer 4,4-A the most plausible. We conclude that formic acid readily deprotonates on Mn4O4+ most likely binding formate in a bidentate bridging geometry on a ring-shaped Mn4O4+ cluster core (isomer 4,4-A).

3.5 Carboxylic acids as ligands for MnO based cluster catalysts

The IR-MPD spectra of MnxOy(HCOOH)+ revealed the deprotonation of the hydroxyl group for all clusters with x ≥ 2, whereas formic acid is absorbed as intact molecule on Mn+, MnO2+, and MnO4+. Similar results were previously found for acetic acid.36 In contrast, for Mn2O4+ and Mn3O5+ indications of an intact acetic acid were found, suggesting deprotonation is kinetically hindered for these cluster sizes. This difference between formic and acetic acid can likely be attributed to the former's higher acidity.

More important is, however, the binding motif of the resulting formate and acetate groups, respectively. For both Mn2O2+ and Mn2O4+ binding of formate and acetate36 in the bidentate chelating configuration is energetically more favorable which is in agreement with the IR-MPD spectra. Only formate can also bind to Mn2O2+ in a bidentate bridging structure. This indicates that bridge-binding of a carboxylic acid residue on a bare di-manganese oxide cluster core is generally rather unfavorable. This seems to conflict with the reported bridge-bound acetate (OAc) groups of molecules synthesized in the liquid phase. However, in all these cases the Mn atoms were co-coordinated to larger ligands, which can potentially block the chelating position (see e.g. structures 51 to 67 in ref. 8 or structures 295 and 297 in ref. 15).

The situation changes for the Mn4O4+ cluster. The lowest energy structure, to which we assign the IR-MPD spectrum, is a 2D ring-like structure, which is stabilized by bridge-binding of a formate group. Additionally, for a 3D cubic core the bridging structure is energetically more favorable than the chelating configuration. So far, only few cubic Mn4O4 clusters with carboxylic residues as stabilizing ligands have been reported; among these are Mn4O4 cores stabilized on vanadium-18 and tungsten-48 based polyoxometalates (POMs), respectively, capped with three acetate groups as well as a Mn4O4 cubic core with three bridge bound acetate groups besides other ligands.17 In all these examples most of the Mn atoms are either co-coordinated to the POM or other ligands, which might potentially inhibit chelating of the acetate. However, with our free cluster model system we show that bridge-binding of a first formate group is considerably more stable even without additional stabilizing ligands. Furthermore, while the Mn2O2+ and Mn2O4+ cluster cores appear to be rather unaffected by formate binding, it induces considerable structural changes to the Mn4O4+ core. The particular structural fluxionality of Mn4O4+ has previously been described with respect to water adsorption, deprotonation, and oxidation.29,32,33 Based on these observations, further structural changes can be expected upon adsorption of multiple ligands, which might eventually also lead to the stabilization of the cubic structure and a change of its activity in the water oxidation reaction. This will be subject of further investigations.

4 Conclusion

We have studied the interaction between manganese oxide clusters and formic acid with the use of IR-MPD spectroscopy in conjunction with DFT calculations. The IR-MPD spectra indicate binding of intact HCOOH to Mn+, MnO2+, and MnO4+, whereas the spectra of all other cluster sizes strongly indicate the deprotonation of formic acid. Supporting DFT calculations for three selected model systems, Mn2O2+, Mn2O4+, and Mn4O4+, show that depending on the cluster size the resulting formate group can bind differently. For Mn2O2+ formate binding in a bidentate bridging and chelating configuration are almost isoenergetic in agreement with the IR-MPD spectrum that can be best described by the co-existence of both isomers. Adding one oxygen molecule to the cluster (Mn2O4+) leads to the stabilization of the bidentate chelating configuration. In contrast, the tetra-manganese oxide cluster Mn4O4+ prefers to bind the formate group in a bidentate bridging configuration. In addition, this is found to lead to the stabilization of the 2D ring-like Mn4O4+ cluster core with respect to the 3D cubic geometry, alebite inducing considerable structural changes. Thus, although all MnxOy+ clusters with x ≥ 2 are able to deprotonate formic acid the binding of the formate appears to be strongly cluster-size dependent. This finding is important with respect to the identification of suitable ligands for MnO based cluster catalysts.

Conflicts of interest

There are no conflicts to declare.

Acknowledgements

S. M. L. and T. M. B. thank the Deutsche Forschungsgemeinschaft for financial support. We gratefully acknowledge the Nederlandse Organisatie voor Wetenschappelijk Onderzoek (NWO) for the support of the FELIX Laboratory. The research leading to these results has received funding from LASERLAB-EUROPE (grant agreement no. 654148, European Union's Horizon 2020 research and innovation programme). Computations were carried out at the Georgia Tech Center for Computational Materials Science, and the work of U. L. and R. N. B. was supported by grant FA9550-21-1-0198 from the Air Force Office of Scientific Research (AFOSR).

Notes and references

  1. D. G. Nocera, Acc. Chem. Res., 2012, 45, 767 CrossRef CAS.
  2. N. Cox, D. A. Pantazis, F. Neese and W. Lubitz, Interface Focus, 2015, 5, 20150009 CrossRef.
  3. J. Barber and P. D. Tran, J. R. Soc., Interface, 2013, 10, 20120984 CrossRef.
  4. S. Paul, F. Neese and D. A. Pantazis, Green Chem., 2017, 19, 2309 RSC.
  5. B. Zhang and L. Sun, Chem. Soc. Rev., 2019, 48, 2216 RSC.
  6. Y. Umena, K. Kawakami, J.-R. Shen and N. Kamiya, Nature, 2011, 473, 55 CrossRef CAS.
  7. M. Yagi and M. Kaneko, Chem. Rev., 2001, 101, 21 CrossRef CAS.
  8. S. Mukhopadhyay, S. K. Mandal, S. Bhaduri and W. H. Armstrong, Chem. Rev., 2004, 104, 3981 CrossRef CAS.
  9. R. Brimblecombe, G. C. Dismukes, G. F. Swiegers and L. Spiccia, Dalton Trans., 2009, 9374 RSC.
  10. G. C. Dismukes, R. Brimblecombe, G. A. N. Felton, R. S. Pryadun, J. E. Sheats, L. Spiccia and G. F. Swiegers, Acc. Chem. Res., 2009, 42, 1935 CrossRef CAS.
  11. M. Yagi, A. Syouji, S. Yamada, M. Komi, H. Yamazaki and S. Tajima, Photochem. Photobiol. Sci., 2009, 8, 139 CrossRef CAS.
  12. J. D. Blakemore, R. H. Crabtree and G. W. Brudvig, Chem. Rev., 2015, 115, 12974 CrossRef CAS.
  13. C. S. Mullins and V. L. Pecoraro, Coord. Chem. Rev., 2008, 252, 416 CrossRef CAS PubMed.
  14. M. Wiechen, H.-M. Berends and P. Kurz, Dalton Trans., 2012, 41, 21 RSC.
  15. M. D. Kärkäs, O. Verho, E. V. Johnston and B. Åkermark, Chem. Rev., 2014, 114, 11863 CrossRef.
  16. S. Mukherjee, J. A. Stull, J. Yano, T. C. Stamatatos, K. Pringouri, T. A. Stich, K. A. Abboud, R. D. Britt, V. K. Yachandra and G. Christou, Proc. Natl. Acad. Sci. U. S. A., 2012, 109, 2257 CrossRef CAS.
  17. J. S. Kanady, E. Y. Tsui, M. W. Day and T. Agapie, Science, 2011, 333, 733 CrossRef CAS.
  18. B. Schwarz, J. Forster, M. K. Goetz, D. Yücel, C. Berger, T. Jacob and C. Streb, Angew. Chem., Int. Ed., 2016, 55, 6329 CrossRef CAS.
  19. E. A. Karlsson, B.-L. Lee, T. Åkermark, E. V. Johnston, M. D. Kärkäs, J. Sun, Ö. Hansson, J.-E. Bäckvall and B. Åkermark, Angew. Chem., Int. Ed., 2011, 50, 11715 CrossRef CAS.
  20. M. R. Fagiani, X. Song, S. Debnath, S. Gewinner, W. Schöllkopf, K. R. Asmis, F. A. Bischoff, F. Müller and J. Sauer, J. Phys. Chem. Lett., 2017, 8, 1272 CrossRef CAS.
  21. T. Zhou, M. Abe, Y. Zhang, S. Kudoh and F. Mafuné, J. Phys. Chem. A, 2022, 126, 8218 CrossRef CAS.
  22. D. Masuzaki, Y. Takehashi and F. Mafuné, Bull. Chem. Soc. Jpn., 2018, 91, 1530 CrossRef CAS.
  23. C. Geng, J. Li, T. Weiske and H. Schwarz, J. Am. Chem. Soc., 2018, 140, 9275 CrossRef CAS.
  24. J. Li, S. Zhou, X.-N. Wu, S. Tang, M. Schlangen and H. Schwarz, Angew. Chem., Int. Ed., 2015, 54, 11861 CrossRef CAS.
  25. L.-S. Chen, Y.-Z. Liu, J.-J. Chen, S.-D. Wang, T.-M. Ma, X.-N. Li and S.-G. He, J. Phys. Chem. A, 2022, 126, 5294 CrossRef CAS.
  26. J.-J. Chen, X.-N. Li, Q.-Y. Liu, G.-P. Wei, Y. Yang, Z.-Y. Li and S.-G. He, J. Phys. Chem. Lett., 2021, 12, 8513 CrossRef CAS.
  27. S. M. Lang, I. Fleischer, T. M. Bernhardt, R. N. Barnett and U. Landman, Nano Lett., 2013, 13, 5549 CrossRef CAS.
  28. S. M. Lang, I. Fleischer, T. M. Bernhardt and R. N. Barnett, J. Phys. Chem. C, 2015, 119, 10881 CrossRef CAS.
  29. S. M. Lang, T. M. Bernhardt, D. M. Kiawi, J. M. Bakker, R. N. Barnett and U. Landman, Angew. Chem., Int. Ed., 2015, 127, 15328 CrossRef.
  30. S. M. Lang, T. M. Bernhardt, D. M. Kiawi, J. M. Bakker, R. N. Barnett and U. Landman, Phys. Chem. Chem. Phys., 2016, 18, 15727 RSC.
  31. S. Mauthe, I. Fleischer, T. M. Bernhardt, S. M. Lang, R. N. Barnett and U. Landman, Angew. Chem., Int. Ed., 2019, 131, 8592 CrossRef.
  32. S. M. Lang, N. Zimmermann, T. M. Bernhardt, R. N. Barnett, B. Yoon and U. Landman, J. Phys. Chem. Lett., 2021, 12, 5248 CrossRef CAS.
  33. S. M. Lang, I. Helzel, T. M. Bernhardt, R. N. Barnett and U. Landman, J. Am. Chem. Soc., 2022, 144, 15339 CrossRef CAS.
  34. X. Xing, X. Li, B. Yoon, U. Landman and J. H. Parks, Int. J. Mass Spectrom., 2015, 377, 393 CrossRef CAS.
  35. S. M. Lang and T. M. Bernhardt, Int. J. Mass Spectrom., 2018, 433, 7 CrossRef CAS.
  36. S. M. Lang, T. M. Bernhardt, J. M. Bakker, R. N. Barnett and U. Landman, J. Phys. Chem. A, 2021, 125, 4435 CrossRef CAS.
  37. J. M. Bakker, V. J. F. Lapoutre, B. Redlich, J. Oomens, B. G. Sartakov, A. Fielicke, G. von Helden, G. Meijer and A. F. G. van der Meer, J. Chem. Phys., 2010, 132, 074305 CrossRef.
  38. M. Haertelt, V. J. F. Lapoutre, J. M. Bakker, B. Redlich, D. J. Harding, A. Fielicke and G. Meijer, J. Phys. Chem. Lett., 2011, 2, 1720 CrossRef CAS.
  39. R. N. Barnett and U. Landman, Phys. Rev. B: Condens. Matter Mater. Phys., 1993, 48, 2081 CrossRef CAS.
  40. J. E. Tackett, Appl. Spectrosc., 1989, 43, 483 CrossRef CAS.
  41. G. B. Deacon and R. J. Phillips, Coord. Chem. Rev., 1980, 33, 227 CrossRef CAS.
  42. F. P. Rotzinger, J. M. Kesselman-Truttmann, S. J. Hug, V. Shklover and M. Grätzel, J. Phys. Chem. B, 2004, 108, 5004 CrossRef CAS.
  43. N. Troullier and J. L. Martins, Phys. Rev. B: Condens. Matter Mater. Phys., 1991, 43, 1993 CrossRef CAS.
  44. J. P. Perdew, K. Burke and M. Ernzerhof, Phys. Rev. Lett., 1996, 77, 3865 CrossRef CAS.
  45. D. Porezag and M. R. Pederson, Phys. Rev. B: Condens. Matter Mater. Phys., 1996, 54, 7830 CrossRef CAS.
  46. P. Giannozzi and S. Baroni, J. Chem. Phys., 1994, 100, 8537 CrossRef CAS.
  47. S. M. Lang, T. M. Bernhardt, J. M. Bakker, R. N. Barnett and U. Landman, J. Phys. Chem. A, 2020, 124, 4990 CrossRef CAS.
  48. X. Fang, M. Speldrich, H. Schilder, R. Cao, K. P. O'Halloran, C. L. Hill and P. Kögerler, Chem. Commun., 2010, 46, 2760 RSC.

This journal is © the Owner Societies 2023