Molecular Seesaw: Intricate Dynamics and ... - ACS Publications

Feb 23, 2017 - Molecular Seesaw: Intricate Dynamics and Versatile Chemistry of. Heteroaromatics on Metal Surfaces. Sergey N. Filimonov,*,†. Wei Liu,...
0 downloads 0 Views 6MB Size
Subscriber access provided by University of Newcastle, Australia

Letter

Molecular Seesaw: Intricate Dynamics and Versatile Chemistry of Heteroaromatics on Metal Surfaces Sergey N. Filimonov, Wei Liu, and Alexandre Tkatchenko J. Phys. Chem. Lett., Just Accepted Manuscript • DOI: 10.1021/acs.jpclett.7b00071 • Publication Date (Web): 23 Feb 2017 Downloaded from http://pubs.acs.org on March 3, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry Letters is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 11

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Molecular Seesaw: Intricate Dynamics and Versatile Chemistry of Heteroaromatics on Metal Surfaces Sergey N. Filimonov,∗,† Wei Liu,‡ and Alexandre Tkatchenko¶ †Department of Physics, Tomsk State University, 634050 Tomsk, Russia ‡Nano Structural Materials Center, School of Materials Science and Engineering, Nanjing University of Science and Technology, Nanjing, Jiangsu 210094, China ¶Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg E-mail: [email protected]

Abstract

Graphical TOC Entry

The design of novel elementary surface processes is important for applications in catalysis, single-molecule junctions, molecular sensors, switches and surface mounted molecular machines. Here we demonstrate by van der Waals inclusive density functional theory calculations that a small and relatively simple heteroaromatic compound s-triazine (C3 H3 N3 ) unexpectedly possesses five metastable states when adsorbed on the Pt(111) surface. This diversity of the adsorption states stems from an interplay between versatile molecule/surface chemical bonding and van der Waals interactions and from “softening” of the aromatic ring by nitrogen substitution, which makes folding of the aromatic ring energetically much less demanding as compared to benzene. The intricate seesaw-like surface dynamics and tunable electronic structure of s-triazine show promise for applications in molecular sensors and switches. The broad implications of our findings are demonstrated for triazine- and pyrimidinebased heteroaromatic compounds and other metal surfaces.

ACS Paragon Plus Environment

1

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 11

e.g. pyridine. 21,22 In the case of bonding solely through the delocalized π-electrons one would expect the aromatic ring to be parallel to the surface, while bonding via the N lone pair favors upright standing orientation of the adsorbed molecule. In a mixed case of bonding involving both delocalized π-electrons and lone electron pairs of N atoms one expects the adsorbed molecule to be tilted relatively to the surface. Remarkably, the mixed case renders different possible configurations at the same adsorption location, as will be presented below. To explore possible adsorption configurations of s-triazine on Pt(111) we performed structural optimization of the molecule at four highsymmetry adsorption positions (atop, bridge, hollow fcc and hollow hcp), starting from two different adsorption heights (2 and 3 ˚ A) and four tilt angles of the molecule to the surface: −30◦ (left), 0◦ (parallel), 30◦ (right), and 90◦ (normal). Our calculations were carried out using the van der Waals (vdW) inclusive DFT+vdWsurf method 23 (see Methods section for computational details), which enables accurate prediction of structure and stability of organic molecules on metal surfaces as has been demonstrated for a range of different adsorption systems. 24,25 The adsorption energy was calculated as Ead = Emol/metal − Emetal − Emol , where Emol/metal is the total energy of the complete adsorption system, Emetal and Emol are the energies of the clean Pt slab and isolated molecule, respectively. According to this definition more negative adsorption energies correspond to a stronger binding of the molecule to the substrate. Relaxed geometry, calculated adsorption energies and simulated constant-current scanning tunneling microscopy (STM) images of the revealed adsorption configurations are shown in Fig. 1. Similar to benzene 26–28 and to some other benzene derivatives, 29–31 s-triazine on Pt(111) preferentially occupies the bri30 position. However, in contrast to benzene, striazine does not lie parallel to the surface in its energetically most favorable adsorption configuration but shows considerable folding of the aromatic ring (Fig. 1(a)). This is originated

