Next Article in Journal
The Antimicrobial Peptide Cathelicidin Exerts Immunomodulatory Effects via Scavenger Receptors
Previous Article in Journal
Integrative Proteome Analysis Revels 3-Hydroxybutyrate Exerts Neuroprotective Effect by Influencing Chromatin Bivalency
Previous Article in Special Issue
Theoretical Analysis of Polynuclear Zinc Complexes Isolobally Related to Hydrocarbons
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploring the Potential Energy Surface of Pt6 Sub-Nano Clusters Deposited over Graphene

1
Departamento de Química Física y Analítica, Universidad de Oviedo, 33006 Oviedo, Spain
2
Departamento de Bioquímica y Biología Molecular y Celular, Facultad de Ciencias, Universidad de Zaragoza, 50009 Zaragoza, Spain
3
Instituto de Biocomputación y Física de Sistemas Complejos (BIFI), Universidad de Zaragoza, 50009 Zaragoza, Spain
4
Departamento de Química Física, Universidad de Zaragoza, 50009 Zaragoza, Spain
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(1), 870; https://doi.org/10.3390/ijms24010870
Submission received: 29 November 2022 / Revised: 22 December 2022 / Accepted: 30 December 2022 / Published: 3 January 2023

Abstract

:
Catalytic systems based on sub-nanoclusters deposited over different supports are promising for very relevant chemical transformations such as many electrocatalytic processes as the ORR. These systems have been demonstrated to be very fluxional, as they are able to change shape and interconvert between each other either alone or in the presence of adsorbates. In addition, an accurate representation of their catalytic activity requires the consideration of ensemble effects and not a single structure alone. In this sense, a reliable theoretical methodology should assure an accurate and extensive exploration of the potential energy surface to include all the relevant structures and with correct relative energies. In this context, we applied DFT in conjunction with global optimization techniques to obtain and analyze the characteristics of the many local minima of Pt6 sub-nanoclusters over a carbon-based support (graphene)—a system with electrocatalytic relevance. We also analyzed the magnetism and the charge transfer between the clusters and the support and paid special attention to the dependence of dispersion effects on the ensemble characteristics. We found that the ensembles computed with and without dispersion corrections are qualitatively similar, especially for the lowest-in-energy clusters, which we attribute to a (mainly) covalent binding to the surface. However, there are some significant variations in the relative stability of some clusters, which would significantly affect their population in the ensemble composition.

1. Introduction

In the last years, supported noble metal sub-nanoclusters (SNC) have proven to be promising heterogeneous catalysts due to their high surface area, which allows to reduce the amount of catalyst loading, and their ability to activate inert chemical bonds [1,2,3]. Hence, a wide range of catalytic processes involving SNC have been reported in the literature [4,5]. For example, Pt SNC are being investigated as catalysts for technologically important processes, such as the oxygen reduction reaction (ORR) [6,7,8,9]. This way, a detailed theoretical understanding of the structure and properties of these Pt-based SNC is required.
Realistic investigations on the electronic structure of Pt SNC should consider its fluxional behavior [10,11] as ensemble effects are of the utmost importance in SNC-catalyzed processes, and the catalytically most active species is likely to be a structure different from the global minimum (GM) [12]. In addition, the support plays an important role in the system activity, as metal–metal and metal–support interactions lead to complex potential energy surfaces (PES). In this regard, we note that there are some reports on the interaction between Pt-clusters and carbon-based supports that reveal that the cluster size plays an important role on the nature of the interaction [13,14,15,16,17]. For example, based on energy decomposition analysis and changes in the electron density upon adsorption of Pt SNC on a graphene model, it was found that smaller clusters are prone to bond covalently, while in larger structures, van der Waals interactions are predominant [16]. Moreover, some authors have reported that there is a charge transfer from the metallic cluster to the support [18,19,20]. Besides, the Pt–C interaction is significantly affected by the support curvature of the carbon support, which was used by Yang et al. to strengthen adsorption of methanol onto Pt7 SNC [15].
Therefore, a thorough and accurate sampling of the PES is required to faithfully model a realistic ensemble. In this context, we performed a detailed analysis of Pt6 SNC supported on graphene by using density functional theory (DFT) combined with global optimization (GO) techniques. Moreover, special attention was paid to analyzing dispersion effects on the cluster ensemble composition, which was made by comparison between results performed with and without dispersion corrections. We chose a catalytically relevant system, as Pt6 clusters deposited over carbon-based supports have proven to be active in such important processes as ORR [20] and hydrogen electro-oxidation [21]. This system has been previously studied by other authors: Nakajima and co-workers studied Pt6 clusters deposited over graphene by pre-optimizing the clusters in the gas phase and then depositing them in the surface [20]; a similar strategy (combined with molecular dynamics simulations) was applied by Da Silva et al., but they found a different geometry for the GM [22]. In addition, to the best of our knowledge, complete GO in which clusters are directly generated over the surface, which is required for modelling the PES of high-fluxional SNC [11], has not been performed yet. Note that, as previously introduced, SNC exhibit a remarkable ability to change shapes and re-adapt their structure in the presence of adsorbates and/or a support [23], and thus, for an accurate exploration of the PES, it is mandatory to generate (hundreds of) the initial structures directly over the support to avoid (or at least mitigate) incomplete sampling effects derived from depositing previously optimized SNC on it.
We note that while a detailed knowledge on the low-lying energy structures of interfaces decorated with Pt6 SNC is required for an accurate representation and understanding of the system, modelling electrocatalytic systems is a very difficult task. As a consequence, a variety of novel computational methodologies has been proposed in the last years [24,25,26,27,28,29,30], and other different effects are also expected to influence the ensemble composition, such as the applied potential [31,32] or the experimental conditions (such as pH) [33,34], whose investigation is out of the scope of this contribution.

