Understanding Heteroatom-Mediated Metal–Support Interactions in Functionalized Carbons: A Perspective Review

Featured Application: Surface engineering for catalysis and energy. Abstract: Carbon-based materials show unique chemicophysical properties, and they have been successfully used in many catalytic processes, including the production of chemicals and energy. The introduction of heteroatoms (N, B, P, S) alters the electronic properties, often increasing the reactivity of the surface of nanocarbons. The functional groups on the carbons have been reported to be effective for anchoring metal nanoparticles. Although the interaction between functional groups and metal has been studied by various characterization techniques, theoretical models, and catalytic results, the role and nature of heteroatoms is still an object of discussion. The aim of this review is to elucidate the metal–heteroatoms interaction, providing an overview of the main experimental and theoretical outcomes about heteroatom-mediated metal–support interactions. Selected studies showing the effect of heteroatom–metal interaction in the liquid-phase alcohol oxidation will be also presented.


Introduction
Since their discovery at the end of the 20th century, nanostructured carbon allotropes such as fullerenes [1,2], carbon nanotubes [1,3], carbon nanofibers [4], and graphene [5] have attracted the interest of several researchers from disparate fields. Nowadays, the family of low-dimensional carbon allotropes includes numerous members [6], supporting a wide range of applications from electronics [7] to medicine [8], from optics [9] to environment protection [10], from energy storage [11][12][13][14] to catalysis [15]. While the synthesis strategies as well as the key structural features have already been extensively investigated, the focus has now shifted toward tuning the textural and electronic properties of these materials to optimize their performances in the pursuit of technological progress and manufacturing efficiency. Due to the nanometric sizes, the behavior of low-dimensional materials is dominated by surface phenomena; therefore, the design of carbon nanomaterials with tailored properties is essentially a matter of surface engineering. A common approach in surface engineering is functionalization, aiming to introduce a functional component that is able to confer novel properties to the material [16,17]. In most cases, the surface functionalization results in the incorporation of heteroatoms (commonly boron, nitrogen, oxygen, phosphorus, and sulfur) in the carbon framework [18][19][20][21]. Functional groups can be attached on the surface of carbon materials conclude with a collection of examples reporting the role of heteroatoms in enhancing metal-support interactions in catalysis.

Basic Concepts for Understanding Heteroatom Mediated Metal-Support Interactions
The deposition of metal nanoparticles or clusters onto the surface of carbon support generates an interface accompanied by a local electrostatic potential perturbation. The entity of this perturbation reflects the strength of the interactions between metal and support atoms at the interface, and it is strictly related to the electronic structures of these two actors [51]. The electronic structure of solids is usually described in terms of electronic bands. The nature and width of the electronic bands can be evaluated considering the density of states (DOS), which are defined as the number of allowed electron states per volume at a given energy. Two parameters related to DOS are of paramount importance for understanding metal-support interactions: the Fermi level and the work function. The Fermi level (EF) is an energy state having a 50% probability of being filled, and it coincides with the electrochemical potential of electrons in the solid. The position of the Fermi level with respect to the D-band is an indicator of the D unoccupied density. The work function, Φ, is the minimum energy that must be furnished to the metal system to extract an electron from the Fermi level, EF, and transport it at infinite distance, when the metal surface is electrically neutral. According to the metal-semiconductor boundary-layer theory, when two phases are contacted, the thermodynamic equilibrium is attained. The equilibrium condition requires that the electrochemical potential and then the Fermi energy level of the electrons of the two solids is the same at the interface. The equalization of Fermi levels implies charge transfers across the interface from one material to the other, whose direction can be predicted from the analysis of the electronic structures. In cases where the work function of the electrons of the metal is larger than in the support, electrons will migrate from the support to the metal, until the Fermi level at the interface is equilibrated, and vice versa. This flow of charge sets up a partial positive charge on the support in the proximity of the interface, and an induced negative charge on the surface of the metal. Moreover, in the presence of the charge at the interface, the rehybridization of electron orbitals occurs, leading to changes in the D-band population and in the density of states at the Fermi level. These phenomena are at the basis of the metal-support interactions of the electronic type.
The introduction of heteroatoms in the support framework causes electron relocalization and local variations in the density of states of the support (LDOS, local density of states), thus producing an increase or decrease of the work function. This modification unavoidably influences the charge

