Pseudocarbynes: Charge-Stabilized Carbon Chains - The Journal of

Apr 14, 2016 - Atomistic simulations of the equation of state and hybridization of liquid carbon at a temperature of 6000 K in the pressure range of 1...
0 downloads 11 Views 3MB Size
Letter pubs.acs.org/JPCL

Pseudocarbynes: Charge-Stabilized Carbon Chains Pilarisetty Tarakeshwar,*,† Peter R. Buseck,†,‡ and Harold W. Kroto§ †

School of Molecular Sciences, Arizona State University, Tempe, Arizona 85287-1604, United States School of Earth and Space Exploration, Arizona State University, Tempe, Arizona 85287-6004, United States § Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States ‡

S Supporting Information *

ABSTRACT: Carbyne is the long-sought linear allotrope of carbon. Despite many reports of solid carbyne, the evidence is unconvincing. A recent report of supposed carbyne shows gold clusters in transmission electron microscopy (TEM) images. In order to determine the effects of such clusters, we performed ab initio calculations of uncapped and capped linear carbon chains and their complexes with gold clusters. The results indicate that gold dramatically alters the electron densities of the CC bonds. The resulting charge-stabilization of the carbon chains leads to pseudocarbynes. These findings are corroborated in calculations of the structures of crystals containing isolated carbon chains and those intercalated with gold clusters. Calculated Raman spectra of these pseudocarbynes with gold clusters are in better agreement with experiment than calculated spectra of isolated carbon chains. The current work opens the way toward the design and development of a new class of metal-intercalated carbon compounds.

T

can be physically isolated and their interesting properties experimentally harnessed.6,40 Although the experimental evidence for such bundles is questionable, finite-sized carbon chains have been observed: (a) in supersonic-beam-deposited carbon films,41,42 (b) in solid argon at 36 K,43 (c) in break junctions involving graphene,44,45 (d) as chains encapsulated in carbon nanotubes,46,47 and (e) on gold/silver nanoparticles or films.48−50 However, these carbon chains are short-lived (a few seconds) and have not been physically isolated Carbon chains in the solid state were reportedly synthesized using gold either as a substrate or catalyst.51,52 High catalytic activity has been correlated with bilayer clusters containing only ∼10 gold atoms, with dimensions of ∼0.5 nm,53 and stabilized in carbon nanotubes or on sheets.54,55 Although there have been investigations of the interactions of gold clusters with ethene, ethyne, carbon nanotubes, and graphene,56−61 the influence of such clusters on the structure and properties of the carbon chains has not been investigated. Such a study becomes important because of a recent report on the isolation of crystalline carbon chains for which the transmission electron microscopy (TEM) image of the samples shows gold nanoparticles or clusters embedded on the carbon nanorods.52 In order to investigate the structure and electronic properties of finite-sized carbon chains in the presence of small gold clusters, we carried out calculations at the density functional level of theory (DFT) using the Becke gradient-corrected exchange functional and Lee−Yang−Parr correlation functional with three parameters (B3LYP).62−65 The cc-pVDZ basis set

he carbon allotropes have long been the subject of intensive investigation. Interest increased following the discovery of the structures and fascinating properties of fullerenes, carbon nanotubes, and graphene.1−6 However, the simplest allotrope of carbon, the one-dimensional carbyne, has defied attempts to synthesize and characterize it in the solid state. Since its first reported synthesis 50 years ago, the quest for carbyne has been marred by repeated synthetic and spectroscopic setbacks, and the trail is littered with red herrings.7−16 The term carbyne has been used both to represent a compound whose structure contains a carbon atom with three nonbonding electrons and for the crystalline carbon allotrope consisting of bundles of linear carbon chains,17−19 the usage we adopt in this paper. The chains can either have alternating single and triple bonds (CC)n or only double bonds (CC)n.17−19 In contrast, polyynes and oligoynes consist of linear carbon chains end-capped with another element, commonly hydrogen, or organic functional groups.17,18 On the basis of the definition of allotrope,20,21 carbyne contains only carbon with no other element or functional group as essential constituents. Moreover, carbyne would be explosively reactive in the condensed phase because the unpaired electrons at either end of each chain would lead to rapid cross-linking or polymerization.15,16,22 Nonetheless, the quests continue. In the attempt to synthesize carbyne, a common approach is to first produce carbyne-like materials by adding massive end groups to linear carbon chains, with the goal of hindering their reactivity.22−28 There have also been numerous theoretical investigations of linear carbon chains that have either been encapsulated in carbon nanotubes or bundled together.29−39 These results led to the notion that bundles of carbon chains © 2016 American Chemical Society

Received: March 24, 2016 Accepted: April 14, 2016 Published: April 14, 2016 1675

DOI: 10.1021/acs.jpclett.6b00671 J. Phys. Chem. Lett. 2016, 7, 1675−1681

Letter

The Journal of Physical Chemistry Letters

Figure 1. Calculated Raman spectra of different sizes of sp carbon chains in the (a) singlet and (b) triplet states.

Figure 2. Comparison of the calculated Raman spectra of (a) different sizes of CnH2 (n = 2, 4, 6, 8, 10, 12) polyynes, (b) C10H2 and C10Au2, (c) C10H2 and C10H2Au20.

