Stabilization and activation of molecular oxygen at biomimetic tetrapyrroles on surfaces: from UHV to near-ambient pressure

Recent advances in the development of surface science methods have allowed bridging, at least partially, the pressure gap between the ultra-high vacuum environment and some applicative conditions. This step has been particularly critical for the characterization of heterogenous catalytic systems (solid–liquid, solid–gas interfaces) and, specifically, of the electronic, structural, and chemical properties of tetrapyrroles at surfaces when arranged in 2D networks. Within a biomimetic picture, in which 2D metalorganic frameworks are expected to model and reproduce in a tailored way the activity of their biochemical proteic counterparts, the fundamental investigation of the adsorption and activation of small ligands at the single-metal atom reaction sites has progressively gained increasing attention. Concerning oxygen, biology offers a variety of tetrapyrrole-based transport and reaction pockets, as e.g. in haemoglobin, myoglobin or cytochrome proteins. Binding and activation of O2 are accomplished thanks to complex charge transfer and spin realignment processes, sometimes requiring cooperative mechanisms. Within the framework of surface science at near-ambient pressure (towards and beyond the mbar regime), recent progress has unveiled novel and interesting properties of 2D metalorganic frameworks and heterostacks based on self-assembled tetrapyrroles, thus opening possible, effective applicative routes in the fields of light harvesting, heterogenous (electro-)catalysts, chemical sensing, and spintronics.


Introduction
Light harvesting and chemical conversion of energy vectors represent relevant technological issues that will deeply affect the future of our energy economy in the view of a sustainable scenario. The development of novel, nano-engineered materials would allow indeed the efficient and sustainable conversion and the storage of solar energy by exploiting several routes like e.g. water splitting and/or carbon dioxide activation and reduction to chemical energy vectors. Photosynthesis in Nature already provides a way to store sunlight in chemical forms by splitting water into oxygen and hydrogen. It's about a thermodynamically disfavoured process, with an energy cost of 2.46 eV. 1 The combustion of H 2 (or other hydrogen-containing chemical vectors) in fuel cells yields large amounts of free energy with only water and heat as side wastes, thus closing ideally the cycle. In the whole process, hydrogen and oxygen evolution reactions (HER and OER, respectively), and the oxygen reduction reaction (ORR) are key catalytic processes that actually determine the efficiency and selectivity of the conversion steps. 1 Thus, a strategic milestone to achieve sufficient insight to design and development of tailored materials to this purpose is a full comprehension and an atomic-level, detailed description of the catalytic activation mechanisms of stable, small molecules like oxygen. [2][3][4][5] In this sense, Nature has a lot to congurations of the metal ion, together with the organic coordination environment, determine the actual catalytic activity and selectivity of the system with respect to specic reactions. In this sense, a biomimetic strategy may consist in the synthesis of 2D, ordered heterostacks of porphyrin-based reaction centres supported by metal surfaces. The latter act both as geometric templates and as electron reservoirs. 10 As we will discuss in detail in the following, the planar conguration of the tetrapyrrolic structure exposes for ligation both sides of the metal ion. In the case of biochemical systems, in addition to the effect of the residues that directly bind to the porphyrin's or corrole's macrocycle, 1 the competition effects in the transcoordination sphere between ligands yield tuning of the reaction centre's electronic conguration. This process, known as molecular trans-effect, corresponds to the surface trans-effect in the case of the biomimetic approach based on 2D crystalline systems (Fig. 1c). 7,8,14,15 Indeed, the trans-effect describes in general the situation in which two ligands in trans position, binding to the same metal centre, compete for coordination. 14 In the case of supported metallorganic molecules, one of the two ligands is replaced by the supporting surface. Thus, the latter mimics the biochemical trans-coordination counterpart, allowing for orbital rehybridization and charge transfer to the adducts across the porphyrin plane. It is indeed charge transfer in general, and electron donation in the specic case of the oxygen molecule, that actually contribute to the binding and activation of the ligand for further reaction (Fig. 1d). 9,16