2. Results and Discussion

We first performed a GO for Pt6 sub-nano clusters deposited over a graphene layer by means of PBE exchange-correlation functional, including a D3BJ scheme for accounting for dispersion interactions [35,36]. During geometry optimizations, we allowed a free relaxation of the magnetic structure of the system, which was further refined by means of single-point calculations by using the tetrahedron method with Blöchl corrections (see Computational Details section). The main motivation for such procedure was obtaining more accurate relative energies, which will affect the ensemble composition as well as magnetic states. In this regard, some reports have put forward the relevance of magnetic interactions within the catalysts (and between the catalysts and reactants) in their activity [37,38,39,40,41,42,43].
There were 123 structures within a relative energy cutoff of 1.0 eV, out of which the 15 most stable systems are provided in Figure 1 (the coordinates of the 50 most stable structures are provided in the Supplementary Materials). The various structures are named as Pt6-#, where # corresponds to the order of the structure in increasing energy relative to the GM (Pt6-I).
As shown in Figure 1, the GM corresponds to a double-square-shaped structure that interacts with the graphene support by means of a bridge coordination to C–C bonds (i.e., Pt atoms lie in the center of the bond). Such kind of Pt–C interaction has been reported to be the most favorable one for the interaction between Pt atoms and graphene-like surfaces [18,22,32,44], which agrees with our observations. On the contrary, to the best of our knowledge, the GM shape had not been previously proposed by other authors but has been reported as a low-lying structure for the gas-phase system, with an energy difference from the GM that ranges between 0.034 eV [22] and ~0.3 eV [45]. Nakajima et al. obtained a GM in which Pt interacts with the support by means of two Pt atoms, forming a square Pt4-core that was completed by two bridge Pt atoms that formed a triangular-like geometry with two different Pt–Pt bonds although they then matched the catalytically relevant active species to a different structure by means of spectroscopy measurements (see ref. [20]). The GM obtained by Da Silva and co-workers is related to the former one [22], but the structure exhibits a more planar character, being closer to a triangle-like geometry that has been reported as the GM for gas-phase Pt6 clusters by several authors by means of pure functionals [20,45,46,47,48]. However, we did not find any structure directly derived from the aforementioned triangular gas-phase cluster in the set of low energy structures, with the partial exception of Pt6-V (and related structures such as Pt6-VI, Pt6-VIII, and Pt6-IX), which resembles a distorted Pt6-triangle, in which the central Pt of a hypothetical three-atom base is displaced upwards, and the apical Pt atom bulges. In order to rule out that the difference is due to an incomplete PES sampling, we searched for similar structures to that reported by the aforementioned authors in the whole set of local minima we obtained (a total of 244 structures). We found relative energies of 0.88 eV and 1.60 eV for two structures that are very related to the GM as proposed by da Silva [22], as they correspond to Pt6-triangles that interact with the surface by two of the three Pt atoms that form the triangle base. Note that the energy difference between them is due to a different pattern of interaction with the support and structural distortions in the Pt6 core. The geometries of these structures, which are identified as Pt6-min da Silva and Pt6-min’ da Silva, are provided in the Supplementary Materials (Table S2). We also found a relative energy of 1.04 eV for the following structure in energy ordering (or a closely related structure) reported by the same authors (identified as Pt6-min2 da Silva in Table S2), which corresponds to a Pt6-triangle that interacts with the support by means of three Pt atoms. Moreover, we obtained a structure that corresponds to a flat Pt6-traingle that interacts with the support by means of van der Waals interactions at a distance of the surface of about 3.1 Å. Such structure presents a relative energy of 1.15 eV, with its geometry provided in Table S2 (Pt6-min6 da Silva). When comparing our local minima with those reported by Nakajima and co-workers, we found that their GM is closely related to a structure that, in our set, has a relative energy of 0.87 eV (Table S2, Pt6-min Nakajima) [20]. They also proposed a low-energy structure with triangular, prismatic shape for which we obtained a relative energy of 0.92 eV (Table S2, Pt6-str7 Nakajima) and the aforementioned van der Waals structure although with a relative energy significantly higher than da Silva: 0.433 eV vs. 0.0668 eV, respectively, while in our case, it was 1.15 eV, and the structure is not exactly the same, as it is partially displaced with respect to the support. Overall, these differences put forward the extreme sensitivity of these systems to the computational approach, including both DFT calculations and PES sampling procedure.
If we turn back to our ensemble (Figure 1), we can see that the four most stable structures, Pt6-I to Pt6-IV, with a maximum energy difference between them of 0.14 eV, exhibit a closely related shape. Both Pt6-I and Pt6-II show a planar double-square shape, which only differs on the interaction mode with the graphene sheet (see top view of Figure 1). Namely, while Pt6-I binds to C–C bonds of the support in a zig-zag manner, Pt6-II binds to parallel C–C bonds, with the energy difference between both structures being only 0.02 eV. The geometrical structure of Pt6-III and Pt6-IV consists of two Pt4 squares that are joined by a Pt–Pt bond, having a shape that resembles a hinge, and as for the previous couple of structures, the main difference between both minima (which translates into an energy difference of only 0.03 eV) is due to the different orientation with respect to the support (top view of Figure 1). The following local minima by relative energy ordering are Pt6-V, (ΔE = 0.17 eV). This system exhibits a significantly different shape, whose geometry (which has already been introduced) might be described as a core of 5 Pt atoms forming two triangles that are joined together by means of a central atom—which has been reported by some of us as a relevant local minimum for Pt5/graphite [32]—and an additional Pt atom on the top that forms an angle that breaks the planarity of the system. As for the previous SNC, Pt6-V interacts with C–C bonds in a bridge fashion, while in this case, it is bonded to the support by two Pt atoms (instead of three). Note that this structure is intimately related to Pt6-VI, Pt6-VIII, and Pt6-IX, which present essentially the same geometrical structure but show a different interaction pattern with the support (see Figure 1, top view). This change in the interaction mode translates into significant energy differences in the relative energy, which is, for example, 0.2 eV higher for Pt6-VIII than for Pt6-V.
We do not comment in detail the geometrical structures of all the other local minima shown in Figure 1, but we can see that all of them interact with the C–C bonds of the support in the bridge positions and by two or three Pt atoms, which agrees with reports from other authors [22], and that there are some other structures whose main difference is the interaction pattern with the support (i.e., Pt6-VII, Pt6-X, and Pt6-XI). To our understanding, this result puts forward the relevant role of the support in SNC stability and thus the importance of an exhaustive sampling of the surface.
Furthermore, the most stable structures are associated to highly planar geometries, which is in line with previous observations that more flat structures are favored by pure functionals (such as PBE), while hybrid ones favor more globular geometries [45,49].
With respect to magnetic states, most of the clusters have magnetic moments close to 2.0 μB. Namely, this is the case for the nine most stable structures (Pt6-I to Pt6-IX), in line with previous reports for both gas-phase and graphene-supported Pt6 clusters [22,45]. Moreover, most of the structures reported herein have different magnetic states, which are very close in energy, which is also consistent with previous spin-state analysis of Pt6 systems [45,50].
In this regard, turning back to the GM, we found that its ground state corresponds to a triplet state with two unpaired electrons overall and a total magnetic moment of 1.9 μB. However, there is a non-magnetic state that is only 0.09 eV higher in energy, which we refer to as Pt6-I’ (see Figure 2), and would involve that both structures may be accessible even at low temperatures (if we do not consider limitations of forbidden spin crossover effects). Note that while Pt6-I corresponds to a magnetic structure with ferromagnetic coupling between all the individual magnetic moments of the various Pt atoms, Pt6-I’ exhibits an antiferromagnetic state, in which the Pt atoms directly bonded to the support show negative individual magnetic moments, while those that are only bonded to other Pt atoms have positive magnetic moments. While the absolute values of the magnetism of atoms bonded to graphene have similar magnitudes (about 0.1 μB), those of Pt atoms that are not bonded to C differ from 0.49–0.58 μB in Pt6-I to 0.09–0.13 μB in Pt6-I’. The magnitude of the individual magnetic moments is in the same range as that reported by Kumar and Kawazoe for gas-phase Ptn clusters [48]. We note that a further discussion of magnetic states is out of the scope of this contribution, but overall, these results support the importance of an adequate consideration of the magnetic state of the metal atom, as it is likely to affect the catalysis [38,51].
We then analyzed the direction of charge transference when the cluster binds to the support. For that, we calculated Bader charges, with the results for the GM shown in Figure 3 (see Figure S1 for the charges of selected local minima with significantly different geometrical structure). We can see that the Pt atoms that are directly bonded to the graphene support have positive charges: 0.31 au for the central Pt atom and 0.10 au for Pt atoms at the ends of the cluster, and this result holds for all the analyzed local minima (Figure S1). With respect to the other three Pt atoms, the top central one exhibits a tiny positive charge (0.02 au), while those at the corners are negatively charged (−0.17 au). Our results indicate that the charge transfer takes place from the cluster—which would have an effective charge of 0.19 au (obtained by summing up the individual charges of the six Pt atoms)—to the support. This result agrees with reports from other authors [18,19,20,22].
Finally, we studied the effect of dispersion corrections in the relative energy and geometry of the clusters. For that, we performed the GO calculations without including dispersion corrections. To differentiate these structures from the dispersion correction-optimized ones, we refer to them as Pt6(no-D)-#, where # indicates the structure position in the relative energy ordering. The 10 most stable structures are shown in Figure 4, in which we also include the equivalent dispersion-corrected structure (in parentheses), and in Table 1, in which we also include the energy values for the analogous dispersion-corrected structure. At first glance, we can see that the geometries of the structures obtained with and without dispersion corrections are very similar. However, there are some significant changes in relative stability that lead to variations in the relative order. For example, there is a switch in order between Pt6(no-D)-IV and Pt6(no-D)-V, which were Pt6-V and Pt6-IV in the dispersion-corrected ensemble (see Figure 4), although such change only involves a minor relative energy difference of up to 0.03 eV (see Table 1). There is also an order transposition in Pt6(no-D)-IX (former Pt6-VII), which switches from position 7 to 9. Contrary to the previous case, this structure is significantly more stable (0.18 eV) in the dispersion-corrected scheme although both geometries are very similar. Another relevant position change corresponds to Pt6(no-D)-X, which was Pt6-XV in the dispersion-corrected scheme, and in this case, the relative energy within the ensemble is 0.07 eV lower in the non-dispersion corrected set of structures.
In order to rationalize this observation, we recurred to the dispersion energy term obtained from D3BJ scheme (Edisp in Table 1). We see that most values are relatively similar, about −7.92 eV. However, for Pt6-VI,I it is −8.115 eV, which explains the significant destabilization when excluding such correction. The opposite trend is also observed, as, for example, structures Pt6(no-D)-VII, VIII, and X (Pt6-VIII, IX, and XV in the dispersion-corrected set, respectively) are comparatively stabilized (by 0.05–0.07 eV, see Table 1) in the non-dispersion ensemble, which correlates with the lower weight of dispersion correction (−7.87 eV on average).
Although overall our findings agree with those reported in ref. [22], in which the authors observed that dispersion correction effects barely affect the cluster structure, we found that, for some cases, dispersion corrections are important in providing correct relative energies (and thus cluster populations), as some structures are significantly more affected than others, and this effect might be relevant in the catalytic activity.
As previously introduced, the interaction of small clusters with carbon-based surfaces has previously been attributed to predominant covalent interactions with the support [16]. In qualitative terms, this result correlates with our findings, as the most stable structures obtained by means of dispersion and non-dispersion-corrected procedures show very similar geometries, and we did not find any low-lying structure bonded to the support in the characteristic parallel manner that would result from predominant van der Waals interactions.