Chemical reactivity and physical properties of organic molecules deposited on inorganic surfaces depend crucially on the geometrical structure adopted by the molecule upon adsorption. This behavior is of paramount importance for selective organic catalysis 1–5 and also has important practical implications in material science and molecular electronics, where controlled transitions between different metastable configurations of adsorbed molecules are utilized to design novel molecular devices from simple molecular switches 6–9 to complex molecular machines mounted on surfaces. 10,11 The diversity of adsorption configurations usually stems from conformational changes of large and flexible molecules, 12–16 whereas the smallest aromatic compounds, like benzene and its simplest derivatives, are commonly considered to be rigid and to be unable to switch the adsorption state by a conformational change. In dense molecular arrays the change of the adsorption configurations of rigid molecules can occur by changing orientation of the molecule with respect to the surface. 17–19 For single adsorbed molecules such orientational transitions are much less anticipated due to the absence of intermolecular interactions, which might stabilize or destabilize different adsorption structures. However, in the present paper we show by using van der Waals-inclusive density functional theory calculations that reduced aromaticity of triazine and pyrimidine compounds combined with subtle interplay of the chemical bonding and dispersive forces can result in both the conformational and orientational transformations of the adsorbed molecules on reactive metal surfaces. The emergence of multiple states of adsorption at the single molecule level offers new perspectives for the smallest heteroaromatic compounds in catalysis and molecular electronics. We begin with the adsorption of the symmetric 1,3,5-triazine (C3 H3 N3 ), also called striazine, on Pt(111). When deposited on a metal surface s-triazine can form chemical bonds with the substrate atoms either via delocalized π electrons of the C-N ring or using lone pairs of electrons of nitrogen atoms, 19 similar to other N-substituted benzene derivatives,

ACS Paragon Plus Environment

2

Page 3 of 11

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Figure 1: Relaxed geometry and simulated constant-current STM images of five representative adsorption configurations of s-triazine on Pt(111): (a) bri30-left, (b) bri30-par, (c) bri30-right, (d) atop30 and (e) hcp30. STM images were calculated using the Tersoff-Hamann approach 20 at a positive bias voltage of 0.05 V. Red lines at the top and right sides of the STM images are the vertical and horizontal line scans across the molecule. from breaking the aromaticity of the C-N ring in the adsorbate. The folding helps to minimize the Pauli repulsion between the molecule and the substrate by lifting up two of three hydrogen atoms to a distance of 2.92 ˚ A from the surface. Thereby two of three nitrogen atoms and one carbon atom are brought very close to the surface. The adsorption height of 1.95 ˚ A of the lowermost N and C atoms suggests formation of chemical bonds with the metal. In the parallel arrangement s-triazine also adsorbs preferentially at the bridge site position. Relaxed geometry of the bri30-par state is characterized by an upward tilt of the C-H bonds and a slight nonplanarity of the aromatic ring with C atoms lying 0.02 ˚ A closer to the surface than N atoms (Fig. 1(b)). The adsorption height of 2.05–2.07 ˚ A in the bri30-par state is typical for chemisorption of small benzene derivatives on Pt(111). 29 The bri30-par state is energetically less favorable by 0.35 eV than the tilted bri30-left state. Another possible configuration of s-triazine, adsorbed at the bridge site, is bri30-right. Here the aromatic ring of the molecule is only slightly deformed, but the relaxed molecule is tilted by 24◦ relative to the substrate surface (Fig. 1(c)). Thereby one nitrogen atom and two carbon atoms reside close enough to the surface to form chemical bonds with the underlying Pt atoms.

The bri30-right state is energetically less favorable by 0.23 eV than the bri30-left state, but it is more favorable by 0.12 eV than the bri30-par state. For s-triazine oriented normal to the surface the preferential adsorption position is the atop site. In the optimal atop30 configuration the plane of the molecule is rotated by 30◦ with respect to the rows of surface atoms and the lowermost N atom of the molecule is located 2.19 ˚ A above Pt (Fig. 1(d)). With the adsorption energy of −1.11 eV, the atop30 state is the second energetically most favorable adsorption state, with bri30-left being the first. Apart from the four chemisorbed states our calculations reveal the existence of a physisorbed state of s-triazine on Pt(111). Among possible physisorbed configurations the most favorable one is hcp30 with the adsorption energy of −0.76 eV and the adsorption height of 2.89 ˚ A (Fig. 1(e)). However, as will be shown below, the energy barriers separating physisorbed and chemisorbed states of s-triazine on Pt(111) are very low, so that the molecule can easily transfer from the physisorbed state to one of the energetically more favorable chemisorbed states. All the revealed adsorption configurations of s-triazine on Pt(111) are strongly stabilized by vdW interactions. In particular, the vdW en-

