Measurement of a linear free energy relationship one molecule at a time B. V. Rao, K.-Y. Kwon, A. Liu, and L. Bartels* Department of Chemistry, University of California, Riverside, CA 92521 Edited by Gerhard Ertl, Max Planck Society for the Advancement of Science, Berlin, Germany, and approved October 15, 2004 (received for review August 24, 2004)
A systematic study of the dehydrogenation of substituted thiophenols by controlled charge injection from the tip of a scanning tunneling microscope (STM) reveals a pronounced dependence of the reaction yield on the position and the chemical nature of the substituent. We evaluate the dehydrogenation rate of para-halosubstituted species within a linear free energy relationship, namely the Hammett equation. The resultant value of 1.4 can faithfully predict the reaction rates of molecules that are meta-halo-substituted or para-methyl-substituted. The positive sign of suggests a negatively charged transition state at the core of the STM-induced process, and the magnitude of the value indicates that the presence of the substrate does not preclude substantial substituent effects. The applicability of the Hammett equation to singlemolecule chemistry offers facile prediction of the rate of STMbased single-molecule chemistry in a field, which so far has been addressed by focusing on involved quantum-mechanical modeling of its underlying processes. Hammett equation 兩 scanning tunneling microscopy 兩 single-molecule chemistry
T
he groundbreaking work of the groups of Eigler (1), Avouris (2), Gimzewski (3), Weiss (4), Rieder (5), and, most recently, Ho (6), has converted the scanning tunneling microscope (STM) from an imaging instrument (7) to a versatile tool capable (i) of control of chemical reactions (8) and (ii) of spectroscopic identification of their reactants and products (9). With the availability of commercial instrumentation for these procedures (10), the need for facile prediction of STM-based singlemolecule chemistry becomes imminent. In conventional organic chemistry, linear free energy relationships such as the Hammett equation† have proven to be powerful tools for the prediction of the reactivity of a broad range of compounds (13–17). Here, we present the transfer of the Hammett equation to the realm where individual molecules are excited one at a time to form specific chemical bonds. The Hammett equation (Eq. 1) predicts the impact of the substituent X on the rate (or equilibrium) constant, k, of the reaction Z 3 Z⬘, where Z is a reactive group attached to the same aromatic ring as X.
Hammett equation: log 共k x兾k 0兲 ⫽ ⫻
[1]
The value describes how a specific substituent affects the acidity of benzoic acid in water (i.e., Z ⫽ COOH, Z⬘ ⫽ COO⫺). These values can be found tabulated in the literature; for the purpose of this manuscript, they are assumed to be universal constants. The value indicates how susceptible a reaction Z 3 Z⬘ is to substitution at the aromatic ring. It compares Z’s susceptibility to the susceptibility of benzoic acid in water, which has by definition a value of 1. Because of its simplicity and strong predictive power, the Hammett equation found applications all over organic chemistry, and their number is still growing. 17920 –17923 兩 PNAS 兩 December 28, 2004 兩 vol. 101 兩 no. 52
We explore the validity of a linear free energy relationship similar to the Hammett equation at the realm of singlemolecule chemistry by a systematic study of thiophenols (TP) rather than benzoic acids, i.e., Z ⫽ SH. TPs are relatively simple molecules that, nevertheless, exhibit a wide range of physical and chemical properties. As an important building block of most molectronic concepts proposed to date, their STM-induced reactivity has potential for future STM-based fine-tuning of molectronic circuitry. In a previous study (18), we showed that the thiol group of 2,5-dichlorothiophenol does not dehydrogenate upon deposition onto Cu(111) at 16 K. Attachment of electrons at a bias of several hundred millivolts of either polarity can cause hydrogen abstraction, and the resultant thiolate is locked rigidly into the surface. We ascertained these findings by STM-based vibrational spectroscopy (18). Before dehydrogenation, TPs rotate around their S substrate anchor even at 16 K, suggesting that the interaction of the aryl moiety with the substrate is very weak. We attribute the high rational mobility of the aryl moiety to the length of the on-top S–Cu bond (19), which suspends the aryl moiety substantially above the substrate. Their weak substrate interactions render TPs highly suitable for the search for a substrate-independent phenomenon. We prepare our Cu(111) single-crystal sample by several cycles of sputtering (Ar⫹, 2 kV) and annealing (700 K) at a base pressure ⬍3 ⫻ 10⫺11 torr (1 torr ⫽ 133 Pa). Low coverages (⬍0.1% of a monolayer) of TPs are deposited after cooling to 16 K. Fig. 1a shows an STM image of TP molecules on Cu(111) at 16 K. The f lower shape of the adsorbates originates from the benzene ring rotating around the thiol bond through six sets of orientations where it is in registry with the top substrate layer (20). Substituted TPs behave in a similar fashion (Fig. 1 b–d). The increase of the feature diameter from para-F-TP (Fig. 1b) to para-Br-TP (Fig. 1d) ref lects the increased size of the halo-substituent (Supporting Text and Figs. 5–7, which are published as supporting information on the PNAS web site). We investigate the reactivity of the surface thiols by exposing their rotation center to a constant tunneling current at ⫹700 mV sample bias while arresting the feedback loop. After ⬇1 ms to 1 s, an abrupt step in the tunneling current appears (Fig. 2a) that corresponds to the formation of the thiolate. The thiolate is locked rigidly into the surface and appears as a single oval protrusion similar to that in ref. 18. We repeated this experiment on hundreds of molecules at identical bias and tunneling current. The resultant distribution of durations follows a decaying exponential (Fig. 2a Inset), from which a decay rate can be calculated at typically ⱕ15% fit uncertainty. By repeating the experiment at different currents [but identical voltage, as realized by slight (⬍10%) variation of the tip This paper was submitted directly (Track II) to the PNAS office. Abbreviations: STM, scanning tunneling microscope; TP, thiophenol. *To whom correspondence should be addressed. E-mail:
[email protected]. †The
Hammett equation (11) is part of any introductory course in organic chemistry (12).
© 2004 by The National Academy of Sciences of the USA
www.pnas.org兾cgi兾doi兾10.1073兾pnas.0406223101
height] the current dependence of the dehydrogenation rate can be established (5, 6). Fig. 2b shows a plot of the dehydrogenation rate versus the applied current for TP, para-F-TP, para-Cl-TP, and para-Br-TP.‡ The linear relationship between rate and current for each of the molecules indicates that a one-electron process is at the core of the thiol activation (5, 6). The resultant rate constants, kpara-X, are indicated in units of molecules per pC. For the sake of testing the predictive power of a linear free energy relationship similar to the Hammett equation for STM-based thiol activation, we use this data set to calculate a corresponding value by fitting our measured activation rates to values from the literature (13). A value of 1.4 ⫾ 0.1 yields the best correlation between observed and literature values (central section of Fig. 3a). The power of this modified Hammett correlation, if it is applicable, lies in its ability to predict from this value the activities of TPs that carry (i) different substituents or (ii) the same substituent at a different (i.e., meta) position of the aromatic ring. We test the applicability of this modified Hammett correlation by repeating the described experiments using para-methyl-substitution and meta-halo-substitution. Fig. 4 a–c shows STM-images of para-methyl-TP, meta-F-TP, and meta-Cl-TP, respectively. Rotation of the substituted aryl moiety is clearly not limited to para-halo-substitution. meta-halo-substitution leads to an asymmetric shape of the reactants that originates from the position of the halogen off the S-aryl axis. The asymmetric nature of the molecule causes the bent shape of features in Fig. 4 b and c. There are two enantiomers of each of the meta-substituted molecules found on the surface. Excitation of these species under identical conditions reveals a linear dependence of the activation rate on the tunneling ‡The
data set obtained in this study comprises several thousand molecular excitations performed during ⬇4 months of measurements with ⬎12 sample preparations separated each by bake-out of the home-built STM system, extensive sample cleaning, and tip reforming. The behavior described has been found in all instances and does not reflect observations limited to a small subset of experiments only.
Rao et al.
Fig. 2. STM-based excitation of unsubstituted and para-halosubstituted thiophenol molecules. (a) Tunneling current, I, during the excitation of a TP molecule at a temperature of 16 K and a bias of 700 mV. The sharp step indicates the dehydrogenation event. Inset shows as an example the distribution of durations, t, until the dehydrogenation occurred for all experiments conducted on TP at ⬇120 pA. Typically, the fit uncertainty is ⱕ⫾15% (here ⫾13%). (b) Dehydrogenation rate, k, of (substituted) TPs versus current at a 700-mV sample bias. The linear dependence (all R2 ⬎ 0.99) indicates that a one-electron process causes the dehydrogenation.
current (Fig. 4d), i.e., similar to the previous set of molecules (Fig. 2). The rate constants for meta-substitution are found to exceed those for para-substitution by a factor of ⬇2. In contrast to halo-substitution, methyl-substitution reduces the activity of the thiol group substantially. The rate constants of this second set of measurements span almost an order of magnitude in value, and they exceed the regime probed earlier both toward smaller and larger rate constants. The shaded portion of Fig. 3a compares literature values (13) for para-methyl-TP and methyl-X-TPs with those derived from the second set of experiments assuming the value of 1.4 that has been obtained from para-halo-substitution. The correlation is very good for para-methyl-substitution, and our application of the Hammett equation overestimates the effect of meta-halo-substitution only by an amount that is comparable to the spread of literature values (compare refs. 13 and 20). Commonly, the value of a new reaction is derived by plotting the logarithm of the ratio of the measured rate or equilibrium constants versus the values obtained on benzoic acids in water (Fig. 3b). When using the entire data set, the best fit has a slope of ⫽ 1.4 at a fit uncertainty of ⫾0.1, i.e., similar to the previously inferred one. We also attempted measurement of the activation rate of meta-Br-TP, which corresponds to the largest value in our set of substances. Unfortunately, meta-Br-TP appears to deprotonate before we can reproducibly identify any rotating species on the surface. Although this confirms the prediction of the Hammett equation qualitatively, no quantitative value could be derived. Conventionally, the value of a reaction is used to infer the reaction mechanism. Positive values indicate that electronwithdrawing groups increase the molecule’s reactivity, which is commonly ascribed to the stabilization of a negatively charged transition state either by dipolar or resonance effects (13). Our PNAS 兩 December 28, 2004 兩 vol. 101 兩 no. 52 兩 17921
CHEMISTRY
Fig. 1. STM images of unsubstituted and para-halosubstituted thiophenol molecules on Cu(111). (a) Three-dimensional representation of TP on Cu(111) at 16 K (80 ⫻ 50 Å, ⫺300 mV, 56 pA). TP has a height of 1.4 Å. The crystallographic directions of the substrate are indicated and the Inset shows the structure of a TP molecule. (b–d) STM images of para-F-TP (⫺250 mV, 42 pA), para-Cl-TP (⫺270 mV, 40 pA), and para-Br-TP (⫺200 mV, 59 pA) taken at 16 K. The flower shape originates from the molecule’s rotation around the S–Cu bond as described in ref. 18. The increase in diameter is caused by the size of the substituents.
Fig. 3. Comparison of our data with Hammett values. (a) Comparison between measured (open bars) and literature (filled bars) Hammett values assuming a value of ⬇1.4, derived from the best fit for the three para-X-TP molecules (unshaded area). The satisfactory correlation between measured and literature values in the shaded areas shows that a value of 1.4 is capable of predicting reaction rates of substituted TPs. (b) Conventional representation of the correlation between the rate constants. The x axis indicates the ratio of the acid dissociation constants of benzoic acids in water (i.e., the value), and the y axis shows the ratio measured for TPs on Cu(111). A linear fit yields a of 1.4 at a fit uncertainty of ⫾0.1.
findings suggest that in the STM-based reaction, a negatively charged transition state plays an important role (which does not conf lict with the electron attachment process used here). The values of phenol in water and TP in 1:1 water/ethanol solution are 2.5 and 2.2, respectively (14) [values were derived for 1:1 water-ethanol solutions (21)]. A value larger than 1 indicates that the reaction is more susceptible to substitution than benzoic acid in water ( ⫽ 1). For TP (and phenol) this can be rationalized by the closer proximity of the abstracted hydrogen to the aryl ring and the substituent than in the case of benzoic acid. Thermal desorption experiments by Meyers and Gellman (17) showed that the rate of Ullmann coupling between substituted iodobenzenes may also be understood in the framework of the Hammett equation even if the reaction proceeds under vacuum on a Cu(111) surface. Their value of 6.0 ⫾ 3.3 has a positive sign as reported here, which the authors interpret similarly as an electron-rich transition state. Although comparison of quantities obtained under very different conditions is problematic, we regard it as remarkable that in both experiments a value larger than 1 was found. This finding indicates that potential screening by the substrate electron density does not preclude or supersede strong interactions between substituents and reactive groups on aromatic rings. Consequently, linear free energy relationships resembling the Hammett equation may be applicable quite generally to surface (single-molecule) chemistry. Previously, STM-based molecular excitations had been analyzed by using either semiquantum mechanical models (5, 6) or sophisticated simulation of the electronic setup of the molecule–substrate system (22) and兾or the tunneling process (10, 23, 24). These approaches yield quantitative results for the 17922 兩 www.pnas.org兾cgi兾doi兾10.1073兾pnas.0406223101
Fig. 4. STM images and STM-based excitation of further thiophenols. (a–c) STM images of para-methyl-TP (⫺288 mV, 72 pA), meta-F-TP (⫺292 mV, 51 pA), and meta-Cl-TP (⫺273 mV, 58 pA) obtained at 16 K. The petals in b and c appear bent because of the asymmetric halo-substitution. (d) The linear dependence of the reaction rate on the tunneling current indicates a oneelectron process similar to Fig. 2.
excitation of simple species such as CO or ammonia. Systematic investigations of the STM-based reactivity of a class of substances (such as TPs in this study) have not been conducted so far. The former set of models do not lend themselves to the evaluation of data on polyatomic organic molecules because they lack the detail to permit predictions that depend on individual substituents. The latter class of models are computationally very involved, and modeling of the dynamics of substituted aromatic systems on metal surfaces tests the limits of the computer power currently available. At this stage, such modeling is not readily applicable for experimentalists. Consequently, simple semiempirical models and linear free energy relationships (such as the one presented here) may be the only available avenue toward the prediction of variations in STMbased reactivity between closely related organic substances. In the solution phase, effects of mobility and solvation often overshadow the intrinsic molecular properties that were at the center of Hammett’s original interest. The single-molecule measurements presented here offer a conceptually direct way to the original theoretical underpinning presented by Hammett.§ Historically, the Hammett equation paved the way for a systematic understanding of reaction rates and mechanisms in organic chemistry. This, in turn, led to the synthesis of a far wider range of reaction products. We hope that the development of simple guiding principles for single-molecule chemistry may similarly spark developments toward a better understanding of its mechanisms and rate dependencies and, ultimately, to its application in the construction of involved surface structures.
§Our
experiments are complementary to gas phase measurements in as much as the substrate provides for nonradiative deexcitation of reactants, facile momentum conservation, and stabilization of the abstracted hydrogen.
This work was supported by National Science Foundation Career Program Grant CHE-0132996 and the Petroleum Research Fund. Rao et al.
14. Jaffe´, H. H. (1953) Chem. Rev. 53, 191–261. 15. Star, A., Han, T. R., Gabriel, J. C. P., Bradley, K. & Gruner, G. (2003) Nano Lett. 3, 1421–1423. 16. Gellman, A. J., Buelow, M. T., Street, S. C. & Morton, T. H. (2000) J. Phys. Chem. A 104, 2476–2485. 17. Meyers, J. M. & Gellman, A. J. (1995) Surf. Sci. 337, 40–50. 18. Rao, B. V., Kwon, K. Y., Liu, A. W. & Bartels, L. (2003) J. Chem. Phys. 119, 10879–10884. 19. Jackson, G. J., Woodruff, D. P., Jones, R. G., Singh, N. K., Chan, A. S. Y., Cowie, B. C. C. & Formoso, V. (2000) Phys. Rev. Lett. 84, 119 – 122. 20. Carey, F. A. & Sundberg, R. J. (2000) Advanced Organic Chemistry (Kluwer Academic兾Plenum, New York), 4th Ed. 21. Schwarzenbach, G. & Rudin, E. (1939) Helv. Chim. Acta 22, 360–376. 22. Lorente, N., Persson, M., Lauhon, L. J. & Ho, W. (2001) Phys. Rev. Lett. 86, 2593–2596. 23. Sautet, P. (1997) Chem. Rev. 97, 1097–1116. 24. Pascual, J. I., Lorente, N., Song, Z., Conrad, H. & Rust, H. P. (2003) Nature 423, 525–528.
CHEMISTRY
1. Eigler, D. M. & Schweizer, E. K. (1990) Nature 344, 524–526. 2. Shen T. C., Wang, C., Abeln, G. C., Tucker, J. R., Lyding, J. W., Avouris, P. & Walkup, R. E. (1995) Science 268, 1590–1592. 3. Gimzewski, J. K. & Joachim, C. (1999) Science 283, 1683–1688. 4. Stranick, S. J., Kamna, M. M. & Weiss, P. S. (1994) Science 266, 99–102. 5. Bartels, L., Meyer, G., Rieder, K. H., Velic, D., Knoesel, E., Hotzel, A., Wolf, M. & Ertl, G. (1998) Phys. Rev. Lett. 80, 2004–2007. 6. Stipe, B. C., Rezaei, M. A., Ho, W., Gao, S., Persson, M. & Lundqvist, B. I. (1997) Phys. Rev. Lett. 78, 4410–4413. 7. Binnig, G., Rohrer, H., Gerber, C. & Weibel, E. (1983) Phys. Rev. Lett. 50, 120–123. 8. Hla, S. W., Bartels, L., Meyer, G. & Rieder, K. H. (2000) Phys. Rev. Lett. 85, 2777–2780. 9. Stipe, B. C., Rezaei, M. A. & Ho, W. (1998) Science 280, 1732–1735. 10. Komeda, T., Kim, Y., Kawai, M., Persson, B. N. J. & Ueba, H. (2002) Science 295, 2055–2058. 11. Hammett, P. L. (1937) J. Am. Chem. Soc. 59, 96–103. 12. Morrison, R. T. & Boyd, R. N. (1992) Organic Chemistry (Prentice–Hall, Englewood Cliffs, NJ), 6th Ed. 13. Johnson, C. D. (1973) The Hammett Equation (Cambridge Univ. Press, Cambridge, U.K.).
Rao et al.
PNAS 兩 December 28, 2004 兩 vol. 101 兩 no. 52 兩 17923