3. Computational Details

Spin-polarized density functional theory (DFT) calculations were performed by means of the Vienna Ab initio Simulation Package (VASP), Version 5.4.4 [52,53,54]. We applied the PBE exchange-correlation functional [55] in conjunction with the projector augmented wave (PAW) method [56,57] to represent interactions between core and valence electrons. Clusters were directly generated and optimized over a p (6 × 6) graphene surface (a and b lattice parameters of 14.777 Å) with a vacuum of 18 Å (between graphene layers) to avoid interactions between parallel images. For geometry optimizations, we considered a Gaussian smearing (width 0.1 eV) and a cutoff of 400 eV for plane waves. For the electronic minimization algorithm, we selected the “ALGO = Fast” option, which selects a mixture of the Davidson and RMM-DIIS algorithms, while for geometry relaxation, we selected “IBRION = 2”, which applies a conjugate-gradient algorithm. The convergence criteria for energy calculations was set to 10−6 eV (“EDIFF = 1e-06”), as the criteria for the geometry optimization were a difference lower than 10−5 eV between two consecutive steps (the default VASP value, which is EDIFF × 10). We also considered a real-space evaluation of projector operators (“LREAL = Auto”). We further performed single-point energy calculations to refine the previous results by using the tetrahedron method with Blöchl corrections, a cutoff of 500 eV and a blocked Davidson algorithm for optimizing the orbitals (“ALGO = Normal”). The Brillouin zone was integrated by a 1 × 1 × 1 K-point mesh for geometry optimizations, which was increased to 5 × 5 × 1 in single-point calculations. Dispersion interactions were accounted for by means of D3BJ scheme developed by Grimme and co-workers [35,36]. In a first step, we considered 250 and 200 initial structures for GO optimization without and with dispersion corrections, respectively. In order to minimize PES incomplete sampling, we also re-optimized the 20 most stable structures of each set with the settings of the other. Then, we added 100 additional structures to make sure that the PES was properly sampled. Structures were generated and filtered thorough PGOPT program suite developed in Alexandrova’s group, which uses a bond length distribution algorithm (BLDA) [58,59].