Oxygen stabilization and activation in Nature
The electronic ground state of O 2 is a stable, low-reactivity triplet, which is 0.93 eV lower in energy with respect to the singlet state with paired spin conguration. 17 In order to promote the molecule's activation for further reaction, biomimetic catalysts need to reproduce electron-donor sites. Here, the molecule can bind, 9 electron donation can take place, and reactive superoxo (O 2 À ) or peroxo (O 2 2À ) species are formed, depending on the system. 18,19 Similar charged metastable transition states were recently observed also at the surface of conventional ORR and OER electrocatalytic materials by means of in situ experimental approaches, conrming electron injection in O 2 as one of the necessary steps for its activation. 20,21 Surface reactive O 2 À , O IÀ , and O IIÀ species were identied by their associated O 1s core level and edge spectroscopic ngerprints. 20,22 There is a long-lasting debate in the literature on this point, even for the "simple" case of oxygen interaction with single crystal metal surfaces in UHV, 23 like for the Al(111) termination. In the latter case, competing oxygen-metal hybridization, spin selection rules, and charge transfer differentially contribute in the O 2 adsorption/activation and dissociation mechanisms. [24][25][26] Spin inversion is surely necessary for O 2 activation by biosystems, like in heme. This is a key point since, even if ORR processes can be largely exothermic, spin-forbidden transitions contribute with high activation barriers. Most transition metals offer (i) unpaired electrons, both in the ground state or in almost degenerate excited states, allowing reaction with triplet O 2 , and (ii) spin-orbit coupled congurations, thus providing a quantum mechanical pathway to promote spin inversion, 17 possibly exploiting diradical intermediates. 12 The overall oxygen reduction process may involve two or four electrons ( Fig. 1e) The former pathway yields H 2 O 2 , a good reactant or possible energy carrier, and is generally mediated by late transition metal centres. The latter path is instead more relevant from the biological point of view, being also energetically more favoured (3.47 eV), 12 and promoted by early transition metal complexes. 1 The origin of these differences is explained on the basis of the selectivity towards the formation of terminal metal-oxygen intermediates as a result of the heterolytic cleavage of the O 2 molecule. The d orbitals of late transition metal elements provide electrons that populate the metal-oxo antibonding states, thus making terminal oxo intermediates unstable. Nevertheless, in some systems, also late transition metals can catalyse four-electron O 2 reduction processes, as in the case of bimetallic reaction centres. Already the very rst one-electron step of oxygen reduction to the superoxide species is endoergonic, suggesting that a rigid, stepwise mechanism is rather unfavourable. Cooperative, multi-participated mechanisms are surely more effective, in which the tetrapyrrole centre plays several functions. This is for example the case of the oxygenation mechanism promoted by mono-oxygenase, where a twoelectron process takes place. One e À reduces the Fe(III) ferric site to the Fe(II) ferrous species and is subsequently transferred to the O 2 adduct that binds forming a ferric-superoxo intermediate. 9 A second e À couples with a proton, yielding a ferric-(hydro)peroxo intermediate, promoting the heterolytic scission of the O-O bond by protonation of the distal oxygen atom. The electron transport chain is therefore a key step: in the case of heme proteins, the iron porphyrin sites are involved in both electron transport and ORR. 1 In the specic case of Fe-O 2 binding in heme complexes, 27 a singlet electronic ground state is obtained upon reaction of triplet molecular oxygen with the high-spin Fe of the porphyrin (S ¼ 2). The conguration of the complex is mainly explained on the basis of the Weiss model: a low-spin (S ¼ 1/2) ferric centre is antiferromagnetically coupled to a doublet superoxide anion (O 2 À ), where the Fe-O 2 ligation is established by a s(Fe-O) bond involving Fe 3d z 2 and p*(O 2 ) orbitals, and by a weak p bonding between Fe d yz and p t *(O 2 ). 12 Not only Fe, but also Co plays relevant roles in the ORR. Cobalt-dioxygen ground states are doublets (S ¼ 1/2) where a Co(III) unpaired electron is donated to the O 2 moiety. A Co(II)-octaethylporphyrin supported on graphite was found to bind dioxygen (0.72 eV), as well as CoTPyP/Au(111) (0.85 eV) thanks to electron donation from the substrate, showing that a biomimetic approach can yield interesting results. 12,16 Indeed, recent developments have shown that metal-organic frameworks (MOFs) based on tetrapyrroles and non-heme M/N/C (M ¼ metal, N ¼ nitrogen, C ¼ carbon) catalysts may provide an optimal playground for oxygen activation. 28,29 However, beyond the local electronic conguration, many other aspects need to be considered in order to tailor the active sites, since transeffects, 14 distortion of the tetrapyrrolic plane, 13 electrostatic interactions, 12 cooperative effects, 30 and steric effects 31 play relevant roles.