ACS Paragon Plus Environment

3

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ergy accounts for more than 50% of the total binding energy in the two energetically most stable bri30-left and atop30 states (see Supporting Information for more details). Moreover, the contribution of vdW interactions determines relative stability of different adsorption configurations of s-triazine on Pt(111). Indeed, the atop30 state becomes equally stable with the bri30-left state when vdW interactions are not taken into account. This finding is in accord with our recent results revealed the critical role of vdW interactions in determining relative stability of the reactive intermediates of complex molecules on surfaces. 32 To check whether the choice of the method for vdW interactions could affect any of the revealed adsorption states of s-triazine on Pt(111) we have performed test calculations with the optimized optB88-vdW functional, 33 which represents a conceptually different approach to treatment of dispersive interactions. Using the PBE+vdWsurf relaxed geometry as an input we have observed only slight changes in the adsorption geometry when the optB88vdW method was utilized. That is all five metastable states of adsorption of s-triazine on Pt(111) were well reproduced when using the optB88-vdW functional. Selection of a particular adsorption configuration upon adsorption is determined by the orientation of the molecule with respect to the surface. Our calculations show that s-triazine can directly adsorb into the chemisorbed atop30 state when the molecule approaches the surface standing upright. In other cases the molecule can be temporary trapped in the physisorbed hcp30 state, from which it can transfer to one of the energetically more favorable chemisorbed states. The states which can be directly occupied by impinging molecules are marked by red arrows in Fig. 2(a). Possible transitions between the physisorbed and chemisorbed states and corresponding activation barriers are shown in Fig. 2(a). As can be seen the activation barrier for the transition from the physisorbed hcp30 state to the bri30-left state is substantially higher than the barriers for the transitions from the physisorbed hcp30 state to the chemisorbed bri30-right and

Page 4 of 11

Figure 2: (a) Adsorption mechanisms of striazine on Pt(111). Red arrows mark channels for direct (barrierless) adsorption from the gas phase, black arrows mark transitions between the physisorbed and chemisorbed states. (b) Potential energy diagram of the transitions between different chemisorbed states of s-triazine on Pt(111). Numbers in (a) and (b) indicate corresponding activation barriers in eV.

ACS Paragon Plus Environment

4

Page 5 of 11

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

bri30-par states. That is, the energetically most favorable bri30-left state is the most difficult one to reach from the physisorbed state, suggesting that filling up this state occurs via transient adsorption into less favorable but more accessible chemisorbed configurations. Bearing in mind the possibility of the direct (barrierless) adsorption into the chemisorbed atop30 state and a close to zero barrier for the transition from the physisorbed hcp30 state to the chemisorbed bri30-right state, we conclude that the most plausible scenario for s-triazine arriving at the Pt(111) surface is immediate chemisorption either into the tilted bri30-right state or into the upright standing atop30 state. We further consider the possibility of striazine trapped in a certain chemisorbed state to change its adsorption configuration and reach the energetically most favorable bri30left state. The potential energy landscape associated with the transitions between different chemisorbed states of s-triazine on Pt(111) is shown in Fig. 2(b). From the atop30 state s-triazine can switch to the bri30-right state via a 30◦ rotation around the N-Pt bond, followed by a 66◦ tilt of the aromatic ring toward the surface. From the bri30-right state the molecule can further switch to the bri30-par state and then from the bri30-par state it can finally reach the energetically most favorable bri30-left configuration. Since the molecule has to pass through a sequence of energetically less favorable states, the overall energy barrier for the atop30→bri30-left transition is 0.36 eV. The activation barrier for the reverse transition from the bri30-left to atop30 state is 0.50 eV. We also searched for possible diffusion pathways of s-triazine on Pt(111). We found that the barrier for surface diffusion of s-triazine in the physisorbed state (0.02 eV) is higher than the barrier for transition from the physisorbed hcp30 state to the chemisorbed bri30-right state. This excludes any long range lateral mobility of physisorbed molecules and rules out the possibility of diffusion of chemisorbed molecules via a transient physisorbed state. The lowest energy diffusion path for chemisorbed s-triazine involves transitions between two neighboring atop30 states. The energy barrier of 0.55 eV