therefore only considered the latter. Because these uncapped chains contain two unpaired electrons, the calculations were done for both the singlet and triplet states. After employing a scaling factor of 0.97 to account for systematic overestimation resulting from neglect of anharmonic effects, the vibrational frequencies obtained at the B3LYP level of theory are close to those obtained experimentally.71,72 To examine how higherlevel correlation effects influence the calculated frequencies of open-shell systems, we also evaluated the vibrational frequencies at the much higher CCSD(T) {coupled-cluster with single, double, and perturbative triple substitutions} level of theory on the C6 systems.73 The vibrational frequencies

was used to represent light elements like carbon and hydrogen. 66 Gold was represented with a LANL2DZ pseudopotential basis set.67 We focus on the Raman spectra because they are extremely sensitive to the microstructure of carbon allotropes.68,69 In prior studies of linear carbon chains, the calculated Raman frequencies at the B3LYP/cc-pVDZ level of theory were in good agreement with experiment.47,49 Simulations were carried out on uncapped linear carbon chains with lengths from 6 to 12 atoms. Prior calculations suggested that systems containing only double bonds (C C)n are energetically less stable than those containing alternating single and triple bonds (CC)n.70 We 1676

DOI: 10.1021/acs.jpclett.6b00671 J. Phys. Chem. Lett. 2016, 7, 1675−1681

Letter

The Journal of Physical Chemistry Letters predicted at the two levels of theory are nearly identical (Figure S1, Supporting Information), indicating that the B3LYP level of calculation is adequate for investigating these systems. Short linear carbon chains (Cn, n ≤ 10) have an energetic preference for cumulenic bonding (:CC···CC:) with nearly equivalent bond lengths rather than acetylenic bonding (·CC···CC·) with alternating bond lengths.74,75 Increase in a Cn chain length results in a transition from cumulenic to acetylenic structure,75 predicted to occur at n = 48−52.75 Even though the (CC)n motif was used in the initial structures for the Cn chains, upon optimization they exhibit signs of cumulenic bonding. Raman spectra of experimental samples of cumulenic systems exhibit a broad feature in the 1100−1700 cm −1 region.14,19,76−79 Our calculated spectra (Figure 1) of Cn systems display peaks with significant intensities in this region, consistent with experiment. These peaks are absent if the carbon chains are capped with hydrogen or gold atoms (Figure 2). The major difference in the calculated Raman spectra of several CnH2 (n = 2, 4, 6, 8, 10, 12) polyynes (Figure 2a) and uncapped linear carbon chains (Figure 1) is the absence of peaks with appreciable intensity in the 750−2000 cm−1 region. The conduction electrons of gold nanoparticles collectively oscillate to produce a surface plasmon resonance (SPR) at the most commonly used excitation wavelengths (∼514 nm, ∼ 633 nm) in Raman spectra.80 The SPR leads to a dramatic enhancement in the intensity of some of the Raman-active modes of a molecule when it is adsorbed on gold clusters.80 Such surface-enhanced Raman scattering (SERS) has been widely used, for example, to investigate the Raman spectra of cumulenes and polyynes.19,78,79,81 The pronounced enhancement in the intensity of the Raman spectra mentioned above presumably reflects changes in the structure and properties of the carbon chains. Gold can replace the terminal hydrogens of the C10H2 polyynes to produce goldpolyyne organometallic complexes,82 or a C10H2 polyyne can be adsorbed onto a gold cluster.56,57 Replacement of the terminal hydrogens by gold leads to an increase in the Raman intensity of the CC stretch around 2100 cm−1 (Figure 2b), but there is no perceptible change in the intensities of the peaks in other regions. However, there is a significant increase in the intensity of the peak at ∼1100 cm−1 if a C10H2 polyyne is adsorbed onto a gold Au20 cluster. Similar effects occur in the calculated Raman spectra of a C10H2 polyyne adsorbed on increasing sizes of gold (Aun, n = 4, 8, 12, 20, 32, 42) clusters (Figure S2, Supporting Information), all of which have been observed in experiments.83−86 The calculations also reveal that the multiple peaks around ∼2100 cm−1 coalesce into a single peak as the size of the gold cluster increases. Calculations of C10H2 with smaller gold clusters were also carried out at the BP86/def2TZVP level,64,87,88 wherein the scalar relativistic effects were considered using the zero-order regular approximation (ZORA).89 Both the B3LYP/(cc-pVDZ+LANL2DZ) and BP86/def2-TZVP spectra exhibit similar trends (Figure S3, Supporting Information). The Raman spectra of C10H2 and its complex with the Au12 and Au20 gold clusters were simulated at excitation wavelengths of 514 and 633 nm (Figure 3a,b) using time-dependent density functional calculations.64,90,91 Although there is little change in the intensity of the peaks around ∼2100 cm−1, the peak around ∼1100 cm−1 only appears in the C10H2−Au12 and C10H2−Au20 complexes. The calculated absorption spectra of the C10H2− Au12 complex, together with the computed electron density

Figure 3. Comparison of the calculated Raman spectra of C10H2, C10H2−Au12, and C10H2−Au20, at (a) 514 nm and (b) 633 nm. Compared to C10H2, the peak around ∼1100 cm−1 (blue arrow) becomes prominent in C10H2−Au12 and C10H2−Au20. (c) Calculated absorption spectra of C10H2−Au12 together with the computed electron density difference between the excited and ground states at ∼514 nm and ∼633 nm. Electron flow is from green to red.

