Next Article in Journal
Effect of Ionic Liquid on the Determination of Aromatic Amines as Contaminants in Hair Dyes by Liquid Chromatography Coupled to Electrochemical Detection
Next Article in Special Issue
Stability Computations for Isomers of La@Cn (n = 72, 74, 76)
Previous Article in Journal
Novel Antimicrobial Organic Thermal Stabilizer and Co-Stabilizer for Rigid PVC
Previous Article in Special Issue
Coordination Modes and Different Hapticities for Fullerene Organometallic Complexes
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 17
Issue 7
10.3390/molecules17077941
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Platinum Clusters on Vacancy-Type Defects of Nanometer-Sized Graphene Patches

by
Takashi Yumura
1,*,
Tatsuya Awano
2,
Hisayoshi Kobayashi
2 and
Tokio Yamabe
2
1
Department of Chemistry and Materials Technology, Kyoto Institute of Technology, Matsugasaki, Sakyo-ku, Kyoto 606-8585, Japan
2
Nagasaki Institute of Applied Science, 536 Aba-machi, Nagasaki 851-0193, Japan
*
Author to whom correspondence should be addressed.
Molecules 2012, 17(7), 7941-7960; https://doi.org/10.3390/molecules17077941
Submission received: 8 May 2012 / Revised: 11 June 2012 / Accepted: 19 June 2012 / Published: 2 July 2012
(This article belongs to the Special Issue Fullerene Chemistry)
Browse Figures
Versions Notes

Abstract

:
Density functional theory calculations found that spin density distributions of platinum clusters adsorbed on nanometer-size defective graphene patches with zigzag edges deviate strongly from those in the corresponding bare clusters, due to strong Pt-C interactions. In contrast, platinum clusters on the pristine patch have spin density distributions similar to the bare cases. The different spin density distributions come from whether underlying carbon atoms have radical characters or not. In the pristine patch, center carbon atoms do not have spin densities, and they cannot influence radical characters of the absorbed cluster. In contrast, radical characters appear on the defective sites, and thus spin density distributions of the adsorbed clusters are modulated by the Pt-C interactions. Consequently, characters of platinum clusters adsorbed on the sp2 surface can be changed by introducing vacancy-type defects.

1. Introduction

Graphitic carbon materials serve as a support material [1] for anode catalysts such as platinum clusters in proton exchange membrane (PEM) fuel cells [2,3,4,5,6,7,8,9]. The supported Pt clusters catalyze the activation of hydrogen molecules to form protons and electrons on the anode of fuel cells. Consequently, the size of Pt clusters is a crucial parameter in determining their catalytic activity. Actually, clusters of less than 3 nm are more effective for catalyzing the H2 dissociation [9]. During the catalytic reactions, adjacent clusters tend to coalescence, forming larger clusters. Accordingly, their catalytic activity decreases as the reaction proceeds. To retain the catalytic activity of supported Pt catalysts, a plausible approach is to strengthen interactions between Pt clusters and underlying carbon sp2 surface.
To come up with a strategy for constructing carbon supports suitable for Pt catalysts, computational simulations are becoming a powerful tool [9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31]. Several computational studies suggest that disrupting the sp2 surface by introducing defects (vacancy-type and Stone-Wales type) [9,18,21,25,26,29,30,31], dopants (nitrogen or boron impurities) [14,15,16], and mechanical strain [25] can enhance the interactions with Pt clusters. With respect to the formation of the vacancy-type defects, recent high-resolution transmission electron microscopy (TEM) studies [32,33,34,35] show that electron irradiation of graphene creates vacancy-type defects by removing a few carbon atoms from the surface. After these events, unsaturated carbon atoms are generated. Some of the unsaturated atoms make a covalent bond with an adjacent atom to form a five-membered ring, whereas the others remain two-coordinated. These carbon atoms, which cannot be seen in pristine graphene, are more chemically reactive, and thus they serve as sites for strong adsorption of Pt clusters. Reactive carbon atoms can be also found on edges of zigzag-nanoribbons and zigzag-graphene patches, because their frontier orbital coefficients are located on edge carbon atoms [36,37,38,39,40,41,42]. Thus, one can utilize such reactive edge atoms to trap well Pt clusters [22,27]. Previously, we investigated by means of density functional theory (DFT) calculations how a Pt cluster is bound to the nanometer-size rhombic sp2 patch with zigzag-edges (C96H26) (Figure 1) [22]. Such H-terminated sp2 patches are contained in activated carbons as condensed-aromatic-ring fractions [43]. The DFT calculations found that a Pt6 cluster preferentially binds into edge atoms of C96H26 rather than into center atoms. In fact, the Pt6 additions to edge atoms were about 50 kcal/mol more stable than those to center atoms [22].
Molecules 17 07941 g001 550
Figure 1. Optimized structure for C96H26graphene patch.
Figure 1. Optimized structure for C96H26graphene patch.
Molecules 17 07941 g001
Another interesting feature of nanometer-size graphenes is that they have radical character in the ground state, depending on their shape and size [44,45,46,47,48,49,50,51]. Thus, we assume that interactions between a Pt cluster and a radical sp2 patch can modulate the catalytic activity of the supported cluster due to the onset of unpaired electrons on Pt atoms. Based on the assumption, the current study will focus on whether spin states of C96H26 support have a power to influence the properties of the adsorbed Pt clusters. Furthermore, we are interested in how introduction of vacancy-type defects on the radical C96H26 support changes the electronic properties of the surface. These changes would have an impact on determining the properties of Pt clusters adsorbed on the sp2 support. To increase our understanding of the interactions between a radical sp2 support and a Pt cluster, we performed density functional theory (DFT) calculations. The main aim in the current DFT study is to clarify how different electronic properties of C96H26 support with or without vacancy-type defects influence the interactions with Pt clusters, and concomitantly the properties of the absorbed clusters.

2. Results and Discussion

2.1. Platinum Clusters on C96H26 in the Triplet State