for this transition is higher than any of the barriers in Fig. 2(b), therefore lateral mobility of chemisorbed s-triazine molecules should be suppressed on the timescale of switching between the atop30 and bri30-left configurations. Suppressed lateral mobility of s-triazine on Pt(111) and relatively high barriers for switching between the bri30-left and atop30 states should make the switching dynamics accessible for direct experimental investigations e.g. by STM. An elementary surface process can be tracked by STM if its rate is lower than the STM scan rate. In the following we identify three different temperature regimes of the adsorbate dynamics by comparing the typical scan rate of 1 frame/min with the transition rates of s-triazine on Pt(111). At low temperatures (T < 50 K) preferential adsorption of s-triazine into the upright standing atop30 configuration is expected. This state can be either directly occupied by impinging molecules, or filled up after short transient adsorption into less favorable bri30-par and bri30right states. Very low barriers for transitions from the bri30-par to bri30-right state and from the bri30-right to atop30 state make these transitions the dominant types of molecular moves at low substrate temperatures. In the middle temperature range (50–125 K), transitions from the atop30 to bri30-right, from bri30-right to bri30-par and from bri30-par to bri30-left states become possible. However, at middle and high temperatures it would be difficult to detect the bri30-par and bri30-right states with conventional STM due to a very short residence time of the molecules in these states. Therefore, the transition to the energetically most favorable bri30-left configuration should be seen as a spontaneous swap between the atop30 and bri30-left states. As a result of such transitions, the atop30 population will decay with time in favor of the bri30-left states. At elevated temperatures (T > 125 K) the net rate of the atop30 to bri30-left transitions becomes higher than the typical STM scan rate of 1 frame/min. Therefore, only the bri30-left state could be imaged with STM at high T . When the the bri30-left to atop30 transitions (via bri30-par and bri30-right states) become

ACS Paragon Plus Environment

5

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 11

Figure 3: Electronic properties of two most stable adsorption configurations of s-triazine on Pt(111): (a) bri30-left, and (b) atop30. Left panel: density of states of s-triazine on Pt(111) projected onto the HOMO and LUMO of a free s-triazine molecule. Middle panel: charge density difference isosurfaces of s-triazine/Pt(111). The isovalues are ±0.04 electron/˚ A3 . The regions of the electron accumulation (depletion) are colored blue (red). Right panel: adsorption induced charge density difference of s-triazine/Pt(111) integrated over the xy plane. activated, the STM image of the molecule might become blurred due to a seesaw-like motion of s-triazine switching back and forth between the bri30-left and atop30 states. Transitions between different adsorption states can be activated not only by thermal fluctuations, but also by some external stimuli. For instance, Borca et al. 34 and Schendel et al., 35 have demonstrated recently that transitions between chemisorbed and physisorbed states of individual organic molecules (anthradithiophene on Cu(111)) can be triggered by charge injection from a STM tip. One can expect that similar switching mechanism can be applied to s-triazine on Pt(111). Transitions between different adsorption states of s-triazine on Pt(111) cause changes in the electronic structure of the adsorbed molecule. The left column in Fig. 3 show the projected density of states (PDOS) of the adsorbed s-triazine onto the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) of the free molecule. 36 The HOMO and LUMO states of s-triazine broaden considerably and strongly overlap upon adsorption in the bri30-left state (Fig. 3(a), left). The former HOMO and LUMO

states of s-triazine in the bri30-left configuration are occupied by 1.7 and 0.98 electrons, respectively, indicating considerable charge redistribution between the metal and the molecule. This is confirmed by the electron density difference plots, shown in the middle and right columns of Fig. 3(a), where accumulation of electrons on the N-Pt and C-Pt bonds just below the N and C atoms and depletion of electrons above the metal surface are seen. Similar character of the charge redistribution along the N-Pt and C-Pt bonds suggests that in both cases bonding is primarily due to the delocalized π electrons of the molecule. The adsorption induced redistribution of charge between the molecule and the substrate, and deformation of the backbone of the adsorbed molecule, generate an interfacial dipole, which changes the work function of the metal surface. Interestingly, despite the apparent charge transfer from the substrate to the molecule, adsorption of s-triazine into the bri30-left state reduces the surface work function by 0.38 eV. We explain this by a significant polarization of the adsorbate (note considerable charge oscillations within the molecule), which overcompensates the charge transfer effect and contributes to the

ACS Paragon Plus Environment

6

Page 7 of 11

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