Tetrapyrroles at the solid-liquid interface
Concerning biomimetic 2D systems based on tetrapyrroles, several efforts have been recently performed in order to characterize the oxygen interaction at the solid-liquid interface of electrolytic environments in situ and operando in real time. 32 Specically, the interest for the single-atom level detail is growing in order to achieve a thorough and fundamental insight. Scanning Tunnelling Microscopy (STM) is the main technique in this specic case, limited however to selected system choices and by environmental constraints. In this way, homolytic O 2 cleavage was observed on Mn porphyrins on Au(111) in Ar-saturated n-tetradecane, yielding reactive Mn(IV)] O species (Fig. 2a-c). 32 The reversible binding of O 2 (0.71 eV) with Co(II)OEP/HOPG molecules has been observed in a Langmuir uptake experiment in phenyloctane ( Fig. 2d and e). 33 Formation of the Fe-O 2 complex was observed for the FePc/ Au(111) system in 0.1 M HClO 4 as a transition state of the ORR to H 2 O 2 via the two-electron mechanism occurring at 50 mV potential ( Fig. 2f and g). 34 In the same environment, also

Specic focus of this review
Despite the remarkable insight achieved at the liquid/solid interface thanks to in situ atomic resolution imaging approaches, still, the inuence of the environment (very selected solvents and systems) and the lack of available complementary characterization techniques prevent a genuine focus on the fundamental aspects of the O 2 adsorption and activation mechanisms. On the other way, Ultra-High Vacuum (UHV) methods typical of surface science approaches can provide deeper insight, but severely suffer from the extreme pressure gap with respect to ambient, biological, or applicative conditions. 8 The recent development of spectroscopic techniques that allow at least a partial mitigation of the pressure gap without loss of insight opened to the possibility of investigating 2D systems up to Near-Ambient Pressure (NAP) conditions. These techniques include Near-Ambient Pressure X-ray Photoelectron Spectroscopy (NAP-XPS) and Infrared-Visible Sum-Frequency Generation spectroscopy (IR-Vis SFG 2. Solid-gas interfaces from UHV to near-ambient pressure: oxygen interaction with 2D metallorganic frameworks Within the framework of ultra-high vacuum surface science, the adsorption of ligands on tetrapyrroles generally requires cryogenic temperatures. This is ascribed to the small binding energy of the ligands to the metal ions, of the order of a few tenths of an eV. As an example, adsorption of carbon monoxide or nitric oxide onto the FePc/Au(111) monolayer is observed at 20 K, 51,52 while pyridine and ammonia bind at 80 K. 51 NH 3 was found to bind to NiTPP/Co(001) aer an uptake performed at 78 K. 53 However, nitrosyl complexes can be more stable in some cases and e.g. form in UHV on M(II)-tetraphenylporphyrins (M ¼ Fe, Co, Zn) on Ag(111) at 140 K, 14 and even at room temperature on CoTPP/O/Ni(001)/Cu(001) 54 and NiTPP/Cu(001). 55 Apart from these few cases, the population of the reactive metal ion centres of tetrapyrroles with ligands at room temperature, thus at applicative conditions, requires pressure regimes beyond the conventional surface science environment, generally in the few mbar range, in line with the biologic counterparts. Recently, the rst in situ experimental observation of the monocarbonylation of FePc/Ir(111) was obtained at room temperature in the mbar pressure range, 46 thus opening to the possibility of investigating these systems at close-to-ambient conditions with atomic-level insight. 8 Soon thereaer, the monocarbonylation of FePc/GR/ Ir(111) at 10 mbar unveiled interesting singlet-ssion mechanisms involved in the efficient absorption of visible light and its conversion in long-lived excitons. 44 Carbon dioxide could be stabilized on the same system by means of a biomimetic tuning of the charge transfer through the metallorganic heterostack. 7 A similar evolution in the fundamental characterization process could be observed in the very recent years also for the study of dioxygen adsorption.

Oxygen adsorption and activation in UHV
2.1.1 Oxygen adsorption on 2D frameworks at surfaces in UHV. Partial stabilization of the dioxygen adduct was observed at room temperature on Co(II)Pc/Ag(111) upon exposure of the layer to 1800 L of O 2 . 56 The combination of Scanning Tunnelling Microscopy (STM), Raman spectroscopy, and Density Functional Theory (DFT) approaches allowed the identication of end-on O 2 /CoPc/Ag(111) species, always coexisting with the dissociation products, namely O/CoPc/Ag(111). In the STM images, bright asymmetric protrusions appear upon dioxygen ligation to Co, associated with the O-O tilting, at variance with the O/CoPc molecules (Fig. 3a middle-top). The O-O internal stretching mode is identied by isotopic labelling from TERS data combined with ab initio calculations (Fig. 3a). The tilted, end-on dioxygen stable ligation geometry at the Co atom is attained thanks to the surface trans-effect for which, upon adsorption on Ag(111), the Co atom of the CoPc receives more charge from both the porphyrin ring and the metal surface (0.98 e À ) with respect to the gas phase tetrapyrrole case (0.77 e À ). As already anticipated, also in the specic case of the CoPc/Ag(111) system it is evident that stabilization of adsorbed O 2 requires charge donation, thus weakening the O-O bond and competing with dissociation through superoxo or peroxo intermediates, yielding the copresence of O/CoPc/Ag(111) species.
2.1.2 Oxygen activation and dissociation. From the computational point of view, a detailed description of the O 2 adsorption and possible decomposition at metal ions in tetrapyrroles is still an outstanding challenge, reecting in the choice of the xc approximation within the DFT framework, in particular for large system sizes. 57 As an example, for the case of O 2 ligation to CoTPP, both in the gas phase and on Au(111), it is shown that the results obtained by means of approaches based on DFT + U and GGA can be dramatically different (Fig. 3b). 57 In the latter case, stabilization of dioxygen is favoured and the dissociation barrier is of the order of 1 eV, while in the former case adsorption is energetically uphill and the O-O bond scission costs about 2 eV. Considering that the experimental evidence proves that, for the CoPc/Ag(111) system, stabilization and activation pathways may coexist (Fig. 3a), 56 this makes the picture quite puzzling. An even more complex example is represented by the interaction of O 2 with a Mg-TPA network on Cu(100), a 2D biomimetic model of the RuBisCO enzyme active site. 58 Interaction of dioxygen with the Mg 2+ centre at room temperature occurs in UHV already aer only 27 L O 2 , yielding the progressive disassembly and collapse of the network for larger exposures (Fig. 3c). A metastable dioxygen adduct binds to Mg 2+ (0.79 eV) undergoing dissociation and inducing oxidation of both Mg and the uncovered Cu(100) substrate. The same TPA molecules on Cu(100) can self-assemble in order to coordinate di-iron sites, forming an ordered 2D superstructure (Fig. 4a), 59 where O 2 readily dissociates at room temperature in UHV (Fig. 4b). The complex mechanism reveals the cooperative contribution of adjacent Fe sites and multiple O 2 molecules, with an activation barrier of only 0.74 eV (Fig. 4c). Also in this case, exposure to molecular oxygen progressively destroys the superstructure (Fig. 4b). Cooperative dissociation of O 2 is observed for Mn(II)TPP/Ag(111), where exposure of the Mn(II) TPP monolayer yields adjacent couples of Mn(II)OTPP molecules. 60 Interestingly, oxygen activation and dissociation was observed in UHV also on metal-free, N-doped graphene grown on Ir(111). The system was exposed to 5 Â 10 À5 mbar of molecular oxygen at 200 C for 30 min, inducing the formation of epoxy and ether species, accompanied by minority carbonyl species. 61 2.1.3 Oxygen at the interface. As discussed in the precious section, following exposure to O 2 , a few, selected metalorganic frameworks and tetrapyrroles can stabilize the dioxygen moiety and induce its dissociation also in UHV. However, while the O atoms can oxidize the metal ion, they may also interact with the underlying supporting surface. This is the case of the FePc/ Ag(110) system. 62   coordinate to on-top Ag sites. Upon exposure to 1700 L O 2 at RT, the Fe centres of the tetrapyrroles undergo strong dimming in their STM imaging appearance ( Fig. 5a and b). By combining experimental information and DFT modelling, it is demonstrated that oxygen intercalates at the tetrapyrrole-surface interface, forming a stable FePc-(h 2 -O 2 )-Ag coordination complex with a saddle-shape, distorted macrocycle (Fig. 5b). Interface oxygen can strongly inuence the oxidation state of the coordinated metal and can be thus exploited as a chemical switch to control the metal ion reaction site. This occurs for the NiTPP/Cu(100) system. 45,63 The strong metal-substrate interaction induces a sizeable charge transfer to the molecule, yielding population of the LUMOs up to the 3+ level and the formation of a Ni(I) ion, at variance to the Ni(II) oxidation state typical of the gas phase porphyrin. On the contrary, in the NiTPP/O/Cu(100) layer, oxygen decouples the molecule from the substrate and the Ni(II) oxidation state is preserved. The latter phenomena were observed at room temperature in UHV by means of a combination of spectroscopy, microscopy, and computational approaches ( Fig. 5c and d). Finally, the combined presence of a supporting and interacting surface and of an oxygen phase is known to lower the metalation barrier for adsorbed tetrapyrroles. In the biologic environment, free-based species undergo metalation in the homogenous liquid phase by reaction with dissolved metal salts. Enzymes catalyse the process by a multistep reaction pathway involving a signicant geometric distortion of the macrocycle. At surfaces, at least two common aspects are shared with the biologic counterpart. The so-called sitting atop complex plays a relevant role as an intermediate metalation step, 64,65 in which the metal centre already forms a coordination bond to the macrocycle, but the central H atoms are not yet released. Secondly, the molecule-surface interaction may induce geometric distortions of the tetrapyrrolic macrocycle, in analogy to the enzymatic process. 31,66 This is obtained thanks to the conformational adaptation of the tetrapyrrole to the local adsorption environment through the exibility of the porphyrin plane and the rotational degrees of freedom of the meso-groups. 67 All these steps are facilitated by the co-presence of oxygen, nally forming water upon reaction with the central H atoms in the tetrapyrrole pocket. 8,68,69 This is the case for the self-metalation reaction of 2H-TPP molecules adsorbed on Cu(111), where the pre-oxidation of the metal termination lowers the reaction temperature by 185 AE 15 K, down to room temperature (Fig. 5e). 69 Same reasoning applies to the Pd(100) substrate, where the temperature-induced self-metalation of 2H-TPP to PdTPP cannot be obtained up to 600 K in UHV, when competing desorption and decomposition of the tetrapyrroles set in. Instead, by depositing the 2H-TPP molecules on the oxygen-saturated Pd(100) surface, full metalation of the layer can be accomplished already at 430 K. 48

Oxygen adsorption and activation at NAP
By progressively increasing the oxygen pressure beyond the UHV limits and approaching the Near-Ambient Pressure (NAP) regime, new phenomena set in, as already evidenced for other small molecules like CO and CO 2 . 8 We will discuss here two specic cases, i.e. oxygen binding, stabilization, and activation at Co tetrapyrroles on Au(111), 16 and the oxygen-induced selfmetalation of 2H-TPP on Pd(100). 49 2.2.1 Binding of molecular oxygen at near-ambient pressure. A 2D metal-organic framework self-assembled at the Au(111) termination is able to mimic the O 2 stabilization and activation mechanisms that are typical of the biochemical environment of proteins and enzymes. 16 A CoTPyP/Au(111) layer binds indeed dioxygen at room temperature above 5 Â 10 À4 mbar, forming a covalent bond at the Co centre. Vibronic measurements performed in situ at O 2 saturation conditions reveal the development of an O-O stretching mode at 1276 cm À1 (Fig. 6a), associated with a ligand conguration laying between the oxo-and the superoxo-species (O 2 dÀ , d ¼ 0.22 e from DFT calculations for the specic case). The injection of charge into the ligand occurs via the surface trans-effect, from Au through the macrocycle, yielding both stabilization of the adduct and the weakening of its internal O-O bond. The O 2 -Co desorption energy is evaluated in a Langmuir uptake experiment (Fig. 6b) on the basis of a simple kinetic model. With the constrain of the strong assumption of a pre-exponential factor of 10 13 s À1 , a value of 0.85 AE 0.03 eV is obtained, compatibly with what obtained by ab initio calculations within the framework of DFT. Upon adsorption, O 2 develops a strong dipole moment (0.11 D) due to the differential charging of the two non-equivalent O atoms, associated with the end-on bonding geometry. As in the biological counterpart, the O 2 molecule sits on-top of the Co atom in an end-on conguration, and the molecular axis is tilted by 118 ( Fig. 6c and d). There is evidence of a Co(3d z 2)-O(2p) overlap yielding covalent bonding at À7.9, À1.0, and À0.4 eV below the Fermi level, with a strong contribution associated with the charge transfer from Au through the macrocycle plane (trans-effect) (Fig. 6d) and the formation of a shared Au-Co-O 2 peak in the DOS at À7.9 eV. By comparing geometric and electronic structures, it is found that the Co-O 2 bond consists of a dative s(Co-O) interaction between vacant 3d z 2 and p*(O 2 ) and a weak p interaction between 3d yz and p t *. The tilting originates from the geometry of the p*(O 2 ), which, in order to vertically align with the Co 3d z 2, requires the bending of the molecular axis (Fig. 6c, DOS isodensity surface at 0.4 eV below the Fermi level). Interestingly, it can be noticed that also one of the N atoms of the macrocycle participates directly in the oxygen bonding. The total O 2 magnetization is lowered from S ¼ 1 by 15% upon adsorption at the Co site. This is a key point in the activation process of O 2 : indeed, its ground state is a triplet, poorly reactive with other small, common molecules that are instead in a singlet state. As described in Section 1.2, activation involves the formation of superoxo or peroxo species, undergoing the 2e À or the 4e À reduction processes, depending on the system. Spin inversion is necessary since spin-forbidden transitions contribute with high activation barriers. Most transition metals offer unpaired electrons, both in the ground state or in almost degenerate excited states, allowing reaction with triplet O 2 , and spin-orbit coupled congurations, thus providing a quantum mechanical pathway to promote spin inversion, possibly exploiting di-radical intermediates for reaction. 2.2.2 Molecular oxygen inducing self-metalation at nearambient pressure. As discussed above, the ordered 2H-TPP/ Pd(100) monolayer can undergo self-metalation in UHV when atomic oxygen is pre-adsorbed at the interface and by heating to 430 K. 48 Oxygen atoms favour the formation of the sitting-atop complex by lowering the metalation barrier, so that the reaction can be thermally activated below the critical temperature of competing 2H-TPP desorption and decomposition processes. However, it has been recently proven that by exposing the 2H-TPP/Pd(100) monolayer to 1 mbar O 2 the self-metalation reaction to PdTPP can be induced already at room temperature. 49 Evidence has been obtained by means of in situ IR-Vis SFG measurements (Fig. 7a). The intensity ratio of the vibronic resonances at 1327 and 1369 cm À1 was observed to be a good metalation marker when the 2H-TPP monolayer is deposited on the oxygen pre-covered Pd(100) termination and subsequently annealed to 410 K, thus inducing metalation and water production. By extending the investigation beyond UHV conditions, it is found that in 1 mbar O 2 the feature at 1369 cm À1 undergoes a blue shi of 6 cm À1 and almost disappears (red deconvolution and dashed line in Fig. 7a). Three modes contribute to the resonance: the C b -C a -NH asymmetric stretching, the H rocking of the macrocycle pyrrolic moieties, and the C a -N-C a 0 stretching. On the contrary, the lowest energy feature (1327 cm À1 ) mainly originates from the phenyl H rocking modes. The spectroscopic information suggests therefore a strong local deformation of the C-N bonds and the loss of the central H atoms, compatibly with the 2H-TPP metalation to PdTPP. NAP-XPS data of the N 1s core level conrm the picture. In the pristine layer in UHV, before exposure to oxygen, two symmetric N 1s spectral components of equal intensity can be resolved at binding energies of 399.8 and 397.8 eV, associated with the pyrrolic (N-H) and iminic N, respectively, of the 2H-TPP molecules. By exposing the layer to oxygen at room temperature, a dramatic evolution of the spectral lineshape is observed (Fig. 7b). Two intense N 1s photoemission peaks grow at 397.5 and 398.2 eV at the expenses of the former doublet, indicating full metalation and, additionally, allowing distinction between intact PdTPP porphyrins and PdTPP whose phenyl groups have undergone partial dehydrogenation due to reaction with adsorbed oxygen. By means of a simple kinetic model, the analysis of the dynamic uptake data of both IR-Vis SFG and NAP-XPS allow the quantitative evaluation of the self-metalation activation barrier (0.4 eV). Within this picture, the impinging O 2 molecules undergo dissociation at the centre of the tetrapyrroles, contributing at the same time to the formation of the metalation transition state (rst O atom, sitting-atop complex) and to the oxidation (second O atom) of the underlying Pd surface.

Summary and perspectives
I have reviewed very recent experimental steps towards the characterization of the interaction mechanisms of dioxygen with 2D heterostructures based on tetrapyrroles or MOFs in general. Important pieces of fundamental information have already been understood. Indeed, the cooperative contribution of several effects is necessary in order to stabilize and activate O 2 at single metal atom sites. These include (i) the rehybridization of the O 2 -metal electronic states, (ii) the trans-effect yielding charge transfer to the ligand through and from the macrocycle, and (iii) the spin reconguration to rearrange the stable oxygen triplet. Biomimetic systems can be designed and tuned to achieve all these steps, but still there is something missing with respect to the biologic counterpart. Specically, enzymatic and, more generally, proteins sites available for O 2 binding and/or reaction are placed within 3D pockets. The 3D environment allows for (i) the control of the in and out ux of both reactants and products, 70 (ii) the induction of steric/ geometric effects, 31,71,72 and (iii) the distal coordination of the ligands to functional terminations of the protein backbone. 73 Thus, the biomimetic modelling of the functions and properties connected with the 3D nature of the binding pockets is still far from reach in the context of 2D heterostructures. Steps (i) to (iii) may represent the next challenges in the perspective of the development of novel 2D materials based on tetrapyrroles with applicative potential in the elds of heterogeneous (electro)catalysis, light harvesting, electronics and spintronics.

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