To obtain a basic insight on how different spin states of C96H26 patch affect the interactions with Pt clusters, we investigated how Pt6 clusters bind into the sp2 surface. Following the previous study [22], two types of Pt6 cluster were considered, denoted by (i) and (ii) in Figure 2. The DFT calculations found that their triplet states are energetically stable relative to the corresponding singlet states. The stability of spin-polarized states in Pt clusters was also reported by other groups [52,53,54]. We obtained three optimized geometries for Pt6(i) or Pt6(ii) clusters adsorbed on C96H26 (Pt6-C96H26) in the singlet and triplet states. Within the three optimized geometries in Figure 3, one is that the Pt6(i) cluster makes four Pt-C bonds with the sp2 surface, and the other two are distinguished by whether the number of Pt-C bonds formed between the Pt6(ii) cluster and the surface is 2 or 4. In these geometries, optimized lengths of Pt-C bonds range from 2.288 to 2.334 Å.
Molecules 17 07941 g002 550
Figure 2. Optimized structures for bare Pt6 cluster, and their spin density distributions. Isosurface α- and β-spins are given by pink and blue, respectively.
Figure 2. Optimized structures for bare Pt6 cluster, and their spin density distributions. Isosurface α- and β-spins are given by pink and blue, respectively.
Molecules 17 07941 g002
We estimated in Table 1 the energy difference between the triplet and singlet spin states in each configuration, ΔEstate(Pt6-C96H26), defined as [Etotal(triplet state) – Etotal(singlet state)] where Etotal(triplet state) or Etotal(singlet state) is the total energy in each state. As shown in Table 1, the three configurations have negative ΔEstate(Pt6-C96H26) values. These negative ΔEstate(Pt6-C96H26) values indicate that the triplet state of a Pt6-C96H26 configuration is energetically favorable relative to the singlet state, irrespective of the cluster shapes. Furthermore, we see from Table 1 more significant ΔEstate values in the Pt6-C96H26 configurations than those in the bare Pt6 clusters. Thus, the Pt-C interactions influence relative stability of the triplet to singlet states of the Pt6 clusters.
Spin density distributions in the triplet Pt6-C96H26 structures are also displayed in Figure 3, where isosurface α- and β-spins are given by pink and blue, respectively. As shown in Figure 3, the Pt-C interactions induce spin densities on C96H26, although the stable triplet state of pristine C96H26 has radical characters only on edge carbon atoms. Likewise, we see spin densities on the adsorbed Pt cluster in the configurations. Basically their spin density distributions are similar to those in the bare Pt6 clusters (Figure 2), but spin densities slightly decrease on Pt atoms that participates the Pt-C bond formation. The similarity between Pt clusters with and without the carbon support is understandable, because underlying carbon atoms do not have radical characters in pristine C96H26 in the stable triplet state, and thus they cannot perturb the spin density distributions of Pt clusters even though they interact substantially.
Molecules 17 07941 g003 550
Figure 3. Spin density distributions of optimized C96H26 and Pt6-C96H26 configurations in the triplet state. Parts of the optimized geometries, corresponding to the region surrounded by pink hashed lines in Figure 1, are given. Isosurface α- and β-spins are given by pink and blue, respectively. Optimized bond lengths are in Å.
Figure 3. Spin density distributions of optimized C96H26 and Pt6-C96H26 configurations in the triplet state. Parts of the optimized geometries, corresponding to the region surrounded by pink hashed lines in Figure 1, are given. Isosurface α- and β-spins are given by pink and blue, respectively. Optimized bond lengths are in Å.
Molecules 17 07941 g003
Table
Table 1. Energy difference between the singlet and triplet states in Pt6-C96H26 (ΔEstate in kcal/mol) a.
Table 1. Energy difference between the singlet and triplet states in Pt6-C96H26 (ΔEstate in kcal/mol) a.
Pt6(i)Pt6(ii)-(1)Pt6(ii)-(2)
Bare clusters–14.5–5.4–5.4
Clusters on C96H26–31.4–23.6–22.1
a ΔEstate(Pt6-C96H26) = Etotal(triplet state) – Etotal(singlet state). Negative ΔEstate values indicate that the triplet state of a Pt6-C96H26 configuration is energetically stable relative to the singlet state.

2.2. Vacancy-Type Defects Formed by Removing Carbon Atoms from C96H26

Prior to discussing Pt clusters adsorbed on the sp2 surface with vacancy-type defects, we look at how introduction of a vacancy-type defect on C96H26 changes its electronic structures. In this study, we considered the number (n) of carbon atoms removed from C96H26, ranging from 1 to 3. Removed carbon atoms are colored in Figure 4. First, we constructed mono-, di-, and tri-vacancy defects by removing the green atom, the green and blue atoms, and the three colored atoms, respectively. The vacancy-type defects will be denoted by C96–nH26. Using the initial geometries, we obtained optimized structures for the vacancy-type defects in the triple and singlet states. Then, the energy difference between the two spin states was evaluated in each vacancy-type defect, given as ΔEstate(C96–nH26) in Table 2. We can see from Table 2 negative ΔEstate(C96–nH26) values irrespective of the number of carbon atoms removed from C96H26. The negative ΔEstate(C96–nH26) values indicate that each C96–nH26 has energetically stable triplet state.
Molecules 17 07941 g004 550
Figure 4. Vacancy type defects (C96–nH26), constructed by removing a few carbon atoms from C96H26 where n ranges from 1 to 3. Parts of the optimized geometries, corresponding to the region surrounded by pink hashed lines in Figure 1, are given. Optimized bond lengths are given in Å. Their spin density distributions in the triplet state are also given. Isosurface α- and β-spins are given by pink and blue, respectively.
Figure 4. Vacancy type defects (C96–nH26), constructed by removing a few carbon atoms from C96H26 where n ranges from 1 to 3. Parts of the optimized geometries, corresponding to the region surrounded by pink hashed lines in Figure 1, are given. Optimized bond lengths are given in Å. Their spin density distributions in the triplet state are also given. Isosurface α- and β-spins are given by pink and blue, respectively.
Molecules 17 07941 g004
Table
Table 2. Energy difference between the singlet and triplet states in C96-nH26 (ΔEstate in kcal/mol) a.
Table 2. Energy difference between the singlet and triplet states in C96-nH26 (ΔEstate in kcal/mol) a.
N
0123
ΔEstate-12.8–14.0–14.8–31.4
a ΔEstate(C96–nH26) = Etotal(triplet state) – Etotal(singlet state). Negative ΔEstate values indicate that the triplet state of a C96–nH26 configuration is energetically stable relative to the singlet state.
Figure 4 also displays their optimized structures in the triplet state as well as corresponding spin density distributions. As shown in Figure 4, all optimized structures for the vacancy-type defects have some five-membered rings formed by connecting two orange atoms. There is one five-member ring in the mono-vacancy defect (C95H26), while there are two five-member rings in the other vacancy-defects (C94H26 and C93H26). Besides, removing odd-numbered carbon atoms from C96H26 generates one coordinatively unsaturated carbon atom, given by red in Figure 4. In fact, they are bound to only two neighboring atoms. The presence of vacancy-type defects perturbs significantly spin density distributions of pristine C96H26. As displayed in Figure 4, we can see radical carbon atoms around defective sites. In particular, significant spin densities were found in the tri-vacancy defect site. More interestingly, we found that the structural features of the defects have a correlation with how spin densities are distributed. On the unsaturated (red) atoms in the mono- and tri-vacancy-defects, spin densities are distributed on the carbon plane, which come from non-bonding orbitals. In contrast, spin densities on the orange atoms, which are a part of five-membered rings, are found perpendicular to the plane. The onset of radical carbon atoms at the center of the patch differentiates the defective surfaces from the pristine in terms of the interactions with Pt clusters, as will be mentioned below.

2.3. Platinum Clusters on Vacancy-Type Defects (C96H26) Depending on Spin States

2.3.1. Singlet State