observed work function decrease. 37 In the vertical atop30 configuration the overlap of the former HOMO and LUMO states is negligible (Fig. 3(b)). Analysis of the projected orbitals occupation indicates that in the atop30 state there are 1.64 and 0.09 electrons in the former HOMO and former LUMO orbitals, respectively. The fact that the LUMO state is almost empty and a noticeable depletion of the HOMO state, suggest a net charge transfer from the molecule to the metal surface, in agreement with the integrated charge difference distribution showing considerable depletion of electrons below the molecule and accumulation of electrons close to the metal surface. The surface work function of Pt(111) is reduced by 0.8 eV due to adsorption of s-triazine in the vertical atop30 configuration. The key features responsible for appearance of multiple adsorption states and intricate dynamics of s-triazine on Pt(111) are (i) relatively high reactivity of the metal substrate, (ii) substantial role of dispersive interactions in the chemisorbed states, (iii) symmetry breaking between the molecule and the substrate and (iv) “softening” of the aromatic ring by nitrogen substitution. The (i) ensures chemical bonding of the molecule to the surface, (ii) and (iii) serve for enhanced stability of nonplanar adsorption configurations and (iv) makes deformation of the aromatic ring energetically much less demanding as compared to benzene. Evidently, neither of the aforementioned features is unique to s-triazine on Pt(111), therefore one would expect similar behavior for a wide range of different adsorption systems. To confirm this we have performed structural relaxation of s-triazine on different metal surfaces and checked out a number of other heterocyclic compounds on Pt(111). Multiple states of adsorption connected by seesaw-like transitions have been found for s-triazine on Pd(111) and Rh(111) and for a number of triazine and pyrimidine derivatives on Pt(111) (see Fig. 4 and Supporting Information for more details). The occurrence of multiple state of adsorption offers new perspectives for the smallest heteroaromatic compounds in catalysis and molecular electronics.

Figure 4: Triazine and pyrimidine compounds with multiple adsorption states on Pt(111). In summary, the structure and electronic properties of s-triazine adsorbed on Pt(111) were investigated with density functional theory calculations. It is shown that being trapped at a certain adsorption position s-triazine on Pt(111) can adopt several distinct adsorption configurations, which differ considerably by the adsorption geometry, bonding configuration and electronic properties of the adsorbed molecule. Noticeably, reversible transitions between different adsorption configurations may occur without changing the location of the molecule on the surface, which makes s-triazine a strong candidate for molecular switch applications.

Methods Density-functional theory calculations with semilocal Perdew-Burke-Ernzerhof (PBE) functional 38 were performed using the FHI-aims all-electron code. 39 The dispersive interactions were treated with the DFT+vdWsurf method, 23 which combines the DFT+vdW method 40 for intermolecular interactions with the Lifshitz-Zaremba-Kohn theory 41,42 for the nonlocal Coulomb screening within the inorganic bulk. The Pt(111) surface was represented by a periodic slab with a (4 × 4) unit cell containing six layers of Pt and 29 ˚ A of vacuum. The molecule was placed on the top surface of the slab. Dipole correction was applied to compensate an artificial dipole field due to

ACS Paragon Plus Environment

7

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 11

Supporting Information Available

non-equivalent slab surfaces. 43 For the geometry optimization, the molecule and the two uppermost metal layers were allowed to relax, while atoms of four bottommost Pt layers were constrained at their bulk positions employing the theoretical PBE+vdWsurf lattice constant of 3.979 ˚ A . 28 The structure was optimized until the maximum residual force component per atom was converged to 10−4 eV. The final single-point calculations of the relaxed structures were carried out with the “tight” settings including the “tier2” standard basis set for light elements and “tier1” for heavy metal atoms. Relativistic effects were treated using the scaled zero-order regular approximation. 44 Convergence criteria of 10−5 electrons per unit volume for the electron density and 10−4 eV for the total energy of the system were applied. The reciprocal space integration was performed with a MonkhorstPack grid of (6 × 6 × 1) kpoints. We also performed test calculations with the optB88-vdW functional, 33 to confirm that the presence of the revealed adsorption states of s-triazine on Pt(111) is not affected by the method for dispersive interactions. The optB88-vdW calculations were performed with the VASP code 45 using PBE-based projector augmented wave potentials 46 with a cut-off energy of 400 eV and a 2 × 2 × 1 MonkhorstPack grid of k-points. The activation energies for lateral diffusion of the adsorbed molecule and for transitions between different adsorption states were calculated using the string method 47 combined with the climbing image technique, 48 as implemented in the aimsChain code from the standard FHI-aims distribution. The rate constants of the surface processes were evaluated using the harmonic transition state theory 49 and the DFT-calculated activation energies and normal vibrational modes for the metastable and transition states. The structural models, simulated STM images and charge density maps are visualized with the VESTA software 50