difference between the excited and ground states, indicate that a charge transfer from the filled orbitals of the gold cluster to the C10H2 unoccupied orbitals associated with the C−C single bonds occurs in the excited states at ∼514 nm and ∼633 nm (Figure 3c). A similar resonant charge transfer also occurs in pyridine adsorbed on a gold cluster or surface.90,92 Although the same wavelengths were employed by Pan et al. to obtain the Raman spectra of their samples,52 the role of embedded gold clusters were not considered,52 suggesting a misinterpretation by them. A comparison of the calculated absorption spectra,65 of different C10 species (Figure S4, Supporting Information) indicates that the absorption maxima of pure carbon clusters (C10 singlet and C10 triplet) appear at much higher wavelengths than either a C10H2 polyyne or a C10H2−Au12 complex. Furthermore, the calculated absorption maxima, which are in good agreement with earlier experiments of a C 10 H 2 polyyne,28,93 indicate a pronounced red shift in methanol. The calculations indicate that the binding energies of an isolated C10H2 polyyne with different-sized gold clusters range from 0.5 eV (Au4) to 0.1 eV (Au42). These binding energies are similar to that between an isolated acetylene and a Au7 cluster.56 In particular, the interaction of C10H2 with the larger gold clusters is characteristic of a π-type interaction, wherein 1677

DOI: 10.1021/acs.jpclett.6b00671 J. Phys. Chem. Lett. 2016, 7, 1675−1681

Letter

The Journal of Physical Chemistry Letters

Figure 4. Contour plots (isovalue: 0.001 au) of electrostatic potential (ESP) mapped onto the electron density of an isolated C10H2 molecule and its complexes with Au20 (central panel) and Au42 (right panel) clusters. ESP values range from −8.0 × 10−3 a.u. (red) to +1.0 × 10−3 a.u. (blue).

the gold atom interacts with the π electron cloud of the polyyne.56 The gold cluster dramatically alters the electron densities of the CC bonds of C10H2 (Figure 4), showing that it is essential for stabilizing the carbon chains. Self-consistent density functional theory calculations were also carried out to examine the effect of gold on the structure of crystalline bundles of carbon chains.94 Atomic-orbital basis sets were used for the valence electrons and norm-conserving pseudopotentials for the atomic cores.94 The pseudopotentials were constructed using the Trouiller−Martins scheme,95 and the vdW functional was employed for the exchange-correlation potential.96 We employed a double-ζ plus polarization (DZP) basis set. The numerical integrals were performed and projected on a real-space grid with an equivalent cutoff of 200 Ry. Energy minimizations on all the structures were carried out using the conjugate gradient method with a force convergence of 0.004 eV/Å. Several sets of a and c values have been reported for the hexagonal lattice parameters of crystalline carbon chains.8,12,52 The a values range from 4.76 to 9.92 Å and the c values from 7.4 to 15.36 Å.8,12,52 Using the most recently reported parameters (a = 5.78 Å, c = 9.92 Å)52 resulted in a kinked structure (Figure 5a). Calculations in which the lattice parameters were relaxed in the geometry optimizations yielded a linear structure (Figure 5b). A model of kinked carbon chains based on the structure of graphite was proposed to explain the huge variations in the observed lattice parameters.12,52 We suggest that use of accurately determined lattice parameters of a genuine crystalline sample of carbyne, should it exist, would lead to a linear structure. However, we are strongly skeptical since we have yet to see convincing evidence of the occurrence of solid carbyne.9,16,97,98 The effects of gold on the structures of crystals containing carbon chains is both profound and relevant. Kinking appears to increase proportionally to the number of gold atoms adjacent to the chains (Figure 5c,d) and can be attributed to changes in the charge-density profiles of the carbon chains as a result of interactions with the gold atoms (Figure 6). A similar kinking was reported in a recent theoretical study of adsorption of carbon chains on metal surfaces.36 Optimization of the structure of a crystal with carbon chains intercalated with Au8 clusters led to the formation of chains of gold along the carbon chains (Figure S5, Supporting Information). An earlier calculation of carbon nanotubes intercalated with small gold clusters resulted in a similar formation of a gold chain alongside the nanotube axis.59 Implications of our results include: (a) uncapped finite-sized carbon chains exhibit a cumulenic character, indicated by the appearance of a peak in the 1100−1700 cm−1 region of the Raman spectra; (b) although spectra of isolated hydrogencapped carbon chains contain a strong peak around ∼2100 cm−1, no peaks appear in the 750−2100 cm−1 region; (c)

Figure 5. Optimized 2 × 2 × 2 structures of hypothetical crystals of carbon chains based on an (a) fixed lattice, (b) variable lattice, and (c,d) similar carbon chains with increasing concentrations of gold atoms between the carbon chains. Unit cells of these pseudocarbyne crystals are indicated by light black lines.