Despite the preferences of the triplet state of C96–nH26 over the singlet state, let us first use the singlet state to increase our understanding of how a Ptk cluster interacts with a defective site on the sp2 surface. Following the previous study on interactions between a Pt6 cluster and pristine C96H26, we have a special interest on how the presence of a vacancy-type defect of the patch affects the interactions with a Pt6 cluster. In addition, we will discuss dependences of the interaction energies on size of clusters whose number of contained Pt atoms (k) being smaller than 6. Figure 5–Figure 7 show optimized structures for a Pt6 cluster adsorbed on the mono-, di-, and tri-vacancy-type defects, respectively.
Several modes for the cluster bindings were considered. For example, we obtained six optimized geometries for a Pt6 cluster binding into the mono-vacancy-type defect in Figure 5. The four Pt6-C95H26 structures displayed in Figure 5 are relatively stable in energy. In general, stable Ptk-C95H26 structures have a Ptk-1 moiety contained in stable Ptk-1-C95H26 structures. Of course, there are other possibilities for Ptk binding modes. However, our computational resource is limited, reluctantly we did not obtain other optimized geometries. We evaluated the binding energy in each configuration defined as [Ebind = Etotal(Ptk-C96–nH26) – Etotal(C96–nH26) – Etotal(Ptk)], where k ranges from 1 to 6 (Table 3–Table 8).
Molecules 17 07941 g005 550
Figure 5. Optimized geometries for Pt6 cluster on the mono-vacancy-type defect in the singlet state (Pt6-C95H26). Optimized bond lengths are given in Table 3.
Figure 5. Optimized geometries for Pt6 cluster on the mono-vacancy-type defect in the singlet state (Pt6-C95H26). Optimized bond lengths are given in Table 3.
Molecules 17 07941 g005
Table
Table 3. Key parameters of Ptk on mono-vacancy defect (C95H26) (k is 1 or 6) in Figure 5. Separations of a Pt atom from orange atoms (Pt-C(orange)) and those from the red atom (Pt-C(red)). Separations of carbon atoms from a Pt atom except for the nearest Pt atom (other Pt-C), and those between the two orange atoms (C-C). Bond lengths are in Å. The Ebind and ΔE(Ptk) values are given in kcal/mol. Their definition was given in the text.
Table 3. Key parameters of Ptk on mono-vacancy defect (C95H26) (k is 1 or 6) in Figure 5. Separations of a Pt atom from orange atoms (Pt-C(orange)) and those from the red atom (Pt-C(red)). Separations of carbon atoms from a Pt atom except for the nearest Pt atom (other Pt-C), and those between the two orange atoms (C-C). Bond lengths are in Å. The Ebind and ΔE(Ptk) values are given in kcal/mol. Their definition was given in the text.
EbindPt-C(orange)Pt-C(red)other Pt-CC-C BondΔE(Ptk)
Pt1–145.71.953, 1.9541.942––2.764––
Pt6(i)-(A)–152.01.986, 1.9881.968––2.7454.2
Pt6(i)-(B)–146.61.947, 1.9681.9582.2152.76816.8
Pt6(ii)-(C)–151.61.980, 1.9811.9912.0742.72817.4
Pt6(ii)-(D)–157.11.983, 1.9832.0102.1142.7298.1
Pt6(ii)-(E)–136.71.978, 1.9761.9912.200, 2.083, 2.2522.80315.3
Pt6(ii)-(F)–143.41.986, 1.9771.9822.1132.7406.3
Molecules 17 07941 g006 550
Figure 6. Optimized geometries for Pt6 cluster on the di-vacancy-type defect in the singlet state (Pt6-C94H26). Optimized bond lengths are given in Table 4.
Figure 6. Optimized geometries for Pt6 cluster on the di-vacancy-type defect in the singlet state (Pt6-C94H26). Optimized bond lengths are given in Table 4.
Molecules 17 07941 g006
Table
Table 4. Key parameters of Ptk on di-vacancy defect (C94H26) (k is 1 or 6) in Figure 6. Separations of a Pt atom from orange atoms (Pt-C(orange)), those of carbon atoms from a Pt atom except for the nearest Pt atom (other Pt-C), and those between the two orange atoms (C-C). Bond lengths are in Å. The Ebind and ΔE(Ptk) values are given in kcal/mol.
Table 4. Key parameters of Ptk on di-vacancy defect (C94H26) (k is 1 or 6) in Figure 6. Separations of a Pt atom from orange atoms (Pt-C(orange)), those of carbon atoms from a Pt atom except for the nearest Pt atom (other Pt-C), and those between the two orange atoms (C-C). Bond lengths are in Å. The Ebind and ΔE(Ptk) values are given in kcal/mol.
EbindPt-C(orange)other Pt-CC-C BondΔE(Ptk)
Pt1–106.11.999, 1.985, 1.999, 1.985––2.846, 2.846––
Pt6(i)-(A)–106.72.010, 2.001, 2.091, 2.1172.088,2.0332.810, 2.94012.7
Pt6(i)-(B)–106.32.014, 2.101, 2.109, 2.0082.087,2.1082.918, 2.92114.8
Pt6(ii)-(C)–94.52.014, 2.119, 2.088, 2.0052.206,2.314, 2.084, 2.1192.928, 2.92117.3
Pt6(ii)-(D)–85.11.994, 2.021, 2.036, 2.1142.0402.836, 2.94117.9
Molecules 17 07941 g007 550
Figure 7. Optimized geometries for Pt6 cluster on the tri-vacancy-type defect in the triplet state (Pt6-C93H26). Optimized bond lengths are given in Table 5.
Figure 7. Optimized geometries for Pt6 cluster on the tri-vacancy-type defect in the triplet state (Pt6-C93H26). Optimized bond lengths are given in Table 5.
Molecules 17 07941 g007
Table
Table 5. Key parameters of Ptk on tri-vacancy defect (C93H26) (k is 1 or 6) in Figure 7. Separations of a Pt atom from orange atoms (Pt-C(orange)) and those from the red atom (Pt-C(red)). Separations of carbon atoms from a Pt atom except for the nearest Pt atom (other Pt-C), and those between the two orange atoms (C-C). Bond lengths are in Å. The Ebind and ΔE(Ptk) values are given in kcal/mol.
Table 5. Key parameters of Ptk on tri-vacancy defect (C93H26) (k is 1 or 6) in Figure 7. Separations of a Pt atom from orange atoms (Pt-C(orange)) and those from the red atom (Pt-C(red)). Separations of carbon atoms from a Pt atom except for the nearest Pt atom (other Pt-C), and those between the two orange atoms (C-C). Bond lengths are in Å. The Ebind and ΔE(Ptk) values are given in kcal/mol.
EbindPt-C(orange)Pt-C(red)other Pt-CC-C Bond ΔE (Ptk)
Pt1–162.62.114, 2.333, 2.097, 2.6412.059––2.637, 2.641––
Pt6(i)-(A)–198.21.992, 1.959, 1.971, 1.9692.1052.0542.753, 2.83012.2
Pt6(i)-(B)–193.42.006, 1.988, 1.974, 1.9382.1442.1302.825, 2.7465.9
Pt6(ii)-(C)–188.41.971, 2.081, 1.970, 1.9692.0462.378, 2.217, 2.1072.924, 2.92220.0
Pt6(ii)-(D)–142.02.006, 2.0152.0212.130, 2.2312.71017.9
Corresponding optimized structures for the smaller Ptk clusters adsorbed (k = 2, 3, 4, and 5) are also displayed in Figure 8–Figure 10.
Molecules 17 07941 g008 550
Figure 8. Optimized geometries for Ptk cluster (k = 2~5) on the mono-vacancy-type defect in the singlet state (Ptk-C95H26). Optimized bond lengths are given in Table 6.
Figure 8. Optimized geometries for Ptk cluster (k = 2~5) on the mono-vacancy-type defect in the singlet state (Ptk-C95H26). Optimized bond lengths are given in Table 6.
Molecules 17 07941 g008
Table
Table 6. Key parameters of Ptk on mono-vacancy defect (C95H26) (k = 2 ~ 5) in Figure 8.
Table 6. Key parameters of Ptk on mono-vacancy defect (C95H26) (k = 2 ~ 5) in Figure 8.
EbindPt-C(orange) Pt-C(red) other Pt-CC-C BondΔE(Ptk)
Pt2-(A)–131.71.964, 2.0061.969––2.75529.6
Pt2-(B)–146.11.965, 1.9651.942––2.7439.0
Pt3-(A)–132.41.971, 1.9721.971––2.7150.9
Pt3-(B)–126.91.952, 1.9761.962––2.7902.5
Pt4(ii)-(A)–139.71.983, 1.9841.967––2.7597.4
Pt4(ii)-(B)–140.01.981, 1.9871.9872.0542.7437.4
Pt4(i)–128.31.953, 1.9661.9602.2272.75719.8
Pt5(ii)-(A)–140.51.976, 1.9761.980––2.7534.3
Pt5(ii)-(B)–147.81.982, 1.9822.0022.1052.7415.4
Pt5(i)–132.41.947, 1.9661.9722.211, 2.2752.7904.6
Molecules 17 07941 g009 550
Figure 9. Optimized geometries for Ptk cluster (k = 2~5) on the di-vacancy-type defect in the singlet state (Ptk-C94H26). Optimized bond lengths are given in Table 7.
Figure 9. Optimized geometries for Ptk cluster (k = 2~5) on the di-vacancy-type defect in the singlet state (Ptk-C94H26). Optimized bond lengths are given in Table 7.
Molecules 17 07941 g009
Table
Table 7. Key parameters of Ptk on di-vacancy defect (C94H26) (k = 2 ~ 5) in Figure 9.
Table 7. Key parameters of Ptk on di-vacancy defect (C94H26) (k = 2 ~ 5) in Figure 9.
EbindPt-C(orange)other Pt-CC-C BondΔE(Ptk)
Pt2-(A)–91.21.999, 2.028, 2.038, 2.066––2.831, 2.90410.6
Pt2-(B)–59.21.934, 1.935, 1.997, 1.997––2.786, 2.7407.1
Pt3-(A)–70.92.025, 2.118, 2.025, 2.1192.045, 2.0462.893, 2.8939.6
Pt3-(B)–63.21.990, 2.013, 2.125, 2.1992.081, 2.0822.791, 2.9767.0
Pt3-(C)–57.61.972, 1.987, 1.971, 1.986––2.690, 2.6901.1
Pt4(i)-(A)–84.82.034, 2.035, 2.101, 2.1002.053, 2.0542.852, 2.8536.4
Pt4(i)-(B)–75.82.004, 2.005, 2.064, 2.1342.067, 2.0662.786, 2.8964.7
Pt5(i)-(A)–93.42.005, 2.006, 2.085, 2.0992.122, 2.1382.917, 2.92312.3
Pt5(i)-(B)–81.31.991, 2.001, 2.113, 2.1232.053, 2.0822.806, 2.93213.4
Molecules 17 07941 g010 550
Figure 10. Optimized geometries for Ptk cluster (k = 2~5) on the tri-vacancy-type defect in the singlet state (Ptk-C93H26). Optimized bond lengths are given in Table 8.
Figure 10. Optimized geometries for Ptk cluster (k = 2~5) on the tri-vacancy-type defect in the singlet state (Ptk-C93H26). Optimized bond lengths are given in Table 8.
Molecules 17 07941 g010
Table
Table 8. Key parameters of Ptk on tri-vacancy defect (C93H26) (k = 2 ~ 5) in Figure 10. Separations of a Pt atom from orange atoms (Pt-C(orange)) and those from the red atom (Pt-C(red)). Separations of carbon atoms from a Pt atom except for the nearest Pt atom (other Pt-C), and those between the two orange atoms (C-C). Bond lengths are in Å. The Ebind and ΔE(Ptk) values are given in kcal/mol. Their definition was given in the text.
Table 8. Key parameters of Ptk on tri-vacancy defect (C93H26) (k = 2 ~ 5) in Figure 10. Separations of a Pt atom from orange atoms (Pt-C(orange)) and those from the red atom (Pt-C(red)). Separations of carbon atoms from a Pt atom except for the nearest Pt atom (other Pt-C), and those between the two orange atoms (C-C). Bond lengths are in Å. The Ebind and ΔE(Ptk) values are given in kcal/mol. Their definition was given in the text.
EbindPt-C(orange)Pt-C(red)other Pt-CC--C BondΔE(Ptk)
Pt2-(A)–172.31.944, 1.968, 1.988, 1.9502.038––2.839, 2.85222.0
Pt2-(B)–132.02.002, 1.9852.0092.0212.874, 1.6308.3
Pt3–167.51.948, 1.964, 1.989, 2.1242.0212.0672.999, 2.8329.8
Pt4–180.31.978, 2.068, 1.973, 1.9812.0462.0502.897, 2.7043.6
Pt5–173.01.950, 1.976, 1.934, 1.9972.099––2.840, 2.75013.9
Here Etotal(Ptk-C96–nH26) is the total energy of an optimized C96–nH26 geometry, Etotal(C96–nH26) is that of the optimized C96–nH26 geometry, and Etotal(Ptk) is that of the optimized Ptk cluster. For the DFT calculations of the binding energies, a counterpoise correction for basis set superposition error (BSSE) was included [55]. When an Ebind value has a negative sign, the binding of a Pt cluster or the Pt atom into C96–nH26 is energetically preferable. As shown in Table 3–Table 5, the calculated Ebind values in the single Pt addition are similar to those reported in [35]. These similarities verify the reliability of our DFT results.
Looking at the Ebind values, the bindings of a Ptk cluster into the sp2 surface are strongly facilitated by introducing vacancy-type defects. In fact, their stabilizing energies (–Ebind) in stable Ptk-C95H26 structures are around 150 kcal/mol. These values are much larger than the pristine cases (about 50 kcal/mol [22]). Similar enhancement in the stabilization energies was also found in the Ptk-C94H26 and Ptk-C93H26 structures. Judging from the Ebind values, reactivity of vacancy-type defects toward Pt clusters declines in the order: tri-vacancy > mono-vacancy > di-vacancy. These results suggest that the tri-vacancy defect is more suitable for binding of Pt clusters into carbon surface rather than the mono- and di-vacancy defects.
From Table 3–Table 8, we see different behaviors between the three types of defect in terms of dependences of the Ebind values on Pt cluster size. Most stable Ptk-C95H26 structures except for Pt3-C95H26 have Ebind values similar to that in Pt1-C95H26. In contrast, the Ebind values in Ptk-C94H26 and Ptk-C93H26 are deviated from those in Pt1-C94H26 and Pt1-C93H26. When a Ptk cluster binds into the di-vacancy-defect, their Ebind values are smaller than the Pt1-C94H26 value. However, these absolute values increase gradually with an increase of the cluster size, and seem to converge to the Pt1-C94H26 value at k = 6. In the Ptk-C93H26 cases, the Ebind values, being around 180 kcal/mol, are always larger than the Pt1-C93H26 value.
To understand the energetics in the optimized Ptk-C96–nH26 structures (Table 3–Table 8), let us first look at in detail geometrical features of Pt1-C96–nH26. These key geometrical parameters in the Pt1-C96-nH26 configurations (lengths of newly formed Pt-C bonds and of lengthening CC bonds) are listed in Table 3–Table 5.
In these tables, we can distinguish two types of the formed Pt-C bond, by whether a Pt atom binds into orange or red atoms. When the single Pt atom binds into the mono-vacancy defect, it inserts between the orange atoms in the five-membered ring, and then two P-C(orange) bonds are formed newly. As a result of the Pt addition, the separation between the orange atoms lengthens from 1.754 to 2.764 Å. At the same time, the Pt atom also coordinates to the unsaturated red atom. The binding Pt atom lifts from the sp2 surface, because the hole is not large enough to accommodate the Pt atom. Similar Pt lifting can be seen in the Pt1-C94H26 configuration, where the Pt atom inserts between orange atoms in both five-membered rings, and it breaks the connections. The degree of Pt lifting in Pt1-C94H26 is less significant than that in Pt1-C95H26, due to relatively larger hole in C94H26.
In contrast, the hole of C93H26, surrounded by ten carbon atoms, can house the Pt atom, and therefore the binding Pt atom is on the sp2 surface. Then, four Pt-C bonds are formed, accompanying the cleavage of the bonds between orange atoms in the five-membered rings. Moreover, the Pt binding into the unsaturated C atom was also seen. When a Ptk cluster binds into a vacancy-type defect, slightly longer separations of the nearest Pt atom from reactive (orange and red) atoms were found. Despite the stabilization operated between a Ptk cluster and C96–nH26, slightly longer Pt-C bonds imply weakening interactions of the nearest Pt atom from the reactive carbon atoms compared with Pt1-C96–nH26 case.
Compensating the weakening of the interactions, remaining Pt atoms of a clusters are additionally bound to carbon atoms of a defective site to maximize the Pt-C interactions. Then their clusters are more or less deformed from the most stable configuration in the gas-phase [56,57,58]. The degree of cluster deformation was estimated by using ΔE(Ptk), defined as [E(Ptk on surface) – E(Ptk)], where E(Ptk on surface) is the total energy of Ptk cluster taken from an optimized Ptk-C96–nH26 structure and E(Ptk) is that of the optimized geometry for the bare Ptk cluster. Positive ΔE(Ptk) values in Table 3–Table 8 suggest destabilization from cluster deformation upon the interactions with a vacancy-type defect. Although we cannot find a clear correlation between Ebind and ΔE(Ptk) values, the balance between the stabilization from the Pt-C bond formation and the destabilization from the cluster deformation is a key in determining the stability. From Figure 5–Figure 7 and Figure 11, we found clear differences between Ptk-C93H26 and Ptk-C95H26 (Ptk-C94H26) in terms of the number of Pt atoms binding directly into orange atoms in the defective site to cleave connections between adjacent orange atoms.
Molecules 17 07941 g011 550
Figure 11. Optimized geometries for the singlet Pt atom on the mono-, di-, and tri-vacancy-type defects in the triplet state (Pt1-C95H26, Pt1-C94H26, and Pt1-C93H26, respectively). Optimized bond lengths are given in Table 3–Table 5.
Figure 11. Optimized geometries for the singlet Pt atom on the mono-, di-, and tri-vacancy-type defects in the triplet state (Pt1-C95H26, Pt1-C94H26, and Pt1-C93H26, respectively). Optimized bond lengths are given in Table 3–Table 5.
Molecules 17 07941 g011
In Ptk-C95H26 (Ptk-C94H26), one Pt atom participates in cleaving the orange connection(s), irrespective of the cluster size. In the tri-vacancy cases containing larger ten-membered ring, two Pt atoms in a cluster bind to the defective site to split two orange connections. The accommodation of two Pt atoms cannot be seen in the Pt1-C94H26 structure, and thus the significant enhanced stabilization in Ptk-C93H26 is understandable. Moreover the acceptability of the ten-membered-ring to trap Pt atoms differentiates C93H26 from C95H26 and C94H26 in terms of their reactivity toward Ptk clusters. Due to the strong interactions between a Pt6 cluster and a vacancy-type defect, we can see unique orbital features, which cannot be seen in C96–nH26 (Figure 12).
In fact, 5d(Pt)-based orbitals, given by blue bars in Figure 12, appear in the frontier orbital regions of the Pt6-C94H26 and Pt6-C93H26 configurations. As the most striking case, we can see in Figure 13 that the Pt6(i)-C94H26(B) configuration has the HOMO and LUMO consisting of 5d(Pt) orbitals. On the other hand, levels of 5d(Pt)-based orbitals in the Pt6-C95H26 strongly depend on their cluster-shape. In the Pt6(i)-C95H26(A) and Pt6(ii)-C95H26(C) configurations, such 5d(Pt)-based orbital lies larger than 1.3 eV above the LUMO, whereas the LUMO+1 consists of 5d(Pt)-based orbitals in the other Pt6-C95H26 configurations.