The following files are available free of charge. Transition rates of s-triazine on Pt(111); adsorption energies of s-triazine on Pt(111), Pd(111) and Rh(111); adsorption energies of triazine and pyrimidine derivatives on Pt(111).

Author information Corresponding Author: [email protected] (S.N.F.) Notes: The authors declare no competing financial interests Acknowledgement S.N.F. gratefully acknowledges financial support from the Russian Science Foundation through the grant 14-1200813 and computational resources provided by the Supercomputing Center of Tomsk State University. W.L. acknowledges support from the National Natural Science Foundation of China (no. 21403113), and the Natural Science Foundation for Distinguished Young Scholars of Jiangsu Province (no. BK20150035).

References (1) Zaera, F. Outstanding Mechanistic Questions in Heterogeneous Catalysis. J. Phys. Chem. B 2002, 106, 4043–4052. (2) Griffin, M. B.; Pang, S. H.; Medlin, J. W. Surface Chemistry of 2-Iodoethanol on Pd(111): Orientation of Surface-Bound Alcohol Controls Selectivity. J. Phys. Chem. C 2012, 116, 4201–4208. (3) Brandt, K.; Chiu, M. E.; Watson, D. J.; Tikhov, M. S.; Lambert, R. M. Adsorption Geometry Determines Catalytic Selectivity in Highly Chemoselective Hydrogenation of Crotonaldehyde on Ag(111). J. Phys. Chem. C 2012, 116, 4605–4611. (4) Wang, S.; Vorotnikov, V.; Vlachos, D. G. Coverage-Induced Conformational Effects

ACS Paragon Plus Environment

8

Page 9 of 11

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

on Activity and Selectivity: Hydrogenation and Decarbonylation of Furfural on Pd(111). ACS Catal. 2015, 5, 104–112.

Induced by Scanning Tunneling Microscopy Manipulation: A Route to Molecular Switching. Phys. Rev. Lett. 2001, 86, 672–675.

(5) Dostert, K.-H.; O’Brien, C. P.; Liu, W.; Riedel, W.; Savara, A.; Tkatchenko, A.; Schauermann, S.; Freund, H.-J. Adsorption of Isophorone and TrimethylCyclohexanone on Pd(111): A Combination of Infrared Reflection Absorption Spectroscopy and Density Functional Theory Studies. Surf. Sci. 2016, 650, 149–160.

(14) Iancu, V.; Deshpande, A.; Hla, S.-W. Manipulating Kondo Temperature via Single Molecule Switching. Nano Lett. 2006, 6, 820–823. (15) Choi, B.-Y.; Kahng, S.-J.; Kim, S.; Kim, H.; Kim, H. W.; Song, Y. J.; Ihm, J.; Kuk, Y. Conformational Molecular Switch of the Azobenzene Molecule: A Scanning Tunneling Microscopy Study. Phys. Rev. Lett. 2006, 96, 156106.

(6) Browne, W. R.; Feringa, B. L. Light Switching of Molecules on Surfaces. Annu. Rev. Phys. Chem. 2009, 60, 407–428.

(16) Weigelt, S.; Busse, C.; Petersen, L.; Rauls, E.; Hammer, B.; Gothelf, K. V.; Besenbacher, F.; Linderoth, T. R. Chiral Switching by Spontaneous Conformational Change in Adsorbed Organic Molecules. Nat. Mater. 2006, 5, 112–117.

(7) Morgenstern, K. Switching Individual Molecules by Light and Electrons: From Isomerisation to Chirality Flip. Prog. Surf. Sci. 2011, 86, 115–161. (8) Pathem, B. K.; Claridge, S. A.; Zheng, Y. B.; Weiss, P. S. Molecular Switches and Motors on Surfaces. Annu. Rev. Phys. Chem. 2013, 64, 605.

(17) Demuth, J.; Christmann, K.; Sanda, P. The Vibrations and Structure of Pyridine Chemisorbed on Ag(111): the Occurrence of a Compressional Phase Transformation. Chem. Phys. Lett. 1980, 76, 201–206.

(9) Zhang, J. L.; Zhong, J. Q.; Lin, J. D.; Hu, W. P.; Wu, K.; Xu, G. Q.; Wee, A. T. S.; Chen, W. Towards Single Molecule Switches. Chem. Soc. Rev. 2015, 44, 2998–3022.