interaction of polyynes with small gold clusters leads to a peak around ∼1100 cm−1, which arises from a C−C stretch; (d) the electron-density profiles indicate that gold atoms have a significant interaction with the polyyne, which we suggest leads to the electronic stabilization of polyynes in crystals; and (e) intercalated gold atoms or clusters in crystals consisting largely of carbon chains leads to kinking that results from interaction of the electronic densities of the carbon chains and gold. On the basis of the above, the evidence suggests that the material synthesized and characterized by Pan et al.52 is more consistent with C/Au clusters than with pure or even hydrogen-capped carbon chains. Crystals such as in Figure 5 are plausible explanations for their experimental results. The current work opens the possibility of a range of chargestabilized carbon chains having different stoichiometries, with gold and also other noble metals, but whose occurrence 1678

DOI: 10.1021/acs.jpclett.6b00671 J. Phys. Chem. Lett. 2016, 7, 1675−1681

Letter

The Journal of Physical Chemistry Letters

Figure 6. Profiles of (a) the electron density (ρ(r), isovalue: 0.08 e/Å3) and (b) the charge density difference (Δρ(r), isovalue: ± 0.004 e/Å3) of the gold−carbon chain complex in Figure 5d. The red regions around the carbon atoms show charge enhancement, and the green regions show depletion. They are, respectively, contracted and distended near the carbon atoms (black arrows) close to gold. (10) Whittaker, A. G. Carbyne Forms of Carbon: Evidence for Their Existence. Science 1985, 229, 485−486. (11) Smith, P. P. K.; Buseck, P. R. Reply to Carbyne Forms of Carbon: Do They Exist? Science 1985, 229, 486−487. (12) Heimann, R. B.; Kleiman, J.; Salansky, N. M. A Unified Structural Approach to Linear Carbon Polytypes. Nature 1983, 306, 164−167. (13) Lagow, R. J.; Kampa, J. J.; Wei, H.-C.; Battle, S. L.; Genge, J. W.; Laude, D. A.; Harper, C. J.; Bau, R.; Stevens, R. C.; Haw, J. F.; Munson, E. Synthesis of Linear Acetylenic Carbon: The sp Carbon Allotrope. Science 1995, 267, 362−367. (14) Kastner, J.; Kuzmany, H.; Kavan, L.; Dousek, F. P.; Kürti, J. Reductive Preparation of Carbyne with High Yield. An in Situ Raman Scattering Study. Macromolecules 1995, 28, 344−353. (15) Baughman, R. H. Dangerously Seeking Linear Carbon. Science 2006, 312, 1009−1010. (16) Kroto, H. Carbyne and Other Myths about Carbon. Chem. World 2010, 7, 37. (17) Chalifoux, W. A.; Tykwinski, R. R. Synthesis of Polyynes to Model the sp-Carbon Allotrope Carbyne. Nat. Chem. 2010, 2, 967− 971. (18) Jevric, M.; Nielsen, M. B. Synthetic Strategies for Oligoynes. Asian J. Org. Chem. 2015, 4, 286−295. (19) Casari, C. S.; Tommasini, M.; Tykwinski, R. R.; Milani, A. Carbon-atom Wires: 1-D Systems with Tunable Properties. Nanoscale 2016, 8, 4414−4435. (20) Jensen, W. B. The Origin of the Term Allotrope. J. Chem. Educ. 2006, 83, 838−839. (21) IUPAC Nomenclature of Inorganic Chemistry; Connelly, N. G.; Damhus, T., Eds.; RSC Publishing: Cambridge, U.K., 2005; p 49. (22) Johnson, T. R.; Walton, D. R. M. Silylation as a Protective Method in Acetylene Chemistry. Tetrahedron 1972, 28, 5221−5236. (23) Eastmond, R.; Johnson, T. R.; Walton, D. R. M. Base-Catalyzed Cleavage of Silyl-Substituted Polyynes. Attenuation of Hydrocarbon Acidity and Transmission of Substituent Electrical Effects in LongChain Conjugated Polyacetylenes. J. Organomet. Chem. 1973, 50, 87− 92. (24) Rubin, Y.; Lin, S. S.; Knobler, C. B.; Anthony, J.; Boldi, A. M.; Diederich, F. Solution-Spray Flash Vacuum Pyrolysis: A New Method for the Synthesis of Linear Poliynes with Odd Numbers of C-C Bonds from Substituted 3,4-Dialkynyl-3-cyclobutene1-,2-diones. J. Am. Chem. Soc. 1991, 113, 6943−6949. (25) Gibtner, T.; Hampel, F.; Gisselbrecht, J.-P.; Hirsch, A. End-Cap Stabilized Oligoynes: Model Compounds for the Linear sp Carbon Allotrope Carbyne. Chem. - Eur. J. 2002, 8, 408−432. (26) Hlavatý, J.; Kavan, L.; Kubišta, J. Carbonaceous Materials from End-capped Alkynes. Carbon 2002, 40, 345−349. (27) Eisler, S.; Slepkov, A. D.; Elliott, E.; Luu, T.; McDonald, R.; Hegmann, F. A.; Tykwinski, R. R. Polyynes as a Model for Carbyne: Synthesis, Physical Properties, and Nonlinear Optical Response. J. Am. Chem. Soc. 2005, 127, 2666−2676.