2.3.2. Triplet State

As shown in Figure 12, the all optimized geometries in the single state have relatively small HOMO-LUMO gaps (0.29 ~ 0.40 eV). Thus, higher spin states can be energetically stable relative to the singlet states. Along the assumption, we obtained their triplet states, and estimated the energy difference between the two spin states (ΔEstate(Pt6-C96–nH26)), as tabulated in Table 9.
Molecules 17 07941 g012 550
Figure 12. Orbital energies (eV) in the frontier orbital region of the optimized Pt6-C95H26, Pt6-C94H26, and Pt6-C93H26 configurations whose structures are given in Figure 5–Figure 7. The HOMO-LUMO gaps are given. Orbitals originated from 5d(Pt) orbitals are denoted by blue bars, and those with no or less 5d(Pt) orbital contribution are denoted by black bars.
Figure 12. Orbital energies (eV) in the frontier orbital region of the optimized Pt6-C95H26, Pt6-C94H26, and Pt6-C93H26 configurations whose structures are given in Figure 5–Figure 7. The HOMO-LUMO gaps are given. Orbitals originated from 5d(Pt) orbitals are denoted by blue bars, and those with no or less 5d(Pt) orbital contribution are denoted by black bars.
Molecules 17 07941 g012
Molecules 17 07941 g013 550
Figure 13. Frontier orbitals (the HOMO and LUMO) in the Pt6(i)-C94H26(B) configuration (Figure 6) are given as a representative Pt6-C96–nH26 configuration.
Figure 13. Frontier orbitals (the HOMO and LUMO) in the Pt6(i)-C94H26(B) configuration (Figure 6) are given as a representative Pt6-C96–nH26 configuration.
Molecules 17 07941 g013
Table
Table 9. Energy difference between the singlet and triplet states in the Pt6-C96–nH26 configurations (ΔEstate in kcal/mol) a.
Table 9. Energy difference between the singlet and triplet states in the Pt6-C96–nH26 configurations (ΔEstate in kcal/mol) a.
Pt6(i)-(A)Pt6(i)-(B)Pt6(ii)Pt6(iii)
Clusters on C95H26–14.3–17.8–14.0b
Clusters on C94H26–23.3–16.2–23.3–38.5
Clusters on C93H26–15.0–23.7–13.7b
a ΔEstate(Pt6-C96–nH26) = Etotal(triplet state) – Etotal(singlet state). Negative ΔEstate values indicate that the triplet state of a Pt6-C96–nH26 configuration is energetically stable relative to the singlet state. ΔEstate in kcal/mol; b we could not obtain the optimized geometry in the triplet state.
Table 9 shows that their triplet spin states are energetically stable relative to the singlet states, as expected. According to DFT calculations, most Pt6-C96–nH26 configurations have radical Pt6 clusters on defective graphene patches As representative cases, spin density distributions on the Pt6(i) cluster binding into C95H26 or C94H26 are shown in Figure 14. Figure 14 shows that substantial spin densities appear on the adsorbed clusters in the three Pt6-C96–nH26 configurations (Pt6(i)-C95H26(B), and two Pt6(i)-C94H26 configurations), whereas the Pt6(i)-C95H26(A) configuration does not have such radical cluster due to the absence of 5d(Pt)-based frontier orbitals (Figure 12). In the three configurations with radical clusters, a variety of the spin density distributions was found. In the Pt6(i)-C95H26(B) and Pt6(i)-C94H26(B), the spin density distributions are strongly deviated from those in the bare Pt6 cluster, whereas the Pt6(i)-C94H26(B) has similar distributions. Furthermore, we found a relationship between spin densities of Pt clusters and those on the defective sp2 surface. When spin density distributions of the Pt cluster of a Pt6-C96-nH26 configuration are (not) similar to the bare case, spin densities are (not) delocalized over the carbon surface. Similar tendencies were found in spin density distributions of the stable triplet Pt13-C96–nH26 configurations [59], because radical Pt clusters exist on the defective graphene patches as shown in Figure 15.
Molecules 17 07941 g014 550
Figure 14. Spin density distributions of representative Pt6-C96–nH26 configurations (Pt6(i) cluster on C95H26 or C94H26 in two binding fashions, displayed in Figure 5 and Figure 6). Isosurface α- and β-spins are given by pink and blue, respectively.
Figure 14. Spin density distributions of representative Pt6-C96–nH26 configurations (Pt6(i) cluster on C95H26 or C94H26 in two binding fashions, displayed in Figure 5 and Figure 6). Isosurface α- and β-spins are given by pink and blue, respectively.
Molecules 17 07941 g014
Finally, let us compare the spin density distributions on Pt6 clusters on defective sp2 surfaces (Figure 14) with the pristine case (Figure 2). As mentioned above, spin density maps on Pt6 cluster on pristine C96H26 are similar to those in the bare cluster. However, such similarity cannot be always found in the defective graphene cases. The different tendencies come from whether underlying carbon atoms have radical characters or not in Figure 2 and Figure 4. In the pristine patch, underlying carbon atoms do not have radical characters, and they cannot perturb spin density distributions of the absorbed cluster. In contrast, unpaired electrons exist on underlying carbon atoms in the defective sites, and thus spin density distributions of the adsorbed clusters are modulated by the Pt-C interactions. Therefore, perturbation of the radical sp2 surface by introduction of vacancy-type defects can change characters of adsorbed Pt clusters.
Molecules 17 07941 g015 550
Figure 15. Spin density distributions of representative Pt13-C96–nH26 configurations. Isosurface α- and β-spins are given by pink and blue, respectively.
Figure 15. Spin density distributions of representative Pt13-C96–nH26 configurations. Isosurface α- and β-spins are given by pink and blue, respectively.
Molecules 17 07941 g015
The DFT findings are important in the catalytic activity of Pt clusters on sp2 surface, because the supported clusters can serve as active site for catalytic reactions. For example, if radical Pt clusters exist on carbon surface, it can cleave the H-H bond of hydrogen molecules in a homolytic manner. Otherwise, the H-H bond is activated by the clusters via a non-radical mechanism in [22]. Thus, the DFT calculations propose that one can change chemical reactivity of Pt clusters on graphene patches by introducing of vacancy-type defects on the surface.