(18) Bader, M.; Haase, J.; Frank, K. H.; Puschmann, A.; Otto, A. Orientational Phase Transition in the System Pyridine/Ag(111): A Near-Edge X-RayAbsorption Fine-Structure Study. Phys. Rev. Lett. 1986, 56, 1921–1924.

(10) Browne, W. R.; Feringa, B. L. Making Molecular Machines Work. Nat. Nanotechnol. 2006, 1, 25–35. (11) Balzani, V.; Credi, A.; Venturi, M. Molecular Machines Working on Surfaces and at Interfaces. Chemphyschem 2008, 9, 202– 220.

(19) Mart´ınez-Galera, A. J.; de la Torre, B.; Ugeda, M. M.; G´omez-Rodr´ıguez, J. M. Adsorption and Growth of 1,3,5-Triazine on Cu(111) at Low Temperature under Ultrahigh Vacuum Conditions. J. Phys. Chem. C 2012, 116, 9568–9574.

(12) Collier, C. P.; Mattersteig, G.; Wong, E. W.; Luo, Y.; Beverly, K.; Sampaio, J.; Raymo, F. M.; Stoddart, J. F.; Heath, J. R. A [2]Catenane-Based Solid State Electronically Reconfigurable Switch. Science 2000, 289, 1172–1175.

(20) Tersoff, J.; Hamann, D. R. Theory and Application for The Scanning Tunneling Microscope. Phys. Rev. Lett. 1983, 50, 1998–2001.

(13) Moresco, F.; Meyer, G.; Rieder, K.-H.; Tang, H.; Gourdon, A.; Joachim, C. Conformational Changes of Single Molecules

(21) Bili´c, A.; Reimers, J. R.; Hush, N. S. Adsorption of Pyridine on the Gold(111) Surface: Implications for “Alligator Clips” for

ACS Paragon Plus Environment

9

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Molecular Wires. J. Phys. Chem. B 2002, 106, 6740–6747.

Page 10 of 11

(29) Liu, W.; Filimonov, S.; Carrasco, J.; Tkatchenko, A. Molecular Switches from Benzene Derivatives Adsorbed on Metal Surfaces. Nat. Commun. 2013, 4, 2569.

(22) Mollenhauer, D.; Gaston, N.; Voloshina, E.; Paulus, B. Interaction of Pyridine Derivatives with a Gold (111) Surface as a Model for Adsorption to Large Nanoparticles. J. Phys. Chem. C 2013, 117, 4470–4479.

(30) Pek¨oz, R.; Johnston, K.; Donadio, D. Tuning the Adsorption of Aromatic Molecules on Platinum via Halogenation. J. Phys. Chem. C 2014, 118, 6235–6241. (31) Johnston, K.; Pekoz, R.; Donadio, D. Adsorption of Polyiodobenzene Molecules on the Pt(111) Surface using van der Waals Density Functional Theory. Surf. Sci. 2016, 644, 113–121.

(23) Ruiz, V. G.; Liu, W.; Zojer, E.; Scheffler, M.; Tkatchenko, A. Screened van der Waals Interactions: The Key to Accurately Model The Structure of Hybrid Inorganic/Organic Systems. Phys. Rev. Lett. 2012, 108, 146103.

(32) Rodriguez-Reyes, J. C. F.; Siler, C. G. F.; Liu, W.; Tkatchenko, A.; Friend, C. M.; Madix, R. J. Van der Waals Interactions Determine Selectivity in Catalysis by Metallic Gold. J. Am. Chem. Soc. 2014, 136, 13333–13340.

(24) Liu, W.; Tkatchenko, A.; Scheffler, M. Modeling Adsorption and Reactions of Organic Molecules at Metal Surfaces. Acc. Chem. Res. 2014, 47, 3369–3377. (25) Maurer, R. J.; Ruiz, V. G.; CamarilloCisneros, J.; Liu, W.; Ferri, N.; Reuter, K.; Tkatchenko, A. Adsorption Structures and Energetics of Molecules on Metal Surfaces: Bridging Experiment and Theory. Prog. Surf. Sci. 2016, 91, 72–100.

(33) Klimeˇs, J.; Bowler, D. R.; Michaelides, A. Van der Waals Density Functionals Applied to Solids. Phys. Rev. B 2011, 83, 195131. (34) Borca, B.; Schendel, V.; Ptuya, R.; Pentegov, I.; Michnowicz, T.; Kraft, U.; Klauk, H.; Arnau, A.; Wahl, P.; Schlickum, U. et al. Bipolar Conductance Switching of Single Anthradithiophene Molecules. ACS Nano 2015, 9, 12506–12512.