remains to be tested experimentally. Such pseudocarbynes could prove to have interesting properties. However, the bottom line is that carbyne, the linear carbon allotrope, remains as elusive and theoretical as ever.



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpclett.6b00671. Calculated Raman spectra of several C10H2−Aun gold cluster complexes at different levels of theory, absorption spectra of several C10 species, and the optimized structure of Au8 gold cluster intercalated carbon chains. (PDF)



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Notes

The authors declare no competing financial interest.

■ ■

ACKNOWLEDGMENTS We thank Dr. Péter Németh for helpful discussions. REFERENCES

(1) Kroto, H. W.; Heath, J. R.; O’Brien, S. C.; Curl, R. F.; Smalley, R. E. C60: Buckminsterfullerene. Nature 1985, 318, 162−163. (2) Iijima, S. Helical Microtubules of Graphitic Carbon. Nature 1991, 354, 56−58. (3) Novoselov, K. S.; Geim, A. K.; Morozov, S. V.; Jiang, D.; Zhang, Y.; Dubonos, S. V.; Grigorieva, I. V.; Firsov, A. A. Electric Field Effect in Atomically Thin Carbon Films. Science 2004, 306, 666−669. (4) Novoselov, K. S.; Jiang, D.; Schedin, F.; Booth, T. J.; Khotkevich, V. V.; Morozov, S. V.; Geim, A. K. Two-Dimensional Atomic Crystals. Proc. Natl. Acad. Sci. U. S. A. 2005, 102, 10451−10453. (5) Cantrill, S. Carbon-Based Curiosities. Nat. Chem. 2009, 1, 170− 171. (6) Banhart, F. Chains of Carbon Atoms: A Vision or a New Nanomaterial? Beilstein J. Nanotechnol. 2015, 6, 559−569. (7) Korshak, V. V.; Kasatochkin, V. I.; Kudriavt, I. P.; Usenbaev, K.; Sladov, A. M. Synthesis and Properties of Polyacetylene. Dokl. Akad. Nauk SSSR 1961, 136, 1342. (8) Kasatochkin, V. I.; Mel’nichenko, V. M.; Elizen, V. M. Electron Diffraction by Single Crystals of Carbyne. Polym. Sci. U.S.S.R. 1975, 17, 2167−2173. (9) Smith, P. P. K.; Buseck, P. R. Carbyne Forms of Carbon: Do They Exist? Science 1982, 216, 984−986. 1679

DOI: 10.1021/acs.jpclett.6b00671 J. Phys. Chem. Lett. 2016, 7, 1675−1681

Letter

The Journal of Physical Chemistry Letters (28) Cataldo, F. Polyynes and Cyanopolyynes Synthesis from the Submerged Electric Arc: About the Role Played by the Electrodes and Solvents in Polyynes Formation. Tetrahedron 2004, 60, 4265−4274. (29) Fan, X.; Liu, L.; Lin, J.; Shen, Z.; Kuo, J.-L. Density Functional Theory Study of Finite Carbon Chains. ACS Nano 2009, 3, 3788− 3794. (30) Luo, W.; Windl, W. First Principles Study of the Structure and Stability of Carbynes. Carbon 2009, 47, 367−383. (31) Hu, Y. H. Stability of sp Carbon (Carbyne) Chains. Phys. Lett. A 2009, 373, 3554−3557. (32) Cahangirov, S.; Topsakal, M.; Ciraci, S. Long-Range Interactions in Carbon Atomic Chains. Phys. Rev. B: Condens. Matter Mater. Phys. 2010, 82, 195444. (33) Lin, Z. Z.; Ning, X. J. Controlling the Electronic Properties of Monatomic Carbon Chains. Europhys. Lett. 2011, 95, 47012. (34) Hu, Y. H. Bending Effect of sp-Hybridized Carbon (Carbyne) Chains on Their Structures and Properties. J. Phys. Chem. C 2011, 115, 1843−1850. (35) Liu, M.; Artyukhov, V. I.; Lee, H.; Xu, F.; Yakobson, B. I. Carbyne from First Principles: Chain of C Atoms, a Nanorod or a Nanorope. ACS Nano 2013, 7, 10075−10082. (36) Yuan, Q.; Ding, F. Formation of Carbyne and Graphyne on Transition Metal Surfaces. Nanoscale 2014, 6, 12727−12731. (37) Shi, L.; Rohringer, P.; Suenaga, K.; Niimi, Y.; Kotakoski, J.; Meyer, J. C.; Peterlik, H.; Ayala, P.; Pichler, T. Confined Linear Carbon Chains: A Route to Bulk Carbyne. 2015, arXiv:1507.04896v2. arXiv.org e-Print archive. http://arxiv.org/abs/1507.04896v2 (accessed December 2015). (38) Timoshevskii, A.; Kotrechko, S.; Matviychuk, Y. Atomic Structure and Mechanical Properties of Carbyne. Phys. Rev. B: Condens. Matter Mater. Phys. 2015, 91, 245434. (39) Cannella, C. B.; Goldman, N. Carbyne Fiber Synthesis within Evaporating Metallic Liquid Carbon. J. Phys. Chem. C 2015, 119, 21605−21611. (40) Davenport, M. Contention Over Carbyne. Chem. Eng. News 2015, 93, 30−31. (41) Ravagnan, L.; Siviero, F.; Lenardi, C.; Piseri, P.; Barborini, E.; Milani, P.; Casari, C. S.; Bassi, A. L.; Bottani, C. E. Cluster-Beam Deposition and in situ Characterization of Carbyne-Rich Carbon Films. Phys. Rev. Lett. 2002, 89, 285506. (42) Hu, A.; Griesing, S.; Rybachuk, M.; Lu, Q.-B.; Duley, W. W. Nanobuckling and X-ray Photoelectron Spectra of Carbyne-Rich Tetrahedral Carbon Films Deposited by Femtosecond Laser Ablation at Cryogenic Temperatures. J. Appl. Phys. 2007, 102, 074311. (43) Szczepanski, J.; Fuller, J.; Ekern, S.; Vala, M. Electronic Absorption and Resonance Raman Spectra of Large Linear Carbon Clusters Isolated in Solid Argon. Spectrochim. Acta, Part A 2001, 57, 775−786. (44) Jin, C.; Lan, H.; Peng, L.; Suenaga, K.; Iijima, S. Deriving Carbon Atomic Chains from Graphene. Phys. Rev. Lett. 2009, 102, 205501. (45) Cretu, O.; Botello-Mendez, A. R.; Janowska, I.; Pham-Huu, C.; Charlier, J.-C.; Banhart, F. Electrical Transport Measured in Atomic Carbon Chains. Nano Lett. 2013, 13, 3487−3493. (46) Zhao, X.; Ando, Y.; Liu, Y.; Jinno, M.; Suzuki, T. Carbon Nanowire Made of a Long Linear Carbon Chain Inserted Inside a Multiwalled Carbon Nanotube. Phys. Rev. Lett. 2003, 90, 187401. (47) Andrade, N. F.; Vasconcelos, T. L.; Gouvea, C. P.; Archanjo, B. S.; Achete, C. A.; Kim, Y. A.; Endo, M.; Fantini, C.; Dresselhaus, M. S.; Filho, A. G. S. Linear Carbon Chains Encapsulated in Multiwall Carbon Nanotubes: Resonance Raman spectroscopy and Transmission Electron Microscopy Studies. Carbon 2015, 90, 172−180. (48) Ott, A. K.; Rechtsteiner, G. A.; Felix, C.; Hampe, O.; Jarrold, M. F.; Van Duyne, R. P.; Raghavachari, K. Raman Spectra and Calculated Vibrational Frequencies of size-selected C16, C18, and C20 clusters. J. Chem. Phys. 1998, 109, 9652−9655. (49) Tabata, H.; Fujii, M.; Hayashi, S.; Doi, T.; Wakabayashi, T. Raman and Surface-Enhanced Raman Scattering of a Series of SizeSeparated Polyynes. Carbon 2006, 44, 3168−3176.