3. Experimental

To investigate interactions of Pt clusters with C96H26, we carried out DFT calculations implemented in the Gaussian 03 and 09 program packages [60,61]. Adsorbed clusters that we considered consist of one, six, and thirteen Pt atoms. To perform the calculations, a hybrid Hartree–Fock/density functional theory method, B3LYP [62,63,64,65,66] was used. The B3LYP method consists of the Slater exchange, the Hartree–Fock exchange, the exchange functional of Becke [62,63,64], the correlation functional of Lee, Yang, and Parr (LYP) [65], and the correlation functional of Vosko, Wilk, and Nusair (VWN) [66]. In general the hybrid B3LYP method has been reported to provide excellent descriptions of various properties. The Gaussian-type basis set we used for the C and H atoms is 6-31G* [67], and that for the Pt atom is the quasi-relativistic effective core potential RECP and valence basis sets recommended by Stuttgart group (SDD) [68]. The SDD RECP is adjusted to total valence energies of a multitude of atomic references states, which are quantum mechanical observables [68]. As indicated in the previous papers [69,70,71,72,73,74,75,76], the B3LYP/6-31G* calculations correctly reproduce experimental data for C60, especially its IR and Raman vibrational frequencies. According to the theoretical report by Nova et al. [77], the method of our choice (B3LYP/SDD + 6-31G*) is appropriate to reproduce experimental values in terms of Pt-C bonds in Pt complexes [77]. The computational method is also suitable to study transition metals adsorbed on graphene. In fact, Pt-C bond lengths obtained from the B3LYP/SDD + 6-31G* calculations fall in the range reported from other theoretical reports [24,30].