(26) Saeys, M.; Reyniers, M. F.; Marin, G. B.; Neurock, M. Density Functional Study of Benzene Adsorption on Pt(111). J. Phys. Chem. B 2002, 106, 7489–7498. (27) Morin, C.; Simon, D.; Sautet, P. Chemisorption of Benzene on Pt(111), Pd(111), and Rh(111) Metal Surfaces : A Structural and Vibrational Comparison from First Principles. J. Phys. Chem. B 2004, 108, 5653–5665.

(35) Schendel, V.; Borca, B.; Pentegov, I.; Michnowicz, T.; Kraft, U.; Klauk, H.; Wahl, P.; Schlickum, U.; Kern, K. Remotely Controlled Isomer Selective Molecular Switching. Nano Lett. 2016, 16, 93– 97.

(28) Liu, W.; Ruiz, V. G.; Zhang, G.X.; Santra, B.; Ren, X.; Scheffler, M.; Tkatchenko, A. Structure and Energetics of Benzene Adsorbed on Transition-Metal Surfaces: Density-Functional Theory with van der Waals Interactions Including Collective Substrate Response. New J. Phys. 2013, 15, 053046.

(36) Rangger, G. M.; Romaner, L.; Heimel, G.; Zojer, E. Understanding the Properties of Interfaces Between Organic SelfAssembled Monolayers and Noble Metals — a Theoretical Perspective. Surf. Int. Anal. 2008, 40, 371–378. (37) Gossenberger, F.; Roman, T.; ForsterTonigold, K.; Groß, A. Change of the

ACS Paragon Plus Environment

10

Page 11 of 11

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Work Function of Platinum Electrodes Induced by Halide Adsorption. Beilstein J. Nanotechnol. 2014, 5, 152–161.

(47) E, W.; Ren, W.; Vanden-Eijnden, E. Simplified and Improved String Method for Computing the Minimum Energy Paths in Barrier-Crossing Events. J. Chem. Phys. 2007, 126, 164103.

(38) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868.

(48) Henkelman, G.; Uberuaga, B. P.; J´onsson, H. A Climbing Image Nudged Elastic Band Method for Finding Saddle Points and Minimum Energy Paths. J. Chem. Phys. 2000, 113, 9901–9904.

(39) Blum, V.; Gehrke, R.; Hanke, F.; Havu, P.; Havu, V.; Ren, X.; Reuter, K.; Scheffler, M. Ab Initio Molecular Simulations with Numeric Atom-Centered Orbitals. Comp. Phys. Comm. 2009, 180, 2175–2196.

´ (49) Campbell, C. T.; Arnad´ ottir, L.; Sellers, J. R. V. Kinetic Prefactors of Reactions on Solid Surfaces. Z. Phys. Chem. 2013, 227, 1435–1454.

(40) Tkatchenko, A.; Scheffler, M. Accurate Molecular Van Der Waals Interactions from Ground-State Electron Density and Free-Atom Reference Data. Phys. Rev. Lett. 2009, 102, 073005.

(50) Momma, K.; Izumi, F. VESTA 3 for Three-Dimensional Visualization of Crystal, Volumetric and Morphology Data. J. Appl. Crystallogr. 2011, 44, 1272–1276.

(41) Zaremba, E.; Kohn, W. Van der Waals Interaction Between an Atom and a Solid Surface. Phys. Rev. B 1976, 13, 2270– 2285. (42) Lifshitz, E. M. The Theory of Molecular Attractive Forces Between Solids. Sov. Phys. JETP 1956, 2, 73–83. (43) Neugebauer, J.; Scheffler, M. AdsorbateSubstrate and Adsorbate-Adsorbate Interactions of Na and K Adlayers on Al(111). Phys. Rev. B 1992, 46, 16067–16080. (44) van Lenthe, E.; Baerends, E. J.; Snijders, J. G. Relativistic Total Energy Using Regular Approximations. J. Chem. Phys. 1994, 101, 9783. (45) Kresse, G.; Furthm¨ uller, J. Efficient Iterative Schemes for Ab Initio Total-Energy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B 1996, 54, 11169–11186. (46) Kresse, G.; Joubert, D. From Ultrasoft Pseudopotentials to the Projector Augmented-Wave Method. Phys. Rev. B 1999, 59, 1758–1775.

ACS Paragon Plus Environment

11