(50) Lucotti, A.; Tommasini, M.; Zoppo, M. D.; Castiglioni, C.; Zerbi, G.; Cataldo, F.; Casari, C. S.; Bassi, A. L.; Russo, V.; Bogana, M.; Bottani, C. E. Raman and SERS Investigation of Isolated sp Carbon Chains. Chem. Phys. Lett. 2006, 417, 78−82. (51) Xue, K.-H.; Chen, S.-P.; Wang, L.-X.; Wei, R.-B.; Xu, S.-M.; Cui, L.; Mao, B.-W.; Tian, Z.-Q.; Zen, C.-H.; Sun, S.-G.; Wu, L.-J.; Zhu, Y.M. Fabrication and Characterization of Oriented Carbon Atom Wires Assembled on Gold. Chem. Phys. Lett. 2009, 469, 284−288. (52) Pan, B.; Xiao, J.; Li, J.; Liu, P.; Wang, C.; Yang, G. Carbyne with Finite Length: The One-dimensional sp Carbon. Sci. Adv. 2015, 1, e1500857. (53) Herzing, A. A.; Kiely, C. J.; Carley, A. F.; Landon, P.; Hutchings, G. J. Identification of Active Gold Nanoclusters on Iron Oxide Supports for CO Oxidation. Science 2008, 321, 1331−1335. (54) La Torre, A.; Giménez-López, M. C.; Fay, M. W.; Rance, G. A.; Solomonsz, W. A.; Chamberlain, T. W.; Brown, P. D.; Khlobystov, A. N. Assembly, Growth, and Catalytic Activity of Gold Nanoparticles in Hollow Carbon Nanofibers. ACS Nano 2012, 6, 2000−2007. (55) Liu, X.; Cui, X.; Liu, Y.; Yin, Y. Stabilization of Ultrafine Metal Nanocatalysts on Thin Carbon Sheets. Nanoscale 2015, 7, 18320− 18326. (56) Ferraro, F.; Pérez-Torres, J. F.; Hadad, C. Z. Selective Catalytic Activation of Acetylene by a Neutral Gold Cluster of Experimentally Known Gas-Phase Geometry. J. Phys. Chem. C 2015, 119, 7755−7764. (57) Tsipis, A. C. Unveiling the Nature of Binding Interactions of Acetylene and Ethylene with Triangular Coinage Metal Clusters: A DFT Computational Study. Organometallics 2010, 29, 354−363. (58) Lyalin, A.; Taketsugu, T. Adsorption of Ethylene on Neutral, Anionic, and Cationic Gold Clusters. J. Phys. Chem. C 2010, 114, 2484−2493. (59) Deepak, J.; Pradeep, T.; Waghmare, U. V. Interaction of Small Gold Clusters with Carbon Nanotube Bundles: Formation of Gold Atomic Chains. J. Phys.: Condens. Matter 2010, 22, 125301. (60) Amft, M.; Sanyal, B.; Eriksson, O.; Skorodumova, N. V. Small Gold Clusters on Graphene, Their Mobility and Clustering: A DFT Study. J. Phys.: Condens. Matter 2011, 23, 205301. (61) Pulido, A.; Boronat, M.; Corma, A. Theoretical Investigation of Gold Clusters Supported on Graphene Sheets. New J. Chem. 2011, 35, 2153−2161. (62) Becke, A. D. Density-Functional Thermochemistry. III. The Role of Exact Exchange. J. Chem. Phys. 1993, 98, 5648−5652. (63) Lee, C.; Yang, W.; Parr, R. G. Development of the Colle-Salvetti Correlation-Energy Formula into a Functional of the Electron Density. Phys. Rev. B: Condens. Matter Mater. Phys. 1988, 37, 785−789. (64) Neese, F. The Orca Program System. WIREs Comput. Mol. Sci. 2012, 2, 73−78. (65) 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, revision A.1; Gaussian, Inc.: Wallingford, CT, 2009. (66) Dunning, T. H. Gaussian Basis Sets for use in Correlated Molecular Calculations. I. The Atoms Boron Through Neon and Hydrogen. J. Chem. Phys. 1989, 90, 1007−1023. (67) Hay, P. J.; Wadt, W. R. Ab initio Effective Core Potentials for Molecular Calculations. Potentials for K to Au Including the Outermost Core Orbitals. J. Chem. Phys. 1985, 82, 299−310. (68) Saito, R.; Hofmann, M.; Dresselhaus, G.; Jorio, A.; Dresselhaus, M. S. Raman Spectroscopy of Graphene and Carbon Nanotubes. Adv. Phys. 2011, 60, 413−550. (69) Wang, Y.; Alsmeyer, D. C.; McCreery, R. L. Raman Spectroscopy of Carbon Materials: Structural Basis of Observed Spectra. Chem. Mater. 1990, 2, 557−563. (70) Rice, M. J.; Phillpot, S. R.; Bishop, A. R.; Campbell, D. K. Solitons, Polarons, and Phonons in the Infinite Polyyne Chain. Phys. Rev. B: Condens. Matter Mater. Phys. 1986, 34, 4139−4149. (71) Villa, M.; Senent, M. L.; Dominguez-Gomez, R.; Alvarez-Bajo, O.; Carvajal, M. CCSD(T) Study of Dimethyl-Ether Infrared and Raman Spectra. J. Phys. Chem. A 2011, 115, 13573−13580. 1680