4. Conclusions

The potential energy surface for Pt6 sub-nano clusters deposited over a graphene layer was explored using DFT methodology and global optimization techniques. The global minimum energies correspond to structures featuring a planar double-square shape for Pt6 atoms interacting with the graphene support by bridge coordination of three Pt atoms to C–C bonds (Pt6-I to IV). We also found other low-lying structures with relatively high planarity as well as some others with more prismatic character. In addition, although Pt atoms always interact with the C–C bond by means of a bridge coordination, we observed some different patterns of interaction with the support due to several interaction modes with the graphene support (involving different combinations of C–C bonds). Therefore, an extensive sampling in global optimization techniques is required to accurately sample the PES of graphene-supported Pt6 SNC. Analysis of magnetic properties shows that the global minimum (GM) structure features a triplet state with two unpaired electrons and a total magnetic moment of 1.9 μB. Inspection of atomic charges of the GM reveals a transfer of electron density from the cluster to the support of 0.19 au. Finally, the effect of dispersion interactions introduced by D3BJ corrections was analyzed by comparison to global optimization performed without dispersion interactions. Although the geometries are not very affected (which we associate to predominant covalent bonding with the support), and the nature of the GM does not change, and the dispersion corrections affect structures in a different manner. While the relative energy does not change much in general, it does for some structures, which would distort ensemble population and puts forward the importance of including dispersion corrections in the calculations.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms24010870/s1.

Author Contributions