Basic Concepts for Understanding Heteroatom Mediated Metal-Support Interactions
The deposition of metal nanoparticles or clusters onto the surface of carbon support generates an interface accompanied by a local electrostatic potential perturbation. The entity of this perturbation reflects the strength of the interactions between metal and support atoms at the interface, and it is strictly related to the electronic structures of these two actors [51]. The electronic structure of solids is usually described in terms of electronic bands. The nature and width of the electronic bands can be evaluated considering the density of states (DOS), which are defined as the number of allowed electron states per volume at a given energy. Two parameters related to DOS are of paramount importance for understanding metal-support interactions: the Fermi level and the work function. The Fermi level (E F ) is an energy state having a 50% probability of being filled, and it coincides with the electrochemical potential of electrons in the solid. The position of the Fermi level with respect to the D-band is an indicator of the D unoccupied density. The work function, Φ, is the minimum energy that must be furnished to the metal system to extract an electron from the Fermi level, E F , and transport it at infinite distance, when the metal surface is electrically neutral. According to the metal-semiconductor boundary-layer theory, when two phases are contacted, the thermodynamic equilibrium is attained. The equilibrium condition requires that the electrochemical potential and then the Fermi energy level of the electrons of the two solids is the same at the interface. The equalization of Fermi levels implies charge transfers across the interface from one material to the other, whose direction can be predicted from the analysis of the electronic structures. In cases where the work function of the electrons of the metal is larger than in the support, electrons will migrate from the support to the metal, until the Fermi level at the interface is equilibrated, and vice versa. This flow of charge sets up a partial positive charge on the support in the proximity of the interface, and an induced negative charge on the surface of the metal. Moreover, in the presence of the charge at the interface, the rehybridization of electron orbitals occurs, leading to changes in the D-band population and in the density of states at the Fermi level. These phenomena are at the basis of the metal-support interactions of the electronic type.
The introduction of heteroatoms in the support framework causes electron relocalization and local variations in the density of states of the support (LDOS, local density of states), thus producing an increase or decrease of the work function. This modification unavoidably influences the charge transfer, the hybridization between the electron states of the metal and those of the support, and definitively the metal-support interactions ( Figure 2 transfer, the hybridization between the electron states of the metal and those of the support, and definitively the metal-support interactions ( Figure 2). Theoretical simulations are a valuable technique to explore the phenomena occurring at metalsupport interface. Binding energy, atomic distances, and relaxed geometries are useful descriptors can be derived from density functional theory (DFT). Molecular dynamics simulations based on a generic force field have been implemented to obtain information on the mobility of adsorbed metal species. Molecular orbital analysis provides a useful indication of the electron density and charge distribution.

Metal Single-Atom Adsorption on Doped Carbon Surfaces
Among the several functionalizations, the substitution of carbon atoms with nitrogen or boron has received great attention. In the case of substitutional B doping, a carbon atom is replaced by an atom with one less electron; consequently, sharp and localized states below the Fermi level (i.e., in the valence band) are generated. The electronic structure modifications induced by nitrogen doping are more complex, and they depend on the possible bonding configurations for N in graphitic networks, as detailed below.
Acharya et al. [57] provided a deep insight into the adsorption of Pt and Ru atoms onto isolated B and N-doped graphene surfaces. C18, C16B2, and C16N2 were selected as the carbon structural models for the adsorption of metal species respectively on pristine, B-doped and N-doped graphene. In addition, the co-presence of boron and nitrogen (co-doping) with different connections (heteroatom-carbon or heteroatom-heteroatom bonds) was also investigated. According to DFT calculations, the incorporation of boron atoms in the graphene lattice produces a lowering of the Theoretical simulations are a valuable technique to explore the phenomena occurring at metal-support interface. Binding energy, atomic distances, and relaxed geometries are useful descriptors can be derived from density functional theory (DFT). Molecular dynamics simulations based on a generic force field have been implemented to obtain information on the mobility of adsorbed metal species. Molecular orbital analysis provides a useful indication of the electron density and charge distribution.

Metal Single-Atom Adsorption on Doped Carbon Surfaces
Among the several functionalizations, the substitution of carbon atoms with nitrogen or boron has received great attention. In the case of substitutional B doping, a carbon atom is replaced by an atom with one less electron; consequently, sharp and localized states below the Fermi level (i.e., in the valence band) are generated. The electronic structure modifications induced by nitrogen doping are more complex, and they depend on the possible bonding configurations for N in graphitic networks, as detailed below.
Acharya et al. [57] provided a deep insight into the adsorption of Pt and Ru atoms onto isolated B and N-doped graphene surfaces. C 18 , C 16 B 2 , and C 16 N 2 were selected as the carbon structural models for the adsorption of metal species respectively on pristine, B-doped and N-doped graphene. In addition, the co-presence of boron and nitrogen (co-doping) with different connections (heteroatom-carbon or heteroatom-heteroatom bonds) was also investigated. According to DFT calculations, the incorporation of boron atoms in the graphene lattice produces a lowering of the Fermi energy (from −2.25 eV to −3.43 eV) and a consequent increase in the work function from 4.31 eV to 5.43 eV. Looking at the projected density of states (PDOS) for the p-orbital (Figure 3), a lower Fermi energy implies an augmented number of states (bonding and antibonding states) close to the Fermi energy and the appearance of more antibonding empty states above the Fermi level, which can help in minimizing the antibonding repulsions during the adsorption of metal species on the surface. In other words, a strong hybridization between the platinum d-orbitals and boron p-orbitals occurs. The essential consequence is the strengthening of the interaction between the B-doped surface and metal atoms. Fermi energy (from −2.25 eV to −3.43 eV) and a consequent increase in the work function from 4.31 eV to 5.43 eV. Looking at the projected density of states (PDOS) for the p-orbital ( Figure 3), a lower Fermi energy implies an augmented number of states (bonding and antibonding states) close to the Fermi energy and the appearance of more antibonding empty states above the Fermi level, which can help in minimizing the antibonding repulsions during the adsorption of metal species on the surface. In other words, a strong hybridization between the platinum d-orbitals and boron p-orbitals occurs. The essential consequence is the strengthening of the interaction between the B-doped surface and metal atoms. Computed values for binding energy of both Pt and Ru atoms confirmed this hypothesis, since the highest values are associated with the adsorption on boron-modified surfaces. In the presence of nitrogen atoms, the increased Fermi energy (−1.27 eV) and the lower work function (3.37 eV) lead to the presence of fewer antibonding empty states close to the Fermi energy, and thus to a weakening of the interaction between metal and the nitrogen atoms. However, although the hybridization between the nitrogen p-orbital and platinum d-orbital is very weak, the presence of nitrogen atoms can result in an enhancement of the interaction with Pt compared to undoped surfaces. The promotion effect of nitrogen can be ascribed to the activation of the neighboring carbon atoms, whose reciprocal π-π interactions are weakened by the lower electron delocalization induced by the Computed values for binding energy of both Pt and Ru atoms confirmed this hypothesis, since the highest values are associated with the adsorption on boron-modified surfaces. In the presence of nitrogen atoms, the increased Fermi energy (−1.27 eV) and the lower work function (3.37 eV) lead to the presence of fewer antibonding empty states close to the Fermi energy, and thus to a weakening of the interaction between metal and the nitrogen atoms. However, although the hybridization between the nitrogen p-orbital and platinum d-orbital is very weak, the presence of nitrogen atoms can result in an enhancement of the interaction with Pt compared to undoped surfaces. The promotion effect of nitrogen can be ascribed to the activation of the neighboring carbon atoms, whose reciprocal π-π interactions are weakened by the lower electron delocalization induced by the close proximity of nitrogen atoms with large electron affinity [58]. An analogous effect can be envisaged in defective carbon surfaces, where the presence of curvature, vacancies, or edges produces a disruption of the delocalized electronic structure of the graphene surface. Consequently, the actual interaction energy is often the balance between electronic effects deriving from both doping and local geometric strain. As an example, comparing the Pt adsorption energy on doped single-walled carbon nanotubes B(N)-SWCNT(6,6) and B(N)-SWCNT (8,0), in the former system, a more moderate enhancement was unveiled compared to the SWCNT(8,0), which was probably due to the reduced curvature in SWCNT (6,6) [59].
In the co-doped surfaces, the effect of nitrogen and boron doping are mixed; therefore, the interaction of Pt and Ru results stronger close to boron atoms and weaker close to nitrogen atoms.
However, it is worth demonstrating that the coordination of nitrogen species in the carbon framework can determine the interaction of metal atoms with nitrogen-doped supports. Indeed, when a nitrogen atom replaces a carbon atom of graphitic plane, it will be bonded with three carbon neighbors, thus generating a so-called graphitic nitrogen: a sp 2 nitrogen species. Nevertheless, graphitic nitrogen is not the unique species that can be created in an N-doped surface. The substitution of a carbon atom with nitrogen can occur according to almost three primary structural models differing in their nitrogen bonding configuration ( Figure 4): pyridine-like (N pyr ), pyrrole-like (N H ), and graphitic (N gr ) [21]. Pyridinic nitrogen refers to sp 2 two-fold coordinated nitrogen atoms contributing one p electron to the π-system. It could be incorporated as a vacancy (N v ) or as atom on an edge site (N e ). Nitrogen atom in pyrrole-like configuration is an sp 3 coordinated (usually in a five-membered ring) species sharing two p electrons in the π system. Doping with pyridine-like, pyrrole-like, and graphitic nitrogen species impacts differently on the electronic band structure of graphitic carbon materials, having unavoidable implications for the interactions with adsorbates. The adsorption of single atoms of Pt-group metals on graphitic surfaces containing nitrogen species with different coordination has been modeled by DFT. The results from theoretical modeling unanimously identify the following ascending order of the binding energies of noble metals: N(graphitic) > N(pyrrolic) > N(pyridinic) [60][61][62]. close proximity of nitrogen atoms with large electron affinity [58]. An analogous effect can be envisaged in defective carbon surfaces, where the presence of curvature, vacancies, or edges produces a disruption of the delocalized electronic structure of the graphene surface. Consequently, the actual interaction energy is often the balance between electronic effects deriving from both doping and local geometric strain. As an example, comparing the Pt adsorption energy on doped singlewalled carbon nanotubes B(N)-SWCNT(6,6) and B(N)-SWCNT(8,0), in the former system, a more moderate enhancement was unveiled compared to the SWCNT(8,0), which was probably due to the reduced curvature in SWCNT(6,6) [59].
In the co-doped surfaces, the effect of nitrogen and boron doping are mixed; therefore, the interaction of Pt and Ru results stronger close to boron atoms and weaker close to nitrogen atoms.
However, it is worth demonstrating that the coordination of nitrogen species in the carbon framework can determine the interaction of metal atoms with nitrogen-doped supports. Indeed, when a nitrogen atom replaces a carbon atom of graphitic plane, it will be bonded with three carbon neighbors, thus generating a so-called graphitic nitrogen: a sp 2 nitrogen species. Nevertheless, graphitic nitrogen is not the unique species that can be created in an N-doped surface. The substitution of a carbon atom with nitrogen can occur according to almost three primary structural models differing in their nitrogen bonding configuration ( Figure 4): pyridine-like (Npyr), pyrrole-like (NH), and graphitic (Ngr) [21]. Pyridinic nitrogen refers to sp 2 two-fold coordinated nitrogen atoms contributing one p electron to the π-system. It could be incorporated as a vacancy (Nv) or as atom on an edge site (Ne). Nitrogen atom in pyrrole-like configuration is an sp 3 coordinated (usually in a five-membered ring) species sharing two p electrons in the π system. Doping with pyridine-like, pyrrole-like, and graphitic nitrogen species impacts differently on the electronic band structure of graphitic carbon materials, having unavoidable implications for the interactions with adsorbates. The adsorption of single atoms of Pt-group metals on graphitic surfaces containing nitrogen species with different coordination has been modeled by DFT. The results from theoretical modeling unanimously identify the following ascending order of the binding energies of noble metals: N(graphitic) > N(pyrrolic) > N(pyridinic) [60][61][62]. Bulushev et al. [63] modeled the adsorption of a Pd atom on four different nitrogen forms: graphitic (Ngr), pyrrolic (NH), pyridinic at the monovacancy boundary (Nv), and pyridinic at the edge (Ne) nitrogen atoms, probing different adsorption sites in the graphene fragment.
The obtained relaxed structures and binding energies revealed that the Pd atom preferentially interacts with pyridinic nitrogen. A similar behavior was observed with other Pt-group metals Bulushev et al. [63] modeled the adsorption of a Pd atom on four different nitrogen forms: graphitic (N gr ), pyrrolic (N H ), pyridinic at the monovacancy boundary (N v ), and pyridinic at the edge (N e ) nitrogen atoms, probing different adsorption sites in the graphene fragment.
The obtained relaxed structures and binding energies revealed that the Pd atom preferentially interacts with pyridinic nitrogen. A similar behavior was observed with other Pt-group metals (Pt, Ru). In general, all Pt-group metal atoms strongly interact with N pyr sites, although the interaction with the edge N pyr sites is significantly weaker than the one with the vacancy (mono or di) N pyr atoms. This evidence can be rationalized considering the number and the stability of the bonds generated upon the interaction between metal and vacancy N pyr atoms. In particular, the highest affinity for divacancy N pyr atoms can be related to the formation of four stable bonds among the central metal atom and the adjacent heteroatom and carbon atoms. The strong interaction between metal atoms and monovacancy N pyr atoms is associated with the breaking of a C-C bond in the pentagon ring, the formation of three new bonds, and the concomitant creation of three more stable six-member ring systems. Differently, just two and one new bonds are formed when metal atoms are adsorbed on N pyr atoms on the armchair or zigzag graphene edge, respectively. As expected, the interaction with graphitic nitrogen is not energetically favorable. Pd and Pt atoms do not bind to this site, whilst they prefer to migrate to the neighboring carbon sites and locate above the center of the C-C bond, forming a bridge with two carbon atoms. Ru is the only one that is able to interact directly with the N gr nitrogen, forming a bond with this atom and two neighboring carbon atoms. The sequence for the binding energy of metal atoms with a particular site was Ru > Pt > Pd, which is in agreement with Acharya's results [57]. Conflicting conclusions were deduced regarding the influence of the delocalized π-electron system of the graphene surface on the electron density redistribution in the case of Pt and Pd atoms. The comparison of the electron density surfaces for N-doped graphene fragments before and after the interaction with the Pd atom suggests that the interaction with Pd induces a more positive charging of N atoms compared to the initial states. Otherwise, looking at natural bond orbital (NBO) analysis, the interaction of the Pt atom with N pyr sites imparts a partial positive charge to the metal center (except for the multiple interactions with two N pyr atoms at the zigzag edge). Conversely, electron density is transferred from the support to Pt, when the latter is bonded close to N gr .
It is obvious that also the electronic structure of the metal plays a central role. According to Hammer and Nørskov [64], the interaction of the transition metals tends to increase from right to left, while it decreases as one progresses down (from top to bottom) a group on the periodic table. The number of D electrons, and then, the position of D-band center related to the Fermi energy, play a crucial role.
This appears more evident when the study is extended to other transition metal atoms. An et al. [59] used density functional theory (DFT) calculations to investigate the adsorption of 11 transition metal atoms (i.e., Au, Pt, Ru, Pd, Ag, Fe, Co, Ni, Cu, W, and Ti) on the sidewalls of boron and nitrogen-doped single-wall carbon nanotubes (B-SWCNTs and N-SWCNTs). In general, the highest binding energies were obtained for Sc (3d 1 4s 2 ), Ti (3d 2 4s 2 ), and Fe(3d 6 4s 2 ), which was likely due to their capability to easily lose one or two electrons for reaching the full or half-filled steady states [65]. Again, an enhanced binding energy was observed for boron-doped surfaces, whereas N-doping produces a slighter increase in the interaction energies. In addition, a quenching to zero of the spin-magnetic moment (µ) of the individual metal atoms was predicted for almost all of the species (except Fe and W) upon their bonding on doped carbon nanotubes [59]. Mulliken population analysis [65] indicated the occurrence of charge transfer ranging from 1.32 e (for Sc) to 0.40 e (for Cu and Pt).
Ab initio molecular dynamics (AIMD) simulations allowed investigating the short-time mobility of Pt and Ru clusters on pristine and boron-doped carbon over a range of different temperatures for a very short period. The evaluation of the mean squared displacement (MSD) provided useful details on the stability of metal atoms at increasing temperature. As expected, metal atoms tend to migrate easier on an undoped surface compared to the boron-doped one. However, the differences in MSD over the two surfaces became even less evident as the temperature increases (up to 873 K).
In any case, AIMD operates in a restricted simulation time and size of the system, which are not sufficient to obtain results with general validity. Classical molecular dynamics simulations could guarantee a longer observation time. However, a proper analytical potential should be selected to adequately take into account the metal-support interactions.
Besides tuning the interaction with the adsorbed metal species, the doping can also address the insertion of oxygen functionalities. The adsorption and the stability of epoxide and hydroxyl groups on undoped and B-doped graphene surfaces has been considered [57]. Epoxide functionalities are created under the adsorption of oxygen atoms on the bridge site connecting the boron and carbon atoms. The preferential site for hydroxyl group adsorption was the boron on-top site. Furthermore, the B-doping positively affects the adsorption of oxygen functionalities. A shift of the Fermi level to lower energy and an intensification of the bonding and antibonding states close to the Fermi level occur, contributing therefore to an increase of the binding energies compared to a pure graphene surface. The functionalization with oxygen-containing groups does not have an univocal effect on the metal adsorption properties: the adsorption energies of Pt or Ru atoms rise in the presence of O-groups for pure graphene, whereas a detrimental effect is observed when epoxide and hydroxyl functionalities are incorporated in B-doped graphene. According to the Bader analysis, the metal atoms assume a positive charge, while the functional groups have a net negative charge.

From Single Metal Atoms to Metal Clusters
The single atom systems described above represent just ideal situations, which can be used as a model for the interpretation and the prediction of the more realistic metal clusters or metal nanoparticles. In multi-atom systems, additional interactions can be established between adsorbed atoms that are topologically close (adatom-adatom interactions). A general relationship approximating the binding energy (BE) of metal nanoparticles or clusters with carbon surfaces has been proposed by Efremenko [66], who proposed a linear dependency described by the equation: BE = an + b, where (n) is the number of metal atoms in the cluster, and (a) is an index of surface unsaturation. The intercept (b) accounts for the interactions between the first-layer metal atoms and the carbon surface, and it is strictly related to the electron and structural properties of the metal-support interface. The metal-support interactions between carbon and metal NPs have been disclosed by density functional theory (DFT) supported by experimental investigation. Acharya et al. examined the adsorption of several different metal clusters (Pt 6 , Pt 10 , Pt 32 , Pt 2 Ru 4 , and Pt 4 Ru 6 ) on pristine, boron (B), nitrogen (N), and B,N-co-doped graphene surfaces [57]. The trend in binding energies reproduced the predicted one for the monoatomic counterparts (Pt and Ru atoms): the greatest interaction involved metal clusters and boron-doped graphene, while a weak interaction was produced on pristine graphene. Monitoring the minimum distance between metal clusters and graphene as well as the bond lengths in the metal clusters and the graphene models before and after adsorption provided precious details on atom displacement and surface rearrangement induced by the interaction between metal atoms and carbon atoms at the surface. In general, if the interaction is stronger, the minimum distance will be smaller. Calculated cluster-surface distances reflected the binding energy trend and therefore the interaction order. The higher stability of metal clusters when adsorbed on boron-doped surfaces compared to pristine graphene was confirmed also by AIMD simulations, which predicted lower MSDs for Pt and Ru clusters on B-doped graphene compared to the ones immobilized on pristine graphene. Moreover, as expected, the diffusion of clusters depended on the size of the clusters: the larger the cluster, the lower the diffusion.
Concerning N-doping, analogously to single atoms, metal clusters interacted preferentially with pyridinic nitrogen species [67]. The strong covalent chemical adsorption between modified surfaces and metal clusters affected also the density of states of the metal, carbon, and heteroatoms, and the charge distribution. As an instance, the total density of states (TDOS) of Pd in Pd@PdO clusters is deeply altered after adsorption on N-doped graphene, and simultaneously also, the TDOS of C and N atoms are changed after interacting with Pd@PdO shell clusters [68]. Consequently, the electron density was also affected by the adsorption. Electron density difference (EDD) diagrams for Pd@PdO-N-graphene revealed an electron density enhancement around an N atom associated with an electron density depletion around a C atom, indicating that C-N is weakened due to the interaction with Pd atoms. In addition, according to the Mulliken population analysis, an electron transfer process occurred from Pd cluster to support. However, these variations are not observed when Pd@PdO clusters are immobilized onto pristine graphene. This highlights the central role of heteroatoms in activating the carbon atom of the support and mediating the interaction with the metal clusters. QUAMBO (QUAsi-atomic Minimal Basis set Orbitals) charge analysis can be alternatively used to gain theoretical evidence for charge transfer [69]. Rao et al. applied periodic DFT and QUAMBO charge analysis to explore the interactions at the metal-support interface in the case of Pd nanoclusters varied in size and deposited onto carbon nanofibers with different graphitization degrees and O-functionality content. As expected, DFT calculations evidenced a significantly stronger interaction between Pd clusters and the carbon support, increasing the functionalization degree and decreasing the graphitization degree. QUAMBO analysis revealed the occurrence of charge transfer from the metal to the support, independently of particle size, for the metal atoms close to the interface, but these effects rapidly dissipated with distance from the surface. Combining these results with the information deduced from characterization experiments, the authors concluded that the Pd δ+ contributions are consistent with a partial charge localized and distributed over the first one to two atomic layers at the palladium-carbon interface, while Pd atoms in the second layer are negatively charged through a compensation effect.
Transition metal-based systems are not the unique ones; their strong interactions with the support have been investigated. The confinement of Li 2 S and polysulfides (Li 2 S x ) on the conductive carbon host framework is an attractive challenge to be faced for the development of Li-S batteries with superior performance [70]. The Li-S battery is an emerging rechargeable device with high specific energy, whose commercialization is currently limited by some drawbacks, such as the so-called shuttle effect. The latter is related to the loss of active mass, and it can be minimized by controlling the phase evolution of S at the interface and the interactions between the conductive host framework (typically carbon-based materials) and S-containing guests. From this point of view, the introduction of N-containing functional groups in the carbon hosts is expected to improve the interaction and stabilize the S guests. To predict and rationalize the effect of N-doping on the conductive carbon host and in general in the Li-S cell, host-guest interaction has been simulated by performing DFT calculations devoted to describing the interaction between three surfaces (pristine carbon nanotubes, CNT, N pyr -doped CNT, and N gr -doped CNT) and several S-containing guest models (S atom and Li 2 S x (1< x <4) clusters with high polarity and positive charged terminal Li atoms). Relaxed geometries and binding energy values revealed a comparable trend with respect to transition metals. Also in this case, S-containing species interacted more strongly with N-CNTs than pure CNT. The most stable configuration for the S-species interacting with N pyr consisted of a direct bond between highly positive charged Li atoms and the electron-rich N pyr center. In the case of N gr , S-species are adsorbed on the neighboring carbon atoms, and the nitrogen center acts only as promoter, activating carbon atoms. These promising results demonstrated that N-doping assures the anchoring and the permanent trapping of S-species, and it strengthens the host-guest interactions, thus reducing the loss of active mass. In addition, N atoms can promote electron transfer in the electrocatalytic processes.
To generalize and summarize, the role of heteroatoms in mediating metal-support interaction is the result of the interplay of different electronic effects. It includes several contributions: (a) the match between the electronic structure of the metal and that of the support (b) the number of bonds that are created with the surface atoms (carbon and heteroatoms) (c) breaking of the long-range electronic delocalization (d) local surface deformation, corrugation, and defect introduction (e) electron density redistribution and charge transfer (f) the presence of atoms near the adsorption sites that are available to saturate the eventually dangling bonds.

Predicting the Effect of Heteroatom-Mediated SMSI on the Catalytic Performances
The catalytic consequences of the electronic modification induced by doping can be essentially ascribed to two factors: (i) modifications of the metal adsorption properties because of electronic structure changes; (ii) changes in particle shape; and (iii) increased metal dispersion.

Effects on Molecule Adsorption and Metal Cluster Shape
The adsorption of reactant on the metal surface is the first step of a catalytic cycle, occurring on supported transition metal nanoparticles. The description of the formation of a chemical bond between a surface and a molecule can be performed in terms of a d-band model. According to this model, the bond formation at a transition metal surface is the result of the interaction between the valence states of the adsorbate molecule and the S and D states of the transition metal surface. The interaction of valence states with the S states of the metal produces a shift and broadening of the adsorbate. These renormalized valence states interact with the metal D states, thus generating separate bonding and antibonding states. The strength of the bond is determined by the relative occupancy of these states. The filling of antibonding states results in a bond weakening, which is a function of the energy of the antibonding states relative to the Fermi level. The higher the antibonding states are in energy relative to the Fermi level, the lower their filling degree, and the stronger the bond. Since D antibonding states are always above the D states, the evaluation of the D state energy (the center of the D band) relative to the Fermi level can give a rough idea of the bond strength.
Consequently, the variation in the position of the metal D band induced by heteroatom mediated metal-support interaction will affect the strength of the adsorbate bonding.
Kim et al. [71] compared the adsorption of H 2 on transition metal (TM) supported on boron and pyridine-like nitrogen-doped graphenes. The binding energies were similar for adsorption on TM/N-graphene and TM/B-graphene with values slightly smaller than the one on TM/graphene. The number of maximal H 2 adsorption is lower for TM/N-graphene compared to B-doped and pristine graphene. These trends can be explained considering that upon H 2 adsorption, electrons are transferred from metal D orbitals to the antibonding state of H 2 (backdonation). As described in the previous section, the interaction between transition metals and doped surfaces is stronger than the interaction with undoped surfaces. Consequently, when metals are supported on N or B-containing graphene, their back-donation ability is reduced, and then the H 2 -metal binding is weakened.
The adsorption and decomposition of formic acid have been studied on transition metal (Pt, Pd, Ru) NPs supported on nitrogen-doped graphene. The interaction with a couple of edged pyridinic nitrogen atoms was the most energetically favorable, and promoted the formic acid decomposition to CO 2 and H 2 (except in the case of Ru). The easy decomposition of formic acid is favored by the proximity of a formed adsorbed hydrogen atom and a hydrogen atom of the carboxyl fragment, which are directed toward each other. This is an ideal configuration for the formation of hydrogen and CO 2 molecules [72].
The position of the D-band is a fundamental parameter to describe the adsorption of molecules on metal surfaces, since it can be correlated with the adsorbate binding energy. However, this is not always true. Kim et al. [73] demonstrated that a direct correlation between the D-band center and the adsorption energy of H 2 and CO is not applicable in the case of Pt nanoparticles supported on nitrogen-doped graphene. Otherwise, in these cases, an affordable correlation has been disclosed between the D-band center and the difference in the averaged adsorption energies between H 2 and CO, ∆Ead = Ead(CO) − Ead(H 2 ). Authors also reported the weakened adsorption strength on Pt NPs supported on doped graphene. This difference in adsorption energy was more pronounced for CO adsorption than H 2 , and it is a consequence of the strong metal-support interactions. The latter was also responsible for the bond rearrangement and change in the cluster morphology to maximize the contact with heteroatoms present on the support surfaces ( Figure 5). In any case, the impact of metal-support interactions decreases as the particles become larger [57].
nitrogen-doped graphene. Otherwise, in these cases, an affordable correlation has been disclosed between the D-band center and the difference in the averaged adsorption energies between H2 and CO, ΔEad = Ead(CO) − Ead(H2). Authors also reported the weakened adsorption strength on Pt NPs supported on doped graphene. This difference in adsorption energy was more pronounced for CO adsorption than H2, and it is a consequence of the strong metal-support interactions. The latter was also responsible for the bond rearrangement and change in the cluster morphology to maximize the contact with heteroatoms present on the support surfaces ( Figure 5). In any case, the impact of metalsupport interactions decreases as the particles become larger [57].

SMSI and Metal Dispersion
The interactions between the metal phase and the support at the interface can also determine the dispersion of metal clusters onto the support surface, and finally the particle sizes. In fact, the nucleation and the growth of metal clusters are governed by a delicate balance between metal cohesion and metal support interactions: if the former prevails over the latter, clustering and agglomeration will preferentially occur. Conversely, when the metal binding energy exceeds the cohesive energy of the bulk metal, metal aggregation is prevented, and a high metal dispersion is achieved. This is the case of some transition metals (Sc, Ti, and V) deposited on nitrogen-doped graphene [65]. When pyridine-like functionalities are introduced on the surface of graphene, highly localized acceptor-like states close to the Fermi level are created, leading to a stronger interaction between graphene and metal atoms, and thus favoring a high metal dispersion. The DFT method has been implemented to determine how far the nucleation and growth of Pt NPs onto a carbon matrix (described by a 4 × 4 supercell of graphene) are affected by the presence of nitrogen (pyridinic, pyrrolic, graphitic nitrogen) and and/or sulfur (thiophenic defects) [74]. For this purpose, the interaction energy of Pt on doped graphene surfaces was evaluated in terms of adsorption energy. In all of the cases, a higher interaction energy was obtained for nitrogen or sulfur singularly doped surfaces compared to pristine graphene, and in particular, the following interaction order was observed: N pyrrolic > N graphitic > thiophenic > N pyridinic > graphene. Surprisingly, a weaker interaction and a lower stabilization of Pt nuclei were predicted on a surface co-doped with nitrogen (pyrrolic defects) and sulfur (thiophenic defects). In particular, probing the interaction energy for N and S co-doped surfaces with different dopant location, it appears evident that the interaction and the Pt stabilization decrease by increasing the closeness of the two hetero-defects (nitrogen and sulfur). It seems like sulfur acts as preferential anchoring site for Pt, which is probably due to the electron transfer to the Pt atom.
Metal phase can be deposited on the support according to several strategies. Ning et al. focused on the influence of a catalyst synthesis strategy on the metal-support interaction of Pt on nitrogendoped carbon nanotubes (Pt/NCNTs), comparing ethylene glycol (EG) reduction, sodium borohydride reduction, and impregnation H2 reduction [75].
In the case of the EG reduction method, neutral Pt NPs preferentially interacted with graphitic nitrogen, whereas in sodium borohydride reduction and impregnation H2 reduction, the adsorption of charged metal ion precursors is favored on pyridinic nitrogen.

Probing the Metal-Support Interactions at the Interface
Theoretical simulations can help predict and understand metal-support interactions, and in many cases, they are aimed at supporting experimental evidence deriving from characterization and catalytic tests.
Gaining a deeper insight on metal-support interactions and on the role of surface doping requires a huge amount of information to be collected. The combination of different complementary

SMSI and Metal Dispersion
The interactions between the metal phase and the support at the interface can also determine the dispersion of metal clusters onto the support surface, and finally the particle sizes. In fact, the nucleation and the growth of metal clusters are governed by a delicate balance between metal cohesion and metal support interactions: if the former prevails over the latter, clustering and agglomeration will preferentially occur. Conversely, when the metal binding energy exceeds the cohesive energy of the bulk metal, metal aggregation is prevented, and a high metal dispersion is achieved. This is the case of some transition metals (Sc, Ti, and V) deposited on nitrogen-doped graphene [65]. When pyridine-like functionalities are introduced on the surface of graphene, highly localized acceptor-like states close to the Fermi level are created, leading to a stronger interaction between graphene and metal atoms, and thus favoring a high metal dispersion. The DFT method has been implemented to determine how far the nucleation and growth of Pt NPs onto a carbon matrix (described by a 4 × 4 supercell of graphene) are affected by the presence of nitrogen (pyridinic, pyrrolic, graphitic nitrogen) and and/or sulfur (thiophenic defects) [74]. For this purpose, the interaction energy of Pt on doped graphene surfaces was evaluated in terms of adsorption energy. In all of the cases, a higher interaction energy was obtained for nitrogen or sulfur singularly doped surfaces compared to pristine graphene, and in particular, the following interaction order was observed: N pyrrolic > N graphitic > thiophenic > N pyridinic > graphene. Surprisingly, a weaker interaction and a lower stabilization of Pt nuclei were predicted on a surface co-doped with nitrogen (pyrrolic defects) and sulfur (thiophenic defects). In particular, probing the interaction energy for N and S co-doped surfaces with different dopant location, it appears evident that the interaction and the Pt stabilization decrease by increasing the closeness of the two hetero-defects (nitrogen and sulfur). It seems like sulfur acts as preferential anchoring site for Pt, which is probably due to the electron transfer to the Pt atom.
Metal phase can be deposited on the support according to several strategies. Ning et al. focused on the influence of a catalyst synthesis strategy on the metal-support interaction of Pt on nitrogen-doped carbon nanotubes (Pt/NCNTs), comparing ethylene glycol (EG) reduction, sodium borohydride reduction, and impregnation H 2 reduction [75].
In the case of the EG reduction method, neutral Pt NPs preferentially interacted with graphitic nitrogen, whereas in sodium borohydride reduction and impregnation H 2 reduction, the adsorption of charged metal ion precursors is favored on pyridinic nitrogen.

Probing the Metal-Support Interactions at the Interface
Theoretical simulations can help predict and understand metal-support interactions, and in many cases, they are aimed at supporting experimental evidence deriving from characterization and catalytic tests.
Gaining a deeper insight on metal-support interactions and on the role of surface doping requires a huge amount of information to be collected. The combination of different complementary analytical techniques that are able to probe the metal-support interface then represents the most suitable approach. Considering the nature and the effects of SMSI, the characterization process should provide details on the electronic phenomena and structural rearrangements occurring at the metal-support interface.
Electron microscopies are the ideal tools to give a direct illustration of the morphology and structure of the metal-support interface and describe the metal dispersion on the support.
Among spectroscopic techniques, X-ray photoelectron spectroscopy (XPS) and X-ray absorption spectroscopy (XAS) are the most used ones. The former elucidates the electronic properties of the atoms exposed at the interface. XAS is usually divided in two categories: X-ray absorption near-edge structure (XANES or NEXAFS, including metal K-edge, metal L-edge, and ligand K-edge) and extended X-ray absorption fine structure (EXAFS). Precious evidence on the local structure, the individual environment of atoms, binding geometry, binding distances, and coordination number can be derived from the analysis of XANES and EXAFS spectra.
Infrared spectroscopy as well as thermal techniques (thermal programmed desorption, TPD or thermogravimetric analysis, TGA) can be used to monitor the adsorption of probe molecules onto the metal surfaces, thus attaining indirect information on the electronic structure of metal phases and on eventual changes induced by SMSI. Similarly, electrochemical characterization techniques (e.g., voltammetry, chronoamperometry) allow investigating the electrical behavior of supported metal species, and then deducing the effect resulting from the interactions with the support.
Far from claiming to give an exhaustive overview of the current scenario in the field of material characterization, herein we report selected examples. The latter point out the potential, the applicability, and the limitations of the main characterization techniques that are usually employed for the study of the metal-support interface.

X-ray Photoelectron Spectroscopy (XPS)
XPS is a powerful tool that is crucial for detecting the surface layer element composition, chemical state, and electronic structure of metal and carbon materials. Successively, it has been widely used to study the electronic structure change and the doping effect in supported metal nanostructures on doped carbon materials. Significant shifts in binding energy values are symptomatic of electronic structure alterations resulting from metal-support interactions. However, a core level binding energy shift can be ascribed to three different factors (changes in the local electronic structure, changes in the relaxation process, and cluster charging), and it is often difficult to discriminate among them.
In some cases, the interaction can produce a new electronic state that is detectable as an additional peak (or a peak component in the case of complex signals) in the XPS spectrum.
Among transition metals, noble metals (Pd, Pt, Au) are the most studied systems, due to their large use as catalysts in the form of supported metal nanoparticles.
The strong interaction between Pd species and nitrogen functionalities embedded in carbon-based support has been reported in several papers. In most cases, an increase in the Pd 2+ /Pd 0 ratio and a shift to the higher energy of the Pd 3d signal have been observed and ascribed to a change in the electronic structure or in the coordination environment [76]. The simultaneous analysis of N 1s signals allowed unraveling a prevalent interaction of Pd species with pyridinic and pyrrolic nitrogen [68]. This interaction leads to a charge transfer between the metal and the carbon support, resulting in the stabilization of electron-depleted palladium atoms (Pd δ+ ), whose signal has been detected in XPS spectra [67,69]. In some cases, a very strong bond between Pd 2+ precursor ions and pyridinic nitrogen is created before and during the synthesis, resulting in the low reducibility of palladium species, and thus in the stabilization of isolated Pd 2+ ions in the final catalyst [63].
Arrigo et al. [60,77] proposed a convincing rationalization of the Pd 3d XPS signal of Pd NPs supported on nitrogen-doped carbon nanotubes. Four components were separated under the broad and complex Pd 3d signal ( Figure 6). The signal at lower binding energy (Pd1, 335 eV) was attributed to bulk-like metallic palladium, and was associated with large Pd particles formed under drastic annealing conditions. The Pd2 signal (335.55−335.7 eV) accounted for metallic Pd nanoparticles electronically interacting with the carbon support. These species were predominant in the samples after calcination and reduction. The Pd3 signal (336.9 eV) was related to partially positive charged Pd atoms located at the interface, while the Pd4 component (337.8 eV) is generated by divalent Pd species, namely atomically dispersed Pd. The latter were stabilized by pyridinic nitrogen and possessed low reducibility. After annealing at 300 • C, Pd3 and Pd4 species evolved toward the formation of Pd2 species.
The interpretation of XPS spectra is not trivial, and often, opposing results can appear in the literature. Specifically, in the case of Pd-N interactions, the presence of reduced Pd species [78] or negatively charged Pd has been reported because of an increased electron density by the charge transferred from N to the Pd species. The fact that these species have been revealed when prepared by the pyrolysis of chitosan and melamine or dicyandiamide [79] and Ketjen black carbon [80] suggests a possible role of the synthetic routes in addressing metal-support interactions. electronically interacting with the carbon support. These species were predominant in the samples after calcination and reduction. The Pd3 signal (336.9 eV) was related to partially positive charged Pd atoms located at the interface, while the Pd4 component (337.8 eV) is generated by divalent Pd species, namely atomically dispersed Pd. The latter were stabilized by pyridinic nitrogen and possessed low reducibility. After annealing at 300 °C, Pd3 and Pd4 species evolved toward the formation of Pd2 species. The interpretation of XPS spectra is not trivial, and often, opposing results can appear in the literature. Specifically, in the case of Pd-N interactions, the presence of reduced Pd species [78] or negatively charged Pd has been reported because of an increased electron density by the charge transferred from N to the Pd species. The fact that these species have been revealed when prepared by the pyrolysis of chitosan and melamine or dicyandiamide [79] and Ketjen black carbon [80] suggests a possible role of the synthetic routes in addressing metal-support interactions. Similar results were obtained in the case of Pt-N and Au-N [81] interactions. The Pt 4f spectra usually present three transitions corresponding to the different oxidation states of Pt (Pt 0 , Pt 2+ , and Pt 4+ ). After deposition on nitrogen-doped supports, the oxidation state of the Pt and the ratio between these three components is altered, even if a univocal trend cannot be derived. In most cases, an increase in the Pt 0 /Pt δ+ ratio was exposed [80,82]. In addition, the BE of the Pt is affected by the electronic interaction with nitrogen, undergoing a shift to higher energy values. Although this shift has been reported to be proportional to the number and proximity of nitrogen atoms [83], concomitant size and extra atomic relaxation effects cannot be ruled out. A shift toward higher values was produced also in N 1s BE in species interacting with Pt. In addition, the appearance of new components in the N 1s spectra was assigned alternatively to chemisorbed nitrogen oxide (405.2 eV) [84] deriving from the conversion of graphitic nitrogen, or to a nitride (398.1 eV) resulting from the presence of a metal-nitrogen bond. Similar results were obtained in the case of Pt-N and Au-N [81] interactions. The Pt 4f spectra usually present three transitions corresponding to the different oxidation states of Pt (Pt 0 , Pt 2+ , and Pt 4+ ). After deposition on nitrogen-doped supports, the oxidation state of the Pt and the ratio between these three components is altered, even if a univocal trend cannot be derived. In most cases, an increase in the Pt 0 /Pt δ+ ratio was exposed [80,82]. In addition, the BE of the Pt is affected by the electronic interaction with nitrogen, undergoing a shift to higher energy values. Although this shift has been reported to be proportional to the number and proximity of nitrogen atoms [83], concomitant size and extra atomic relaxation effects cannot be ruled out. A shift toward higher values was produced also in N 1s BE in species interacting with Pt. In addition, the appearance of new components in the N 1s spectra was assigned alternatively to chemisorbed nitrogen oxide (405.2 eV) [84] deriving from the conversion of graphitic nitrogen, or to a nitride (398.1 eV) resulting from the presence of a metal-nitrogen bond.
In any case, the main limitations of XPS are the minimal detectable concentrations and the difficult signal interpretation and deconvolution.

X-ray Absorption Spectroscopy (XAS)
XAS measurements describe intermolecular interactions, probing the valence and local structures of investigated elements. The carbon and nitrogen K-edge and metal L-edge near-edge X-ray absorption fine structure (NEXAFS) spectra together with the extended X-ray absorption fine structure (EXAFS) data provide a detailed picture of chemical and electronic structure speciation at the metal-support interface.
Concerning carbon K-edge spectra, the introduction of nitrogen in graphene sheet has been reported to produce a shift in the C-C σ* peak (286.5 eV), a shift to higher energy and a broadening of the C-C π* peak (292.8 eV), and the presence of an additional sharp peak in the region between the π* and σ* resonance (289.1 eV), which is related to the introduction of structural defects. Moreover, other peaks can appear due to π* transitions deriving from the introduction of functionalities (e.g., amine or imide groups) [89].
Nitrogen K-edge spectra is usually dominated by π*-transitions at 398.5 eV, 400 eV, and 401.4 eV, which can be credited to pyridinic, amino, and pyrrolic/graphitic type N-groups. The immobilization of metal on the surface on N-doped surfaces results in the decrease of the π* resonance intensity in N K-edge NEXAFS spectra (Figure 7). This evidence can be explained by invoking the Stöhr model: a metal-support π back-bonding interaction involves the N1σ-π* transition. Thus, a σ donor is formed by the lone pair orbital of the pyridine-like N, which is stabilized by the back-bonding of metal dπ-pπ orbitals with the N π* antibonding orbitals. Due to the enhanced bond strength, the metal dπ-pπ contributes in a greater extent to the π* orbital, thus producing the decrease of π* resonance intensity. In addition, a shoulder to the lower energy side of the π* resonance appears, indicating a charge localization on the N atom and a charge transfer between donor/acceptor sites [60]. The effect of heteroatom-mediated electron transfer at the metal-support interface has been also reported for other transition metals (Ru [85], Fe [49,86], Co [86,87]) and for ZnCo2O4 perovskite quantum dots [88].
In any case, the main limitations of XPS are the minimal detectable concentrations and the difficult signal interpretation and deconvolution.

X-ray Absorption Spectroscopy (XAS)
XAS measurements describe intermolecular interactions, probing the valence and local structures of investigated elements. The carbon and nitrogen K-edge and metal L-edge near-edge Xray absorption fine structure (NEXAFS) spectra together with the extended X-ray absorption fine structure (EXAFS) data provide a detailed picture of chemical and electronic structure speciation at the metal-support interface.
Concerning carbon K-edge spectra, the introduction of nitrogen in graphene sheet has been reported to produce a shift in the C-C σ* peak (286.5 eV), a shift to higher energy and a broadening of the C-C π* peak (292.8 eV), and the presence of an additional sharp peak in the region between the π* and σ* resonance (289.1 eV), which is related to the introduction of structural defects. Moreover, other peaks can appear due to π* transitions deriving from the introduction of functionalities (e.g., amine or imide groups) [89].
Nitrogen K-edge spectra is usually dominated by π*-transitions at 398.5 eV, 400 eV, and 401.4 eV, which can be credited to pyridinic, amino, and pyrrolic/graphitic type N-groups. The immobilization of metal on the surface on N-doped surfaces results in the decrease of the π* resonance intensity in N K-edge NEXAFS spectra (Figure 7). This evidence can be explained by invoking the Stoḧr model: a metal-support π back-bonding interaction involves the N1σ-π* transition. Thus, a σ donor is formed by the lone pair orbital of the pyridine-like N, which is stabilized by the back-bonding of metal dπ-pπ orbitals with the N π* antibonding orbitals. Due to the enhanced bond strength, the metal dπ-pπ contributes in a greater extent to the π* orbital, thus producing the decrease of π* resonance intensity. In addition, a shoulder to the lower energy side of the π* resonance appears, indicating a charge localization on the N atom and a charge transfer between donor/acceptor sites [60].  In some cases, a new peak at~398.9 eV has been observed after Pd deposition, which is ascribable to the nitrogen species interacting with Pd. [63].
Also, the K-edge and L-edge NEXAFS spectra of the metal are affected by the interaction with the heteroatoms incorporated in the support, inducing changes in electron occupancy in the valence orbital and ligand field environments. As an instance, monitoring Co K-edge XANES provided information on the oxidation state of Co and allowed unraveling strong covalent interactions between Co 3 O 4 and N-graphene, resulting in a partial reduction of Co 3+ [90]. A shift toward higher energy in K-edge of the X-ray absorption near-edge structure spectra (XANES) for Pd@PdO on nitrogen-doped graphene revealed a change of the metal's electronic structure and strong interface interaction between the Pd@PdO particle and nitrogen-doped graphene [80].
In many cases, evidence of XANES is corroborated by EXAFS data, reporting lower metal-metal coordination numbers compared with bulk materials (i.e., high metal dispersion in the form of small nanoparticles), shorter metal-heteroatom and longer metal-metal distances, which are indicative of the strong character of metal-support interaction and demonstrate that these interactions are mediated by the heteroatom present in the support surface [80].

Transmission Electron Microscopies and Other Imaging Techniques
Transmission electron microscopy is a powerful technique that is able to disclose the structural and morphological features of the metal-support interface at the atomic scale. TEM images combined with energy-dispersive X-ray spectroscopy (EDS) and electron energy loss spectroscopy (EELS) allowed appreciating the metal dispersion on support and individuating the presence of defects, unsaturation, heteroatoms, and functional groups. Therefore, TEM represents the ideal tool to unravel the effects of heteroatom-mediated metal-support interactions [91].
In fact, the comparison of TEM images of metal nanoparticles supported on undoped and doped surfaces evidenced the fundamental role of homogeneously distributed heteroatoms in providing nucleation and anchoring sites and promoting a higher metal dispersion. High metal dispersion was detected for different metal species (Au [92], Pd [93], Pt [94,95], Co [87], MoS 2 [96]) on several doped carbon surfaces (N and O-doped carbon nanotubes, N-doped mesoporous carbon, N-doped carbon shells, carbon nanosphere, polyaniline (PANI) fibers). The presence of nitrogen has been demonstrated to assure excellent metal dispersion, even at very high metal loading [93].
Electron micrographs (Figure 8) unveiled the drastic difference between the intimate contact at the interface of Pd particles on defective carbon surfaces, including heteroatom and topological defects, and the undamaged grain boundary of Pd NP supported on an undoped graphitic surface [69]. In some cases, a new peak at ~398.9 eV has been observed after Pd deposition, which is ascribable to the nitrogen species interacting with Pd. [63].
Also, the K-edge and L-edge NEXAFS spectra of the metal are affected by the interaction with the heteroatoms incorporated in the support, inducing changes in electron occupancy in the valence orbital and ligand field environments. As an instance, monitoring Co K-edge XANES provided information on the oxidation state of Co and allowed unraveling strong covalent interactions between Co3O4 and N-graphene, resulting in a partial reduction of Co 3+ [90]. A shift toward higher energy in K-edge of the X-ray absorption near-edge structure spectra (XANES) for Pd@PdO on nitrogen-doped graphene revealed a change of the metal's electronic structure and strong interface interaction between the Pd@PdO particle and nitrogen-doped graphene [80].
In many cases, evidence of XANES is corroborated by EXAFS data, reporting lower metal-metal coordination numbers compared with bulk materials (i.e., high metal dispersion in the form of small nanoparticles), shorter metal-heteroatom and longer metal-metal distances, which are indicative of the strong character of metal-support interaction and demonstrate that these interactions are mediated by the heteroatom present in the support surface [80].

Transmission Electron Microscopies and Other Imaging Techniques
Transmission electron microscopy is a powerful technique that is able to disclose the structural and morphological features of the metal-support interface at the atomic scale. TEM images combined with energy-dispersive X-ray spectroscopy (EDS) and electron energy loss spectroscopy (EELS) allowed appreciating the metal dispersion on support and individuating the presence of defects, unsaturation, heteroatoms, and functional groups. Therefore, TEM represents the ideal tool to unravel the effects of heteroatom-mediated metal-support interactions [91].
In fact, the comparison of TEM images of metal nanoparticles supported on undoped and doped surfaces evidenced the fundamental role of homogeneously distributed heteroatoms in providing nucleation and anchoring sites and promoting a higher metal dispersion. High metal dispersion was detected for different metal species (Au [92], Pd [93], Pt [94,95], Co [87], MoS2 [96]) on several doped carbon surfaces (N and O-doped carbon nanotubes, N-doped mesoporous carbon, N-doped carbon shells, carbon nanosphere, polyaniline (PANI) fibers). The presence of nitrogen has been demonstrated to assure excellent metal dispersion, even at very high metal loading [93].
Electron micrographs (Figure 8) unveiled the drastic difference between the intimate contact at the interface of Pd particles on defective carbon surfaces, including heteroatom and topological defects, and the undamaged grain boundary of Pd NP supported on an undoped graphitic surface [69].  Besides TEM, other imaging techniques have been successfully used to characterize supported metal nanoparticles. For instance, atomic force microscopy (AFM) revealed that Au nucleation occurs selectively at step edges in a freshly cleaned highly oriented pyrolytic graphite (HOPG) surface, whereas a homogeneous dispersion of Au nanoparticles in the middle of terraces has been attained on HNO3-treated surfaces due to the introduction of O functionalities and defects acting as preferential anchoring sites [81].
Scanning transmission X-ray microscopy (STXM) uses a nanoscaled focused soft X-ray beam to provide microscopic measurements with high spatial resolution. It has been successfully applied to obtain the spatial mapping of reduced and stoichiometric Co3O4 on nitrogen-doped graphene, demonstrating the nitrogen and oxygen sites are favored sites for the immobilization of Co3O4 nanocrystals [90].

Other Techniques: Electrochemical and Adsorption Properties Characterization
The electronic promotional effects originating from metal-support interactions can result in an enhanced electrocatalytic activity of the supported metal species. A useful parameter that can be used as indicator of the electrocatalytic activity is the electrochemical surface area (ECSA). ECSA can be calculated from the charges of H2 adsorption/desorption peaks in cyclic voltammetry (CV) curves, according to the equation: where QH is the amount of charge for the electrodesorption of hydrogen atoms on the metal surface, M is the mass of metal, and Qref is the metal surface density [97].
Pt-based nanoparticles dispersed on several N-doped carbon surfaces have been demonstrated to possess a higher ECSA compared to the Pt nanoparticles supported on undoped carbon surfaces. This effect could be associated with the enhanced adsorption of protons and electrolyte ions in the electrical double layer, due to the reinforcement in π bonding and higher wettability at the electrode/electrolyte interface associated with nitrogen species. Besides TEM, other imaging techniques have been successfully used to characterize supported metal nanoparticles. For instance, atomic force microscopy (AFM) revealed that Au nucleation occurs selectively at step edges in a freshly cleaned highly oriented pyrolytic graphite (HOPG) surface, whereas a homogeneous dispersion of Au nanoparticles in the middle of terraces has been attained on HNO 3 -treated surfaces due to the introduction of O functionalities and defects acting as preferential anchoring sites [81].
Scanning transmission X-ray microscopy (STXM) uses a nanoscaled focused soft X-ray beam to provide microscopic measurements with high spatial resolution. It has been successfully applied to obtain the spatial mapping of reduced and stoichiometric Co 3 O 4 on nitrogen-doped graphene, demonstrating the nitrogen and oxygen sites are favored sites for the immobilization of Co 3 O 4 nanocrystals [90].

Other Techniques: Electrochemical and Adsorption Properties Characterization
The electronic promotional effects originating from metal-support interactions can result in an enhanced electrocatalytic activity of the supported metal species. A useful parameter that can be used as indicator of the electrocatalytic activity is the electrochemical surface area (ECSA). ECSA can be calculated from the charges of H 2 adsorption/desorption peaks in cyclic voltammetry (CV) curves, according to the equation: where Q H is the amount of charge for the electrodesorption of hydrogen atoms on the metal surface, M is the mass of metal, and Q ref is the metal surface density [97].
Pt-based nanoparticles dispersed on several N-doped carbon surfaces have been demonstrated to possess a higher ECSA compared to the Pt nanoparticles supported on undoped carbon surfaces. This effect could be associated with the enhanced adsorption of protons and electrolyte ions in the electrical double layer, due to the reinforcement in π bonding and higher wettability at the electrode/electrolyte interface associated with nitrogen species.
While the electrocatalytic activity is determined by ECSA, the catalytic activity is correlated to the adsorption properties of metal surfaces. CO and H 2 are classical probe molecules that are usually adsorbed on metal surfaces to determine their binding properties. High-pressure infrared (HP-FTIR) spectroscopy, as well as attenuated total reflection (ATR) spectroscopy, have been used to monitor the adsorption of probe gases on the surface of metal supported on doped surfaces.

Effect of Heteroatoms-Metal Interaction in the Liquid Phase Oxidation of Benzyl Alcohol
In the last decades, the liquid phase oxidation of alcohol using noble metals (Pd, Au, Pd, Ru) deposited on carbonaceous supports has been the subject of several studies [98][99][100]. Many reports demonstrated that the presence of heteroatoms on the carbon surface can modify the reactivity of the metal. In this section, we will limit the discussion to some examples of the effect of metal-heteroatom interaction in the liquid phase alcohol oxidation. Our group was among the first to show that Au, Pd, and AuPd immobilized on nitrogen functionalized carbon nanofibers exhibited better performance than the same metals supported on pristine carbon in the benzyl alcohol (entry 1-4) [101,102]. The enhanced performances of N-containing materials were ascribed to the increased metal-support interaction, limiting the leaching and the coarsening of metal NPs [101,103]. In successive studies, we demonstrated that the amount and the nature of the nitrogen groups have influenced both activity and selectivity [52,104]. Pd deposited on pyridine-containing carbon nanofibers CNFs (entry 5) resulted more efficient than pyrrole/pyridine-containing ones (entry 6-7). The presence of pyridine resulted in a higher surface basicity compared to pyrrole/pyridine, facilitating the H-abstraction, which is proposed to be the rate-determining step in alcohol oxidation [100]. Moreover, a strong metal-N interaction prevented the agglomeration and sintering of small Pd particles. The metal-support interaction was also verified by the enhanced metal wetting observed. TEM studies showed that small particles better wet the carbon surface, resulting in oval-like particles rather than the round particle shape observed in colloids [52]. Ravat et al. demonstrated that B-doped carbon structures are better support that pristine carbon for Pd nanoparticles [105]. The presence of B induced a strong Pd-B interaction, increasing the activity and the durability of the catalyst during the benzyl alcohol oxidation. The effect of carbon dopant on the activity of non-noble metal was also reported. Sun et al. demonstrated the highest activity of cobalt nanoparticles deposited on nitrogen-doped graphitic carbon compared to Co on bare carbon in the veratryl alcohol oxidation [106]. The presence of N showed a positive synergistic effect on Co, influencing the electronic properties of the metal.
Su et al. were among the first to clearly demonstrate that the role of nitrogen on the enhancing activity of cobalt in benzyl alcohol oxidation is due to an electron transfer at the metal/nitrogen-doped carbon interface [107]. Co/N-Carbon catalysts with different nitrogen amounts, in the form of pyridinic and graphitic N (Figure 9b) (Co@NC-0.3, Co@NC-2, Co@NC-4) were prepared using different amounts of C 3 N 4 as the nitrogen source. The catalytic results clearly showed that the activity increases by increasing the nitrogen amount (entry 10-12, Table 1). By a combination of X-ray photoelectron spectroscopy (XPS) and extended X-Ray absorption fine structure (EXAFS) studies, they have proved that the presence of nitrogen can modify the electron density of Co particles, boosting the oxidative properties of the metal (Figure 9f). XPS showed that by increasing the nitrogen amount, the electron density on Co decreases ( Figure 9b). As cobalt is in direct contact with N (Figure 9a), there is a flow of electrons (Mott-Schottky effect) at the metal/nitrogen interface from Co to N (Figure 9d) that depends on the nitrogen amount. Nitrogen-containing carbons with a higher flat band potential will accept electrons from cobalt nanoparticles till the Fermi level reaches an equilibrium.

Conclusions
The interaction between metal nanoparticles and heteroatoms (N, B, P, S) present on the carbon surface has been a matter of debate in recent years. The functional groups on the carbons have been reported to be (1) effective for anchoring metal nanoparticles, and (2) alter the electronic properties. Consequently, the reactivity of the metal nanoparticles is often improved. A combination of theoretical studies, characterization by spectroscopy (i.e., XPS and EXAFS), microscopy (i.e., TEM), and catalytic experiments has been utilized to understand the nature of metal-heteroatom interaction. Although theoretical models and predictions provide useful tools to rationalize experimental evidence and trends, a huge gap still separates the most common models from the actual systems studied under realistic conditions. Computational approaches suffer from several limitations in terms of time, power, cost, and knowledge. Concerning the rational study of metal-support interactions, great efforts should be made to develop new models. These models should consider some important features, such as the role of long-range dispersion forces at the interface or the reciprocal influence of defects, heteroatoms, and impurity or molecules adsorbed in the proximity of the interface.