DOI: 10.1021/acs.jpclett.6b00671 J. Phys. Chem. Lett. 2016, 7, 1675−1681

Letter

The Journal of Physical Chemistry Letters (72) Merrick, J. P.; Moran, D.; Radom, L. An Evaluation of Harmonic Vibrational Frequency Scale Factors. J. Phys. Chem. A 2007, 111, 11683−11700. (73) Pople, J. A.; Head-Gordon, M.; Raghavachari, K. Quadratic Configuration Interaction - A General Technique for Determining Electron Correlation Energies. J. Chem. Phys. 1987, 87, 5968−5975. (74) Van Orden, A.; Saykally, R. J. Small Carbon Clusters: Spectroscopy, Structure and Energetics. Chem. Rev. 1998, 98, 2313− 2358. (75) Yang, S.; Kertesz, M. Linear Cn Clusters: Are They Acetylenic or Cumulenic? J. Phys. Chem. A 2008, 112, 146−151. (76) Milani, A.; Tommasini, M.; Del Zoppo, M.; Castiglioni, C.; Zerbi, G. Carbon Nanowires: Phonon and π-electron Confinement. Phys. Rev. B: Condens. Matter Mater. Phys. 2006, 74, 153418. (77) Tommasini, M.; Milani, A.; Fazzi, D.; Lucotti, A.; Castiglioni, C.; Januszewski, J. A.; Wendinger, D.; Tykwinski, R. R. π-Conjugation and End Group Effects in Long Cumulenes: Raman Spectroscopy and DFT Calculations. J. Phys. Chem. C 2014, 118, 26415−26425. (78) Milani, A.; Lucotti, A.; Russo, V.; Tommasini, M.; Cataldo, F.; Bassi, A. L.; Casari, C. S. Charge Transfer and Vibrational Structure of sp-Hybridized Carbon Atomic Wires Probed by Surface Enhanced Raman Spectroscopy. J. Phys. Chem. C 2011, 115, 12836−12843. (79) Milani, A.; Tommasini, M.; Russo, V.; Bassi, A. L.; Lucotti, A.; Cataldo, F.; Casari, C. S. Raman Spectroscopy as a Tool to Investigate the Structure and Electronic Properties of Carbon-Atom Wires. Beilstein J. Nanotechnol. 2015, 6, 480−491. (80) Schlücker, S. Surface-Enhanced Raman Spectroscopy: Concepts and Chemical Applications. Angew. Chem., Int. Ed. 2014, 53, 4756− 4795. (81) Compagnini, G.; Patanè, G.; D’Urso, L.; Puglisi, O.; Cataliotti, R. S.; Pignataro, B. On the Interaction of Carbon Nanowires with Noble Metals through a Study of Their Surface-Enhanced Raman Spectra. J. Phys. Chem. C 2008, 112, 20301−20306. (82) Liu, H.-T.; Xiong, X.-G.; Dau, P. D.; Wang, Y.-L.; Huang, D.-L.; Li, J.; Wang, L.-S. Probing the Nature of Gold−Carbon Bonding in Gold−Alkynyl complexes. Nat. Commun. 2013, 4, 2223. (83) Li, J.; Li, X.; Zhai, H.-J.; Wang, L.-S. Au20: A Tetrahedral Cluster. Science 2003, 299, 864−867. (84) Gao, Y.; Zeng, X. C. Au42: An Alternative Icosahedral Golden Fullerene Cage. J. Am. Chem. Soc. 2005, 127, 3698−3699. (85) Karttunen, A. J.; Linnolahti, M.; Pakkanen, T. A.; Pyykkö, P. Icosahedral Au72: A Predicted Chiral and Spherically Aromatic Golden Fullerene. Chem. Commun. 2008, 465−467. (86) Wang, Z. W.; Palmer, R. E. Direct Atomic Imaging and Dynamical Fluctuations of the Tetrahedral Au20 Cluster. Nanoscale 2012, 4, 4947−4949. (87) Weigend, F.; Ahlrichs, R. Balanced Basis Sets of Split Valence, Triple Zeta Valence and Quadruple Zeta Valence Quality for H to Rn: Design and Assessment of Accuracy. Phys. Chem. Chem. Phys. 2005, 7, 3297−3305. (88) Weigend, F. Accurate Coulomb-Fitting Basis Sets for H to Rn. Phys. Chem. Chem. Phys. 2006, 8, 1057−65. (89) van Wüllen, C. Molecular Density Functional Calculations in the Regular Relativistic Approximation. Method, Application to Coinage Metal Diatomics, Hydrides, Fluorides and Chlorides, and Comparison with First-Order Relativistic Calculations. J. Chem. Phys. 1998, 109, 392−399. (90) Aikens, C. M.; Schatz, G. C. TDDFT Studies of Absorption and SERS Spectra of Pyridine Interacting with Au20. J. Phys. Chem. A 2006, 110, 13317−13324. (91) Mullin, J. M.; Autschbach, J.; Schatz, G. C. Time-dependent Density Functional Methods for Surface Enhanced Raman Scattering (SERS) Studies. Comput. Theor. Chem. 2012, 987, 32−41. (92) Katayama, K.; Shibamoto, K.; Sawada, T. Direct Observation of Ultrafast Charge Transfer in Relation to the Surface Enhanced Raman Scattering Activation Detected by Transient Refecting Grating Spectroscopy. Chem. Phys. Lett. 2001, 345, 265−271.