Conceptualization, D.B.-E. and J.M.; methodology, D.B.-E. and S.B.; software, S.B.; validation, formal analysis, data curation, D.B.-E. and J.M.; writing—original draft preparation, V.P. and J.M.; writing—review and editing, D.B.-E., V.P. and J.M.; supervision, V.P. and J.M.; funding acquisition, V.P. and J.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Spanish “Ministerio de Ciencia e Innovación” (MICINN), grant numbers PID2021-122763NB-I00 and PID2021-126212OB-I00, and “Fundación para el Fomento en Asturias de la Investigación Científica Aplicada y Tecnológica (FICyT)”, grant number IDI-2021-000054. S.B. was funded by a predoctoral contract from “Gobierno de Aragón” (Spain).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Resources from the supercomputer “Memento” and the technical expertise and assistance provided by the “Instituto de Biocomputación y Física de Sistemas Complejos” (BIFI) are gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tyo, E.C.; Vajda, S. Catalysis by clusters with precise numbers of atoms. Nat. Nanotechnol. 2015, 10, 577–588. [Google Scholar] [CrossRef] [PubMed]
  2. Chakraborty, I.; Pradeep, T. Atomically Precise Clusters of Noble Metals: Emerging Link between Atoms and Nanoparticles. Chem. Rev. 2017, 117, 8208–8271. [Google Scholar] [CrossRef] [PubMed]
  3. Jin, R.; Zeng, C.; Zhou, M.; Chen, Y. Atomically Precise Colloidal Metal Nanoclusters and Nanoparticles: Fundamentals and Opportunities. Chem. Rev. 2016, 116, 10346–10413. [Google Scholar] [CrossRef]
  4. Liu, L.; Corma, A.A. Metal Catalysts for Heterogeneous Catalysis: From Single Atoms to Nanoclusters and Nanoparticles. Chem. Rev. 2018, 118, 4981–5079. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Zhao, X.; Kong, X.; Wang, F.; Fang, R.; Li, Y. Metal Sub-nanoclusters Confined within Hierarchical Porous Carbons with High Oxidation Activity. Angew. Chem. Int. Ed. 2021, 60, 10842–10849. [Google Scholar] [CrossRef]
  6. Halder, A.; Curtiss, L.A.; Fortunelli, A.; Vajda, S. Perspective: Size selected clusters for catalysis and electrochemistry. J. Chem. Phys. 2018, 148, 110901. [Google Scholar] [CrossRef] [Green Version]
  7. Von Weber, A.; Anderson, S.L. Electrocatalysis by Mass-Selected Ptn Clusters. Acc. Chem. Res. 2016, 49, 2632–2639. [Google Scholar] [CrossRef]
  8. Proch, S.; Wirth, M.; White, H.S.; Anderson, S.L. Strong Effects of Cluster Size and Air Exposure on Oxygen Reduction and Carbon Oxidation Electrocatalysis by Size-Selected Ptn (n ≤ 11) on Glassy Carbon Electrodes. J. Am. Chem. Soc. 2013, 135, 3073–3086. [Google Scholar] [CrossRef]
  9. Ohnuma, A.; Takahashi, K.; Tsunoyama, H.; Inoue, T.; Zhao, P.; Velloth, A.; Ehara, M.; Ichikuni, N.; Tabuchi, M.; Nakajima, A. Enhanced oxygen reduction activity of size-selected platinum subnanocluster catalysts: Ptn (n = 3–9). Catal. Sci. Technol. 2022, 12, 1400–1407. [Google Scholar] [CrossRef]
  10. Zhai, H.; Alexandrova, A.N. Fluxionality of Catalytic Clusters: When It Matters and How to Address It. ACS. Catal. 2017, 7, 1905–1911. [Google Scholar] [CrossRef]
  11. Zhang, Z.; Zandkarimi, B.; Alexandrova, A.N. Ensembles of Metastable States Govern Heterogeneous Catalysis on Dynamic Interfaces. Acc. Chem. Res. 2020, 53, 447–458. [Google Scholar] [CrossRef] [PubMed]
  12. Sun, G.; Sautet, P. Metastable Structures in Cluster Catalysis from First-Principles: Structural Ensemble in Reaction Conditions and Metastability Triggered Reactivity. J. Am. Chem. Soc. 2018, 140, 2812–2820. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Poidevin, C.; Paciok, P.; Heggen, M.; Auer, A.A. High resolution transmission electron microscopy and electronic structure theory investigation of platinum nanoparticles on carbon black. J. Chem. Phys. 2019, 150, 041705. [Google Scholar] [CrossRef] [Green Version]
  14. Campos-Roldaén, C.A.; Ramos-Saénchez, G.; Gonzalez-Huerta, R.G.; Vargas García, J.R.; Balbuena, P.B.; Alonso-Vante, N. Influence of sp3−sp2 Carbon Nanodomains on Metal/Support Interaction, Catalyst Durability, and Catalytic Activity for the Oxygen Reduction Reaction. ACS Appl. Mater. Interfaces 2016, 8, 23260–23269. [Google Scholar] [CrossRef] [PubMed]
  15. Yang, J.; Kim, S.H.; Kwak, S.K.; Song, H.-K. Curvature-Induced Metal−Support Interaction of an Islands-by-Islands Composite of Platinum Catalyst and Carbon Nano-onion for Durable Oxygen Reduction. ACS Appl. Mater. Interfaces 2017, 9, 23302–23308. [Google Scholar] [CrossRef] [PubMed]
  16. Schneider, W.B.; Benedikt, U.; Auer, A.A. Interaction of Platinum Nanoparticles with Graphitic Carbon Structures: A Computational Study. ChemPhysChem 2013, 14, 2984–2989. [Google Scholar] [CrossRef] [PubMed]
  17. Jayabal, S.; Saranya, G.; Geng, D.; Lin, L.-Y.; Meng, X. Insight into the correlation of Pt–support interactions with electro-catalytic activity and durability in fuel cells. J. Mater. Chem. A 2020, 8, 9420–9446. [Google Scholar] [CrossRef]
  18. Ramos-Sanchez, G.; Balbuena, P.B. Interactions of platinum clusters with a graphite substrate. Phys. Chem. Chem. Phys. 2013, 15, 11950–11959. [Google Scholar] [CrossRef]
  19. Verga, L.G.; Aarons, J.; Sarwar, M.; Thompsett, D.; Russell, A.E.; Skylaris, C.-K. Effect of graphene support on large Pt nanoparticles. Phys. Chem. Chem. Phys. 2016, 18, 32713–32722. [Google Scholar] [CrossRef] [Green Version]
  20. Tsunoyama, H.; Ohnuma, A.; Takahashi, K.; Velloth, A.; Ehara, M.; Ichikuni, N.; Tabuchi, M.; Nakajima, A. Enhanced oxygen reduction activity of platinum subnanocluster catalysts through charge redistribution. Chem. Commun. 2019, 55, 12603–12606. [Google Scholar] [CrossRef]
  21. Wang, X.; Zhao, L.; Li, X.; Liu, Y.; Wang, Y.; Yao, Q.; Xie, J.; Xue, Q.; Yan, Z.; Yuan, X.; et al. Atomic-precision Pt6 nanoclusters for enhanced hydrogen electro-oxidation. Nat. Commun. 2022, 13, 1596. [Google Scholar] [CrossRef] [PubMed]
  22. Rêgo, C.R.C.; Tereshchuk, P.; Oliveira, L.N.; Da Silva, J.L.F. Graphene-supported small transition-metal clusters: A density functional theory investigation within van der Waals corrections. Phys. Rev. B 2017, 95, 235422. [Google Scholar] [CrossRef]
  23. Lavroff, R.H.; Morgan, H.W.T.; Zhang, Z.; Poths, P.; Alexandrova, A.N. Ensemble representation of catalytic interfaces: Soloists, orchestras, and everything in-between. Chem. Sci. 2022, 13, 8003–8016. [Google Scholar] [CrossRef]
  24. Nørskov, J.K.; Rossmeisl, J.; Logadottir, A.; Lindqvist, L.; Kitchin, J.R.; Bligaard, T.; Jónsson, H. Origin of the Overpotential for Oxygen Reduction at a Fuel-Cell Cathode. J. Phys. Chem. B 2004, 108, 17886–17892. [Google Scholar] [CrossRef]
  25. Sundararaman, R.; Goddard, W.A., III; Arias, T.A. Grand canonical electronic density-functional theory: Algorithms and applications to electrochemistry. J. Chem. Phys. 2017, 146, 114104. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Hutchison, P.; Rice, P.S.; Warburton, R.E.; Raugei, S.; Hammes-Schiffer, S. Multilevel Computational Studies Reveal the Importance of Axial Ligand for Oxygen Reduction Reaction on Fe−N−C Materials. J. Am. Chem. Soc. 2022, 144, 16524–16534. [Google Scholar] [CrossRef]
  27. Melander, M.M.; Kuisma, M.J.; Christensen, T.E.K.; Honkala, K. Grand-canonical approach to density functional theory of electrocatalytic systems: Thermodynamics of solid-liquid interfaces at constant ion and electrode potentials. J. Chem. Phys. 2019, 150, 041706. [Google Scholar] [CrossRef] [Green Version]
  28. Basdogan, Y.; Maldonado, A.M.; Keith, J.A. Advances and challenges in modeling solvated reaction mechanisms for re-newable fuels and chemicals. WIREs Comput. Mol. Sci. 2019, 10, e1446. [Google Scholar]
  29. Abidi, N.; Lim, K.R.G.; Seh, Z.W.; Steinmann, S.N. Atomistic modeling of electrocatalysis: Are we there yet? WIREs Comput. Mol. Sci. 2021, 11, e1499. [Google Scholar] [CrossRef]
  30. Gauthier, J.A.; Ringe, S.; Dickens, C.F.; Garza, A.J.; Bell, A.T.; Head-Gordon, M.; Nørskov, J.K.; Chan, K. Challenges in Modeling Electrochemical Reaction Energetics with Polarizable Continuum Models. ACS Catal. 2019, 9, 920–931. [Google Scholar] [CrossRef] [Green Version]
  31. Zhang, Z.; Zandkarimi, B.; Munarriz, J.; Dickerson, C.E.; Alexandrova, A.N. Fluxionality of Subnano Clusters Reshapes the Activity Volcano of Electrocatalysis. ChemCatChem 2022, 14, e202200345. [Google Scholar] [CrossRef]
  32. Munarriz, J.; Zhang, Z.; Sautet, P.; Alexandrova, A.N. Graphite-supported Ptn Cluster Electrocatalysts: Major Change of Active Sites as a Function of the Applied Potential. ACS Catal. 2022, 12, 14517–14526. [Google Scholar] [CrossRef]
  33. Duan, Z.; Henkelman, G. Atomic-Scale Mechanisms of Electrochemical Pt Dissolution. ACS Catal. 2021, 11, 14439–14447. [Google Scholar] [CrossRef]
  34. Hu, X.; Chen, S.; Chen, L.; Tian, Y.; Yao, S.; Lu, Z.; Zhang, X.; Zhou, Z. What is the Real Origin of the Activity of Fe−N−C Electrocatalysts in the O2 Reduction Reaction? Critical Roles of Coordinating Pyrrolic N and Axially Adsorbing Species. J. Am. Chem. Soc. 2022, 144, 18144–18152. [Google Scholar] [CrossRef]
  35. Grimme, S. Semiempirical GGA-type Density Functional Constructed with a Long-Range Dispersion Correction. J. Comput. Chem. 2006, 27, 1787–1799. [Google Scholar] [CrossRef] [PubMed]
  36. Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A Consistent and Accurate ab initio Parametrization of Density Functional Dispersion Correction (DFT-D) for the 94 Elements H-Pu. J. Chem. Phys. 2010, 132, 154104–154119. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Gracia, J.; Sharpe, R.; Munarriz, J. Principles determining the activity of magnetic oxides for electron transfer reactions. J. Catal. 2018, 361, 331–338. [Google Scholar] [CrossRef]
  38. Biz, C.; Fianchini, M.; Gracia, J. Strongly Correlated Electrons in Catalysis: Focus on Quantum Exchange. ACS Catal. 2021, 11, 14249–14261. [Google Scholar] [CrossRef]
  39. Biz, C.; Fianchini, M.; Polo, V.; Gracia, J. Magnetism and Heterogeneous Catalysis: In Depth on the Quantum Spin-Exchange Interactions in Pt3M (M = V, Cr, Mn, Fe, Co, Ni, and Y)(111) Alloys. ACS Appl. Mater. Interfaces 2020, 12, 50484–50494. [Google Scholar] [CrossRef]
  40. Li, J.; Ma, J.; Ma, Z.; Zhao, E.; Du, K.; Guo, J.; Ling, T. Spin Effect on Oxygen Electrocatalysis. Adv. Energy Sustain. Res. 2021, 2, 2100034. [Google Scholar] [CrossRef]
  41. Sun, Y.; Ren, X.; Sun, S.; Liu, Z.; Xi, S.; Xu, Z.J. Engineering High-Spin State Cobalt Cations in Spinel Zinc Cobalt Oxide for Spin Channel Propagation and Active Site Enhancement in Water Oxidation. Angew. Chem. Int. Ed. 2021, 60, 14536–14544. [Google Scholar] [CrossRef] [PubMed]
  42. Pan, L.; Ai, M.; Huang, C.; Yin, L.; Liu, X.; Zhang, R.; Wang, S.; Jiang, Z.; Zhang, X.; Zou, J.-J.; et al. Manipulating spin polarization of titanium dioxide for efficient photocatalysis. Nat. Commun. 2020, 11, 418. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Li, X.; Cheng, Z.; Wang, X. Understanding the Mechanism of the Oxygen Evolution Reaction with Consideration of Spin. Electrochem. Energy Rev. 2021, 4, 136–145. [Google Scholar] [CrossRef]
  44. Groves, M.N.; Malardier-Jugroot, C.; Jugroot, M. Improving Platinum Catalyst Durability with a Doped Graphene Support. J. Phys. Chem. C 2012, 116, 10548–10556. [Google Scholar] [CrossRef]
  45. Li, R.; Odunlami, M.; Carbonnière, P. Low-lying Ptn cluster structures (n = 6–10) from global optimizations based on DFT potential energy surfaces: Sensitivity of the chemical ordering with the functional. Comput. Theor. Chem. 2017, 1107, 136–141. [Google Scholar] [CrossRef]
  46. Chaves, A.S.; Rondina, G.G.; Piotrowski, M.J.; Tereshchuk, P.; Da Silva, J.L.F. The Role of Charge States in the Atomic Structure of Cun and Ptn (n = 2−14 atoms) Clusters: A DFT Investigation. J. Phys. Chem. A 2004, 118, 10813–10821. [Google Scholar] [CrossRef]
  47. Xiao, L.; Wang, L. Structures of Platinum Clusters: Planar or Spherical? J. Phys. Chem. A 2004, 108, 8605–8614. [Google Scholar] [CrossRef]
  48. Kumar, V.; Kawazoe, Y. Evolution of atomic and electronic structure of Pt clusters: Planar, layered, pyramidal, cage, cubic, and octahedral growth. Phys. Rev. B 2008, 77, 205418. [Google Scholar] [CrossRef]
  49. Zandkarimi, B.; Alexandrova, A.N. Dynamics of Subnanometer Pt Clusters Can Break the Scaling Relationships in Catalysis. J. Phys. Chem. Lett. 2019, 10, 460–467. [Google Scholar] [CrossRef] [Green Version]
  50. Ignatov, S.K.; Razuvaev, A.G.; Loginova, A.S.; Masunov, A.E. Global Structure Optimization of Pt Clusters Based on the Modified Empirical Potentials, Calibrated using Density Functional Theory. J. Phys. Chem. C 2019, 123, 29024–29036. [Google Scholar] [CrossRef]
  51. Munarriz, J.; Polo, V.; Gracia, J. On the Role of Ferromagnetic Interactions in Highly Active Mo-Based Catalysts for Ammonia Synthesis. ChemPhysChem 2018, 19, 2843–2847. [Google Scholar] [CrossRef]
  52. Kresse, G.; Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. Rev. B Condens. Matter Mater. Phys. 1993, 47, 558–561. [Google Scholar] [CrossRef] [PubMed]
  53. Kresse, G.; Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal–amorphous-semiconductor transition in germanium. Phys. Rev. B Condens. Matter Mater. Phys. 1994, 49, 14251–14269. [Google Scholar] [CrossRef] [PubMed]
  54. Kresse, G.; Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B Condens. Matter Mater. Phys. 1996, 54, 11169–11186. [Google Scholar] [CrossRef] [PubMed]
  55. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Blöchl, P.E. Projector augmented-wave method. Phys. Rev. B Condens. Matter Mater. Phys. 1994, 50, 17953–17979. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Kresse, G.; Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B Condens. Matter Mater. Phys. 1999, 59, 1758–1775. [Google Scholar] [CrossRef]
  58. Zhai, H.; Alexandrova, A.N. Ensemble-average representation of Pt clusters in conditions of catalysis accessed through GPU accelerated deep neural network fitting global optimization. J. Chem. Theory Comput. 2016, 12, 6213–6226. [Google Scholar] [CrossRef] [Green Version]
  59. Zhai, H.; Alexandrova, A.N. Local Fluxionality of Surface-Deposited Cluster Catalysts: The Case of Pt7 on Al2O3. J. Phys. Chem. Lett. 2018, 9, 1696–1702. [Google Scholar] [CrossRef]