4. Conclusions

Density functional theory (DFT) B3LYP calculations were employed to investigate the adsorption of Ptk cluster (k is 1–6, and 13) into a nanometer-size graphene patch (C96H26) with or without vacancy-type defects. According to the DFT calculations, removing a few carbon atoms (n) from C96H26 results in the formation of five-membered rings as well as coordinatively unsaturated carbon atoms. Introduction of a vacancy-type defect on C96H26 strongly affects spin density distributions in its stable triplet state. Although spin densities appear only on edge carbon atoms of the triplet C96H26 structure, defective graphene patches have radical carbon atoms at the center where the reactive carbon atoms exist. These spin density distributions differentiate characters of Pt clusters adsorbed on defective graphene patches from those on the pristine. According to the DFT calculations, spin density maps of Pt clusters on C96H26 are similar to those of the corresponding bare clusters. In contrast, Pt clusters interact strongly with radical carbon atoms in defective graphene patches, and thus spin density distributions of the adsorbed Pt clusters are usually deviated from the bare cases. Consequently, DFT calculations propose that characters of Pt clusters adsorbed on the sp2 carbon surface can be modulated by introducing vacancy-type defects.

Acknowledgments

Support by the Japan Society for the Promotion of Science (JSPS) for T.Y., and by the Ministry of Culture, Sports, Science and Technology of Japan (MEXT) for T.Y. and H.K.

References and Notes

  1. Serp, P.; Figueiredo, J.L. Carbon Materials for Catalysis; Wiley: Hoboken, NJ, USA, 2009. [Google Scholar]
  2. Auer, E.; Freund, A.; Pietsch, J.; Tacke, T. Carbons as supports for industrial precious metal catalysts. Appl. Catal. A 1998, 173, 259–271. [Google Scholar] [CrossRef]
  3. Lee, S.-A.; Park, K.-W.; Choi, J.-H.; Kwon, B.-K.; Sung, Y.-E. Nanoparticle synthesis and electrocatalytic activity of Pt alloys for direct methanol fuel cells. J. Electrochem. Soc. 2002, 149, A1299–A1304. [Google Scholar]
  4. Lu, Q.; Yang, L.; Zhuang, L.; Lu, J. Anodic activation of PtRu/C catalysts for methanol oxidation. J. Phys. Chem. B 2005, 109, 1715–1722. [Google Scholar] [CrossRef]
  5. Holstein, W.L.; Rosenfeld, H.D. In-situ X-ray absorption spectroscopy study of Pt and Ru chemistry during methanol electrooxidationt. J. Phys. Chem. B 2005, 109, 2176–2186. [Google Scholar] [CrossRef]
  6. Wang, L.L.; Khare, S.V.; Chirita, V.; Johnson, D.D.; Rockett, A.A.; Frenkel, A.I.; Mack, N.H.; Nuzzo, R.G. Origin of bulklike structure and bond length disorder of Pt37 and Pt6Ru31 clusters on carbon: Comparison of theory and experiment. J. Am. Chem. Soc. 2006, 128, 131–142. [Google Scholar]
  7. Bowker, M. Surface science—The going rate for catalysts. Nat. Mater. 2002, 1, 205–206. [Google Scholar] [CrossRef]
  8. Campbell, C.T.; Parker, S.C.; Starr, D.E. The effect of size-dependent nanoparticleenergetics on catalyst sintering. Science 2002, 298, 811–814. [Google Scholar] [CrossRef]
  9. Kim, S.J.; Park, Y.J.; Ra, E.J.; Kim, K.K.; An, K.H.; Lee, Y.H.; Choi, J.Y.; Park, C.H.; Doo, S.K.; Park, M.H.; et al. Defect-induced loading of Pt nanoparticles on carbon nanotubes. Appl. Phys. Lett. 2007, 90, 023114. [Google Scholar]
  10. Maiti, A.; Ricca, A. Metal-nanotube interactions—Binding energies and wetting properties. Chem. Phys. Lett. 2004, 395, 7–11. [Google Scholar] [CrossRef]
  11. Okamoto, Y. Density-functional calculations of icosahedral M13 (M = Pt and Au) clusters on graphene sheets and flakes. Chem. Phys. Lett. 2006, 420, 382–386. [Google Scholar] [CrossRef]
  12. Chi, D.H.; Cuong, N.T.; Tuan, N.A.; Kim, Y.-T.; Bao, H.T.; Mitani, T.; Ozaki, T.; Nagao, H. Electronic structures of Pt clusters adsorbed on (5,5) single wall carbon nanotube. Chem. Phys. Lett. 2006, 432, 213–217. [Google Scholar]
  13. Kong, K.; Choi, Y.; Ryu, B.-H.; Lee, J.; Chang, H. Investigation of metal/carbon-related materials for fuel cell applications by electronic structure calculations. Mater. Sci. Eng. C 2006, 26, 1207–1210. [Google Scholar]
  14. Acharya, C.K.; Turner, C.H. Stabilization of platinum clusters by substitutional boron dopants in carbon supports. J. Phys. Chem. B 2006, 110, 17706–17710. [Google Scholar] [CrossRef]
  15. Acharya, C.K.; Turner, C.H. Effect of an electric field on the adsorption of metal clusters on boron-doped carbon surfaces. J. Phys. Chem. C 2007, 111, 14804–14812. [Google Scholar] [CrossRef]
  16. Acharya, C.K.; Sullican, D.I.; Turner, C.H. Characterizing the interaction of Pt and PtRu clusters with boron-doped, nitrogen-doped, and activated carbon: Density functional theory calculations and parameterization. J. Phys. Chem. C 2008, 112, 13607–13622. [Google Scholar]
  17. Cuong, N.T.; Fujiwara, A.; Mitani, T.; Chi, D.H. Effects of carbon supports on Pt nano-cluster catalyst. Comput. Mater. Sci. 2008, 44, 163–166. [Google Scholar]
  18. Wang, J.; Lv, Y.; Li, X.; Dong, M. Point-defect mediated bonding of Pt clusters on (5,5) carbon nanotubes. J. Phys. Chem. C 2009, 113, 890–893. [Google Scholar]
  19. Akturk, O.U.; Tomak, M. AunPtn clusters adsorbed on graphene studied by first-principles calculations. Phys. Rev. B 2009, 80, 085417. [Google Scholar] [CrossRef]
  20. Cuong, N.T.; Sugiyama, A.; Fujiwara, A.; Mitani, T.; Chi, D.H. Density functional study of Pt4 clusters adsorbed on a carbon nanotube support. Phys. Rev. B 2009, 79, 235417. [Google Scholar]
  21. Lu, Y.-H.; Zhou, M.; Zhang, C.; Feng, Y.-P. Metal-embedded graphene: A possible catalyst with high activity. J. Phys. Chem. C 2009, 113, 20156–20160. [Google Scholar]
  22. Yumura, T.; Kimura, K.; Kobayashi, H.; Tanaka, R.; Okumura, N.; Yamabe, T. The use of nanometer-sized hydrographene species for support material for fuel cell electrode catalysts; a theoretical proposal. Phys. Chem. Chem. Phys. 2009, 11, 8275–8284. [Google Scholar]
  23. Okazaki-Maeda, K.; Morikawa, Y.; Tanaka, S.; Kohyama, M. Structures of Pt clusters on graphene by first-principles calculations. Surf. Sci. 2010, 604, 144–154. [Google Scholar] [CrossRef]
  24. Cabria, I.; Lopez, M.J.; Alonso, J.A. Theoretical study of the transition from planar to three-dimensional structures of palladium clusters supported on graphene. Phys. Rev. B 2010, 81, 035403. [Google Scholar] [CrossRef]
  25. Zhou, M.; Zhang, A.; Dai, Z.; Feng, Y.P.; Zhang, C. Strain-enhanced stabilization and catalytic activity of metal nanoclusters on graphene. J. Phys. Chem. C 2010, 114, 16541–16546. [Google Scholar]
  26. Zhou, M.; Zhang, A.; Dai, Z.; Zhang, C.; Feng, Y.P. Greatly enhanced adsorption and catalytic activity of Au and Pt clusters on defective graphene. J. Chem. Phys. 2010, 132, 194704. [Google Scholar]
  27. Vedala, H.; Sorescu, D.C.; Kotchey, G.P.; Star, A. Chemical sensitivity of graphene edges decorated with metal nanoparticles. NanoLett. 2011, 11, 2342–2347. [Google Scholar]
  28. Błoński, P.; Hafner, J. Geometric and magnetic properties of Pt clusters supported on graphene: Relativistic density-functional calculations. J. Chem. Phys. 2011, 134, 154705. [Google Scholar] [CrossRef]
  29. Lim, D.-H.; Wilcox, J. DFT-Based Study on oxygen adsorption on defective graphene-supported Pt nanoparticles. J. Phys. Chem. C 2011, 115, 22742–22747. [Google Scholar]
  30. Fampiou, I.; Ramasubramaniam, A. Binding of Pt nanoclusters to point defects in graphene: Adsorption, morphology, and electronic structure. J. Phys. Chem. C 2012, 116, 6543–6555. [Google Scholar]
  31. Lim, D.-H.; Wilcox, J. Mechanisms of the oxygen reduction reaction on defective graphene-supported Pt nanoparticles from first principles. J. Phys. Chem. C 2012, 116, 3653–3660. [Google Scholar] [CrossRef]
  32. Iijima, S. High resolution electron microscopy of phase objects: Observation of small holes and steps on graphite crystals. Optik 1977, 47, 437–452. [Google Scholar]
  33. Hashimoto, A.; Suenaga, K.; Gloter, A.; Urita, K.; Iijima, S. Direct evidence for atomic defects in graphene layers. Nature 2004, 430, 870–873. [Google Scholar]
  34. Banhart, F.; Kotakoski, J.; Krasheninnikov, A.V. Structural defects in graphene. ACS Nano 2011, 5, 26–41. [Google Scholar] [CrossRef]
  35. Wang, H.; Wang, Q.; Cheng, Y.; Li, K.; Yao, Y.; Zhang, Q.; Dong, C.; Wang, P.; Schwingenschlögl, U.; Yang, W.; et al. Doping monolayer graphene with single atom substitutions. NanoLett. 2012, 12, 141–144. [Google Scholar] [CrossRef]
  36. Fujita, M.; Wakabayashi, K.; Nakada, K.; Kusakabe, K. Peculiar localized state at zigzag graphite edge. J. Phys. Soc. Jpn. 1996, 65, 1920–1923. [Google Scholar]
  37. Whangbo, M.-H.; Hoffmann, R.; Woodward, R.B. Conjugated one and 2-dimenstional polymers. Proc. R. Soc. Lond. A Math. Phys. Sci. 1979, 366, 23–46. [Google Scholar] [CrossRef]
  38. Yamabe, T.; Tanaka, K.; Ohzeki, K.; Yata, S. Electronic-structure of polyacenasene—A one-dimensional graphite. Solid State Commun. 1982, 44, 823–825. [Google Scholar]
  39. Kertesz, M.; Hoffmann, R. Higher-order Peierls distortion of one-dimensional carbon skeletons. Solid State Commun. 1983, 47, 97–102. [Google Scholar] [CrossRef]
  40. Stein, S.E.; Brown, B.L. П-electron properties of large condensed polyaromatic hydrocarbons. J. Am. Chem. Soc. 1987, 109, 3721–3729. [Google Scholar]
  41. Yoshizawa, K.; Yahara, K.; Tanaka, K.; Yamabe, T. Bandgap oscillation in polyphenanthrenes. J. Phys. Chem. B 1998, 102, 498–506. [Google Scholar]
  42. Nakada, K.; Fujita, M.; Dresselhaus, G.; Dresselhaus, M.S. Edge state in graphene ribbons: Nanometer size effect and edge shape dependence. Phys. Rev. B 1996, 54, 17954–17961. [Google Scholar]
  43. Marsh, H.; Martinez-Escandell, M.; Rodriguez-Reinoso, F. Semi-cokes from pitch pyrolysis: Mechanisms and kinetics. Carbon 1999, 37, 363–390. [Google Scholar] [CrossRef]
  44. Kusakabe, K.; Maruyama, M. Magnetic nanographite. Phys. Rev. B 2003, 67, 092406. [Google Scholar] [CrossRef]
  45. Wakabayashi, K, Sigrist; Wakabayashi, K; Sigrist, M; Fujita, M. Spin wave mode of edge-localized magnetic states in nanographite zigzag ribbons. J. Phys. Soc. Jpn. 1998, 67, 2089–2093. [Google Scholar]
  46. Hod, O.; Barone, V.; Peralta, J.E.; Scuseria, G.E. Enhanced half-metallicity in edge-oxidized zigzag graphene nanoribbons. Nano Lett. 2007, 7, 2295–2299. [Google Scholar]
  47. Zheng, H.; Duley, W. First-principles study of edge chemical modifications in graphene nanodots. Phys. Rev. B 2008, 78, 045421. [Google Scholar]
  48. Lee, G.; Cho, K. Electronic structures of zigzag graphenenanoribbons with edge hydrogenation and oxidation. Phys. Rev. B 2009, 79, 165440. [Google Scholar]
  49. Hu, X.Y.; Tian, H.W.; Zheng, W.T.; Yu, S.S.; Qiao, L.; Qu, C.Q.; Jiang, Q. Metallic-semiconducting phase transition of the edge-oxygenated armchair graphene nanoribbons. Chem. Phys. Lett. 2010, 501, 64–67. [Google Scholar] [CrossRef]
  50. Yumura, T.; Kobayashi, H.; Yamabe, T. Roles of radical characters of pristine and nitrogen-substituted hydrographene in dioxygen bindings. J. Chem. Phys. 2010, 133, 174703. [Google Scholar]
  51. Yumura, T.; Kobayashi, H.; Yamabe, T. Energetics of dioxygen binding into graphene patches with various sizes and shapes. Sci. China Chem. 2012, 55, 787–795. [Google Scholar]
  52. Kua, J.; Goddard, W.A., III. Chemisorption of organics on platinum. 1. The interstitial electron model. J. Phys. Chem. B 1998, 102, 9481–9491. [Google Scholar]
  53. Lin, X.; Ramer, N.J.; Rappe, A.M.; Hass, K.C.; Schneider, W.F.; Trout, B.L. Effect of particle size on the adsorption of O and S atoms on Pt: A density-functional theory study. J. Phys. Chem. B 2001, 105, 7739–7747. [Google Scholar]
  54. Xiao, L.; Wang, L. Structures of Platinum clusters: planar or spherical? J. Phys. Chem. A 2004, 108, 8605–8614. [Google Scholar] [CrossRef]
  55. Boys, S.F.; Bernardi, F. The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors. Mol. Phys. 1970, 19, 553–566. [Google Scholar] [CrossRef]
  56. Yumura, T.; Takeuchi, M.; Kobayashi, H.; Kuroda, Y. Effect of ZSM-5 zeolite confinement on reaction intermediates during dioxygen activation by enclosed dicopper cations. Inorg. Chem. 2009, 48, 508–517. [Google Scholar]
  57. Wu, J.; Ong, S.W.; Kang, H.C.; Tok, E.S. Hydrogen adsorption on mixed platinum and nickel nanoclusters; The influence of cluster composition and graphene support. J. Phys. Chem. C 2010, 114, 21252–21261. [Google Scholar]
  58. Yamashita, H.; Yumura, T. The role of weak bonding in determining the structure of thiopheneoligomers inside carbon nanotubes. J. Phys. Chem. C 2012, 116, 9681–9690. [Google Scholar]
  59. The Pt13-C95H26, Pt13-C94H26, and Pt13-C93H26 structures in the triplet state are more stable by 20.9, 20.3, and 9.5 kcal/mol relative to the corresponding singlet states, respectively. Similarly, the triplet state in bare Pt13 cluster is more stable by 13.7 kcal/mol relative to the singlet state.
  60. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Montgomery, J.A., Jr.; Vreven, T.; Kudin, K. N.; Burant, J.C.; et al. Gaussian 03; Gaussian, Inc.: Pittsburgh, PA, USA, 2003. [Google Scholar]
  61. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09; Gaussian, Inc.: Wallingford, CT, USA, 2009. [Google Scholar]
  62. Becke, A.D. Density-functional exchange-energy approximation with correct asymptotic behavior. Phys. Rev. A 1988, 38, 3098–3100. [Google Scholar] [CrossRef]
  63. Becke, A.D. Density-functional thermochemistry III: The Role of exact exchange. J. Chem. Phys. 1993, 98, 5648–5652. [Google Scholar] [CrossRef]
  64. Stephens, P.J.; Devlin, F.J.; Chabalowski, C.F.; Frisch, M.J. Ab initio calculation of vibrational absorption and circular dichroism spectra using density functional force fields. J. Phys. Chem. 1994, 98, 11623–11627. [Google Scholar] [CrossRef]
  65. Lee, C.; Yang, W.; Parr, R.G. Development of the Colle-Salvetti correlation-energy formula into functional of the electron-density. Phys. Rev. B 1988, 37, 785–789. [Google Scholar]
  66. Vosko, S.H.; Wilk, L.; Nusair, M. Accurate spin-dependent electron liquid correlation energies for local spin-density calculations—A critical analysis. Can. J. Phys. 1980, 58, 1200–1211. [Google Scholar]
  67. Hehre, W.J.; Ditchfield, R.; Pople, J.A. Self-consistent molecular orbital methods. XII. Further extensions of Gaussian-type basis sets for use in molecular orbital studies of organic molecules. J. Chem. Phys. 1972, 56, 2257–2261. [Google Scholar] [CrossRef]
  68. Andrae, D.; Häuβermann, U.; Dolg, M.; Stoll, H.; Preuβ, H. Energy-adjusted ab initio pseudopotentials for the second and third row transition elements. Theor. Chim. Acta 1999, 77, 123–141. [Google Scholar]
  69. Yumura, T.; Kertesz, M.; Iijima, S. Confinement effects on site-preferences for cycloadditions into carbon nanotubes. Chem. Phys. Lett. 2007, 444, 155–160. [Google Scholar]
  70. Yumura, T.; Kertesz, M. Roles of conformational restrictions of a bismalonate in the interactions with a carbon nanotube. J. Phys. Chem. C 2009, 113, 14184–14194. [Google Scholar] [CrossRef]
  71. Yumura, T.; Hirahara, K.; Bandow, S.; Yoshizawa, K.; Iijima, S. A theoretical study on the geometrical features of finite-length carbon nanotubes capped with fullerene hemisphere. Chem. Phys. Lett. 2004, 386, 38–43. [Google Scholar]
  72. Yumura, T.; Bandow, S.; Yoshizawa, K.; Iijima, S. The role of fullerene hemispheres in determining structural features of finite-length carbon nanotubes. J. Phys. Chem. B 2004, 108, 11426–11434. [Google Scholar]
  73. Yumura, T.; Nozaki, D.; Bandow, S.; Yoshizawa, K.; Iijima, S. End-cap effects on vibrational structures of finite-length carbon nanotubes. J. Am. Chem. Soc. 2005, 127, 11769–11776. [Google Scholar]
  74. Yumura, T.; Kertesz, M.; Iijima, S. Local modifications of single-wall carbon nanotubes induced by bond formation with encapsulated fullerenes. J. Phys. Chem. B 2007, 111, 1099–1109. [Google Scholar] [CrossRef]
  75. Yumura, T.; Sato, Y.; Suenaga, K.; Urita, K.; Iijima, S. Which do endohedral Ti2C80metallofullerenes prefer energetically: Ti2@C80 or Ti2C2@C78? A theoretical study. J. Phys. Chem. B 2005, 109, 20251–20255. [Google Scholar]
  76. Yumura, T.; Sato, Y.; Suenaga, K.; Urita, K.; Iijima, S. Gate effect of vacancy-type defect of fullerene cages on metal-atom migrations in metallofullerenes. NanoLett. 2006, 6, 1389–1395. [Google Scholar] [CrossRef]
  77. Nova, A.; Erhardt, S.; Jasim, N.A.; Perutz, R.N.; Macgregor, S.A.; McGrady, E.; Whitwood, A.C. Competing C−F activation pathways in the reaction of Pt(0) with fluoropyridines: Phosphine-assistance versus oxidative addition. J. Am. Chem. Soc. 2008, 130, 15499–15511. [Google Scholar]
  • Sample Availability: Not Available.

Share and Cite

MDPI and ACS Style

Yumura, T.; Awano, T.; Kobayashi, H.; Yamabe, T. Platinum Clusters on Vacancy-Type Defects of Nanometer-Sized Graphene Patches. Molecules 2012, 17, 7941-7960. https://doi.org/10.3390/molecules17077941

AMA Style

Yumura T, Awano T, Kobayashi H, Yamabe T. Platinum Clusters on Vacancy-Type Defects of Nanometer-Sized Graphene Patches. Molecules. 2012; 17(7):7941-7960. https://doi.org/10.3390/molecules17077941

Chicago/Turabian Style

Yumura, Takashi, Tatsuya Awano, Hisayoshi Kobayashi, and Tokio Yamabe. 2012. "Platinum Clusters on Vacancy-Type Defects of Nanometer-Sized Graphene Patches" Molecules 17, no. 7: 7941-7960. https://doi.org/10.3390/molecules17077941

APA Style

Yumura, T., Awano, T., Kobayashi, H., & Yamabe, T. (2012). Platinum Clusters on Vacancy-Type Defects of Nanometer-Sized Graphene Patches. Molecules, 17(7), 7941-7960. https://doi.org/10.3390/molecules17077941

Article Metrics

Back to TopTop