(93) Matsutani, R.; Inoue, K.; Wada, N.; Kojima, K. Wavelength Dependence of Polyyne Preparation by Liquid-phase Laser Ablation Using Pellet Targets. Chem. Commun. 2011, 47, 5840−5842. (94) Soler, J. M.; Artacho, E. M.; Gale, J. D.; Garcia, A.; Junquera, J.; Ordejón, P.; Sánchez-Portal, D. The SIESTA Method for ab Initio Order-N Materials Simulation. J. Phys.: Condens. Matter 2002, 14, 2745−2779. (95) Troullier, N.; Martins, J. L. Efficient Pseudopotentials for PlaneWave Calculations. Phys. Rev. B: Condens. Matter Mater. Phys. 1991, 43, 1993−2006. (96) Román-Pérez, G.; Soler, J. M. Efficient Implementation of a van der Waals Density Functional: Application to Double-Wall Carbon Nanotubes. Phys. Rev. Lett. 2009, 103, 096102. (97) Buseck, P. R.; Tarakeshwar, P.; Németh, P.; Kroto, H. W. Carbyne Remains Elusive. Sci. Adv. 2016, http://advances.sciencemag. org/content/1/9/e1500857.e-letters. (98) Attempts to secure original material from the authors of ref 52 for independent checking were unsuccessful, but we plan similar syntheses to test the models. (99) Shi, L.; Rohringer, P.; Suenaga, K.; Niimi, Y.; Kotakoski, J.; Meyer, J. C.; Peterlik, H.; Wanko, M.; Cahangirov, S.; et al. Confined Linear Carbon Chains as a Route to Bulk Carbyne. Nat. Mater. 2016, DOI: 10.1038/nmat4617.



NOTE ADDED IN PROOF A recent study reports the production of exceptionally long and stable carbon chains encapsulated in carbon nanotubes.99 Based on the Raman spectra, it is evident that a nanotube is essential for stabilizing the encapsulated carbon chain. Therefore, such an encapsulated chain provides another example of a pseudocarbyne.

1681

DOI: 10.1021/acs.jpclett.6b00671 J. Phys. Chem. Lett. 2016, 7, 1675−1681