Figure 1. Frontal and schematic top views for the fifteen lowest-in-energy local minimum for Pt6/graphene optimized with dispersion-corrected DFT. ΔE corresponds to the energy difference with respect to the GM; the total magnetic moment of the system (per unit cell) is provided in Bohr magnetons. The green dots correspond to the Pt atoms at the top, while the black crosses correspond to the atoms at the bottom.
Figure 1. Frontal and schematic top views for the fifteen lowest-in-energy local minimum for Pt6/graphene optimized with dispersion-corrected DFT. ΔE corresponds to the energy difference with respect to the GM; the total magnetic moment of the system (per unit cell) is provided in Bohr magnetons. The green dots correspond to the Pt atoms at the top, while the black crosses correspond to the atoms at the bottom.
Ijms 24 00870 g001
Figure 2. Spin density, total and individual atomic magnetic moments, and energy different for the two low-lying spin states of Pt6-I. Note that α spin density is shown in light blue and β spin density in light pink (isovalue = 0.03 au).
Figure 2. Spin density, total and individual atomic magnetic moments, and energy different for the two low-lying spin states of Pt6-I. Note that α spin density is shown in light blue and β spin density in light pink (isovalue = 0.03 au).
Ijms 24 00870 g002
Figure 3. Bader charges in (in au) for relevant atoms of the GM (Pt6-I). Positive charges are depicted in.
Figure 3. Bader charges in (in au) for relevant atoms of the GM (Pt6-I). Positive charges are depicted in.
Ijms 24 00870 g003
Figure 4. Frontal view for the ten lowest-in-energy local minimum for Pt6/graphene optimized without dispersion corrections. ΔE corresponds to the energy difference with respect to the GM; the total magnetic moment of the system (per unit cell) is provided in Bohr magnetons.
Figure 4. Frontal view for the ten lowest-in-energy local minimum for Pt6/graphene optimized without dispersion corrections. ΔE corresponds to the energy difference with respect to the GM; the total magnetic moment of the system (per unit cell) is provided in Bohr magnetons.
Ijms 24 00870 g004
Table 1. Correspondence of the 10 most stable structures obtained from the ensemble without dispersion (Pt6(no-D)-#) with those obtained with dispersion corrections (Pt6-#’). All energy values are provided in eV.
Table 1. Correspondence of the 10 most stable structures obtained from the ensemble without dispersion (Pt6(no-D)-#) with those obtained with dispersion corrections (Pt6-#’). All energy values are provided in eV.
Pt6(no-D)-#Pt6-#’∆E(Pt6(no-D)-#)∆E(Pt6-#’)∆(∆E)Edisp
II0.000.000.00−7.920
IIII0.030.020.01−7.930
IIIIII0.120.110.01−7.928
IVV0.140.17−0.03−7.884
VIV0.150.140.01−7.936
VIVI0.210.23−0.02−7.900
VIIVIII0.320.37−0.05−7.877
VIIIIX0.380.43−0.05−7.872
IXVII0.520.340.18−8.115
XXV0.530.60−0.07−7.858
Note that ∆(∆E) corresponds to the difference between the relative energy of equivalent structures obtained without dispersion and with dispersion: ∆(∆E) = ∆E(Pt6(no-D)-#)—∆E(Pt6-#’). Thus, a positive value indicates that the structure is comparatively most stable when computed with dispersion corrections, and a negative value indicates that the structure is comparatively less stable without including dispersion. Edisp is the dispersion energy term obtained from D3BJ scheme.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Barrena-Espés, D.; Boneta, S.; Polo, V.; Munárriz, J. Exploring the Potential Energy Surface of Pt6 Sub-Nano Clusters Deposited over Graphene. Int. J. Mol. Sci. 2023, 24, 870. https://doi.org/10.3390/ijms24010870

AMA Style

Barrena-Espés D, Boneta S, Polo V, Munárriz J. Exploring the Potential Energy Surface of Pt6 Sub-Nano Clusters Deposited over Graphene. International Journal of Molecular Sciences. 2023; 24(1):870. https://doi.org/10.3390/ijms24010870

Chicago/Turabian Style

Barrena-Espés, Daniel, Sergio Boneta, Victor Polo, and Julen Munárriz. 2023. "Exploring the Potential Energy Surface of Pt6 Sub-Nano Clusters Deposited over Graphene" International Journal of Molecular Sciences 24, no. 1: 870. https://doi.org/10.3390/ijms24010870

APA Style

Barrena-Espés, D., Boneta, S., Polo, V., & Munárriz, J. (2023). Exploring the Potential Energy Surface of Pt6 Sub-Nano Clusters Deposited over Graphene. International Journal of Molecular Sciences, 24(1), 870. https://doi.org/10.3390/ijms24010870

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop